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INTRODUCTION 


These  volumes  contain  the  papers  presented  at  the  Aerodynamics 
Range  Symposium,  sponsored  by  the  Ballistic  Research  Laboratories . 
The  meeting  was  held  at  Aberdeen  Proving  Ground  on  2k  and  25  January 

1957- 


The  purposes  of  the  symposium  were  to  introduce  the  free  flight 
range  technique  to  a  wider  audience  and  to  discuss  applications  of 
the  technique  to  problems  beyond  the  shell  flight  mechanics  area. 

Until  recently,  it  has  been  generally  felt  that  free  flight 
ballistic  ranges  are  limited  principally  to  measurement  of  the 
small  amplitude  motion  of  shell.  The  papers  show  much  broader 
applications  of  ranges.  These  include  investigations  of  :  nonlinear 
dynamics  of  shell,  motion  of  aircraft  models  and  of  burning  rockets, 
basic  problems  in  fluid  mechanics,  and  the  design  of  hypersonic  re¬ 
entry  vehicles . 

It  now  appears  that  the  applicability  of  free  flight  spark 
photography  ranges  and  of  the  older  yaw  card  technique  is  limited 
only  by  the  Imagination  of  the  user.  Among  the  still  incompletely 
explored  possibilities  are  the  application  of  telemetry  to  heat 
transfer  measurements  and  to  measurement  of  fluid  properties  at 
very  high  Mach  Numbers . 

The  success  of  the  meeting  suggests*  that  a  similar  meeting 
should  he  held  within  the  next  two  years  to  help  exploit  the  rapid 
expansion  in  free  flight  range  technology. 


C. lL.  POOR 
General  Chairman 
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LIST  OF  CONFEREES 


In  addition  to  representatives  of  the  Ballistic  Research 
Laboratories,  the  following  attended  the  Aerodynamics  Range 
Symposium. 


ADALIA  LIMITED 
Montreal,  P.  Q. 

Reid,  J.  B. 

AEROFHYSICS  DEVELOPMENT  CORP. 
Pacific  Palisades,  California 
Bitondo,  D. 

AIRCRAFT  ARMAMENT,  INC. 
Cockeysville ,  Maryland 
La  Costa,  R.  S. 

AIR  FORCE  ARMAMENT  CENTER 
Eglin  Air  Force  Base,  Florida 
Haltom,  C .  M. 

Jacobs,  R. 

Stevens,  J.  E. 

AIR  FORCE  MISSILE  TEST  CENTER 
Patrick  Air  Force  Base,  Florida 
Burgess,  J.  A. 

Mason,  H.  P. 

APPLIED  PHYSICS  LABORATORY 
Silver  Spring,  Maryland 
Cramer,  R.  H. 

Cronvich,  L.  L. 

Wallskog,  H.  A. 

ARMY  BALLISTIC  MISSILE  AGENCY 
Huntsville,  Alabama 
Reed,  T.  G. 

ARO,  INC. 

Tullahoma,  Tennessee 
Kingery,  M.  K. 

Zazzi,  A.  J. 

AVCO  MANUFACTURING  COMPANY 
Lawrence,  Massachusetts 
Munson,  T,  R. 

Stephenson,  W.  B. 


BELL  AIRCRAFT  CORPORATION 
Buffalo,  New  York 
Meullen,  N.  F. 

Pearce,  K. 

Whalen,  R.  J. 

BOEING  AIRPLANE  COMPANY 
Seattle,  Washington 
Bateman,  R.  E. 

Martin,  D. 

Wichita,  Kansas 

Wallace,  R.  E. 

BRITISH  JOINT  SERVICES  MISSION 
Washington,  D.  C. 

Clarke,  A.  E. 

Walker,  N.  K. 

BUDD  COMPANY 

Philadelphia,  Pennsylvania 
Bomberger,  E.  E. 

Kessler,  E.  L. 

Simpson,  E.  L. 

Zettlemoyer,  E.  A. 

BUREAU  OF  ORDNANCE,  NAVY  DEPARTMENT 
Washington,  D.  C. 

Baker,  W.  K. 

Donoghue,  F.  D. 

Griffin,  T.  F. 

Nicolaides ,  J.  D . 

CANADAIR 
Montreal,  P.  0. 

Jackman,  D.  A. 

Luckert,  H.  J. 

CANADIAN  WESTINGHOUSE  COMPANY 
Hamilton,  Ontario 
Leavitt ,  J.  N . 


T 


CANADIAN  ARMAMENT  RESEARCH  AND 
DEVELOPMENT  ESTABLISHMENT 
Quebec,  P.  Q. 

Bull,  G.  V. 

Gre  envood ,  E .  W. 

Tidy,  G.  H. 

Waldock,  D.  A.  G. 

CHAMBERLAIN  CORPORATION 
Waterloo,  Iowa 
Caponi,  E.  R. 

CHICAGO  MIDWAY  LABORATORIES 
Chicago,  Illinois 
Mails,  M.  F. 

COMPUTING  DEVICES  OF  CANADA  LTD. 
Ottawa,  Ontario 
Howland,  J.  L, 

Jeffery,  C.  B. 

Smith,  J.  E. 

CONVAIR 

Dangerfield,  Texas 
Goin,  K.  L. 

San  Diego,  California 
Frick,  C.  W. 

Katz,  E. 

Romig,  M.  F. 

Starr,  S.  V. 

CORNELL  AERONAUTICAL  LABORATORY, 
Buffalo,  New  York 
Russo,  A.  L. 

Whiting,  A.  A. 

DEFENCE  RESEARCH  BOARD 
Ottawa,  Ontario 
Oatway,  J.  L. 

DeHAVILLAND  AIRCRAFT  OF  CANADA, 
Toronto,  Ontario 
Pounder,  D.  W. 

Warren,  H.  R. 


DIAMOND  ORDNANCE  FUZE  LABORATORIES 
Washington,  D.  C. 

Bowles,  R. 

Brody,  P. 

Cur check,  H. 

Dellasanta,  0. 

Harris,  F. 

Keto,  F. 

Van  der  Linden,  R. 

DIRECTORATE  OF  ARMAMENT  DEVELOPMENT 
Ottawa,  Ontario 
Snarr,  W.  B. 

DOUGLAS  AIRCRAFT  COMPANY,  INC. 

El  Segundo,  California 
(Newton,  F.  C. 

Santa  Monica,  California 
Cohen,  W.  S. 

Klein,  H. 

FIRESTONE  TIRE  AND  RUBBER  COMPANY 
Akron,  Ohio 

LucaB,  V.  E. 

FRANKFORD  ARSENAL 
Philadelphia,  Pennsylvania 
Brady,  J. 

Bushey,  B. 

Ciancioal,  A. 

Dickey,  C. 

INC,  Hall  or  an,  J. 

Kymer,  J. 

LipiaBki,  H. 

Otis,  W. 

Steinberg,  J. 

Tanik,  L. 

Yannuzzi,  A. 


LTD. 

GENERAL  ELECTRIC 
Philadelphia,  Pennsylvania 
Kebely,  V. 

Peck,  R.  F. 

Powers ,  J. 

Smith,  A.  M. 

Warren,  W.  R. 

Yoler,  Y.  A. 
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GENERAL  MILLS,  INC. 

Minneapolis,  Minnesota 
Hakoraaki ,  R .  I . ' 

GLENN  L.  MARTIN  COMPANY 
Baltimore,  Maryland 
Bidwell,  J.  M. 

Coon,  M.  L. 

Jackson,  L.  L. 

Reisert,  T. 

GRUMMAN  AIRCRAFT  ENGINEERING  CORP, 
Bethpage,  New  York 
Gustafson,  R.  L. 

Munier,  A.  E. 

Scheuing,  R.  A. 

HUGHES  AIRCRAFT  COMPANY 
Culver  City,  California 
Kalman,  I. 

Phillips,  W.  L. 

JET  PROPULSION  LABORATORY 
Pasadena,  California 
Elmer,  M. 

Schunnier,  H.  R. 

Wegener,  P.  P  . 

LOCKHEED  AIRCRAFT  CORP. 

San  Jose,  California 
Marsh,  B.  W. 

Tucker,  M. 

Wilson,  L.  H. 

Van  Nuys,  California 
Bershaded,  E. 

Cannon,  E.  T. 

Swanson,  R.  S. 

McDonnell  aircraft  corp. 

St.  Louis,  Missouri 
Rohtert,  R.  E. 

MASSACHUSETTS  INSTITUTE 
OF  TECHNOLOGY 
Cambridge,  Massachusetts 
Burgin,  F.  H. 

Wilkie,  I.  E. 


NATIONAL  ADVISORY  COMMITTEE  FOR 
AERONAUTICS 

Headquarters,  Washington,  D.  C. 
Abbott,  I.  H. 

May,  R.  E. 

Ames  Laboratory 
Moffett  Field,  California 
Allen,  H.  J. 

Canning,  T. 

Neice,  L. 

Langley  Aeronautical  Laboratory 
Langley  Field,  Virginia 
Bird,  J.  D. 

Hopko,  R. 

Huber,  P. 

I'  Sabol,  A . 

Lewis  Flight  Propulsion  Laboratory 
Cleveland,  Ohio 
Reshotko,  E . 

NATIONAL  RESEARCH  COUNCIL  OF  CANADA 
Ottawa,  Ontario 

Lukasiewicz,  J. 

Templin,  R.  J. 

NAVAL  AIR  MISSILE  TEST  CENTER 
Point  Mugu,  California 
Fitzgerald,  M. 

NAVAL  ORDNANCE  LABORATORY 
Silver  Spring,  Maryland 
Auld,  C.  D. 

Brady,  J.  J. 

Carter,  H. 

Chamberlin,  A. 

Crogan,  L.  E. 

Dawson,  V.  C. 

Diggins,  J.  L. 

Fedenia,  J.  N. 

Gallagher,  J.  J. 

Greenwald,  A. 

Hall,  D. 

Jusino,  J.  B. 

Kurzweg,  H.  H. 

Lightfoot,  G.  R. 

Long,  J.  E. 

May,  A. 

Roberts,  J. 

Seigal,  A..  E. 
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NAVAL  ORDNANCE  LABORATORY  (Cant**) 
Sheppard,  B.  S. 

Slawsky,  Z.  I. 

Thurston,  P,  A. 

White,  C.  E. 

Winkler,  E. 

Witt,  W.  R. 

NAVAL  ORDNANCE  TEST  STATION 
China  Lake,  California 
Allan,  W.  H. 

Dunn,  E.  L. 

Haseltine,  W.  R. 

Mayfield,  E.  B. 

Nevkirk,  H.  L. 

NAVAL  PROVING  GROUND 
Dahlgren,  Virginia 
Anderle,  R..  J. 

Brown,  D.  R. 

Cohen,  C.  J. 

Frick,  C.  H. 

Hubbard,  E.  C. 

Jones,  A.  L. 

Mitchell,  J. 

Mulligan,  J.  E. 

NAVAL  RESEARCH  LABORATORY 
Washington,  C. 

Atkins,  W.  W. 

Fuller,  R.  H. 

NORTH  AMERICAN  AVIATION,  INC. 
Columbus,  Ohio 
Simon,  W.  E. 

Downey,  California 
Briggs,  E. 

Dev,  J.  K. 

Oliver,  R.  B. 

OFFICE  ASSISTANT  SECRETARY 
OF  DEFENCE  (R&D) 

Washington,  D.  C. 

Nest ingen,  I. 

OFFICE  OF  NAVAL  RESEARCH 
Washington,  D.  C. 

Sherman,  F.  S. 


PICATINNY  ARSENAL 
Dover,  Nev  Jersey 
Carson,  W. 

Costa,  D. 

Loeb,  A.  A. 

Pallingston,  A. 

Spec tor,  J. 

Trevorrov,  D. 

RAMO-WOOLDRIDGE  CORP. 

Los  Angeles,  California 
Ambrosio,  A. 

Hayes,  W.  D. 

REDSTONE  ARSENAL 
Huntsville,  Alabama 
Murphree,  W.  D. 

Newby,  D.  H. 

Northrop,  C.  L. 

Shapiro,  N.  M. 

REPUBLIC  AVIATION  CORP. 

Farmingdale,  Nev  York 
Mcllroy,  W. 

Rennemann,  C . 

Sanator,  R. 

MJLneola,  Nev  York 
Kennedy,  R.  M. 

A.  V.  ROE  CANADA  LTD. 

Malton,  Ontario 

Chamberlin,  J.  A. 

Taylor,  W. 

ROSEMOUNT  AERONAUTICAL  LABORATORIES 
Rosemount,  Minnesota 
Deleo,  R.  V. 

Herman,  R. 

Moynihan,  F.  A. 

£ANDIA  CORPORATION 
Albuquerque,  Nev  Mexico 
Vaughn,  H.  R. 

UNITED  AIRCRAFT  CORPORATION 
East  Hartford,  Connecticut 
Taylor,  W. 
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UNIVERSITY  OF  SOUTHERN  CALIFORNIA 
Los  Angeles,  California 
Dailey,  C.  L. 

WATERTOWN  ARSENAL 
Watertown,  Massachusetts 
Muldoon,  R. 

WRIGHT  AIR  DEVEL0R4ENT  CENTER 
Dayton,  Ohio 

Huber,  F.  J. 

Liu,  Tung-Sheng 
Shorr,  M. 
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AGENDA* 


Thursday,  2k  January  1957 


Registration 

8:15 

-  9:00  am 

Welcome 

9:00 

-  9:15  am 

Col.  C.  L.  Register 

Director,  Ballistic  Research  Laboratories 

Remarks  by  General  Chairman 

9:15 

-  9:50  am 

Mr.  C.  L.  Poor 

Chief,  Exterior  Ballistics  Laboratory 

Session  I  9:50  -  11:10  am 

Chairman,  Dr.  H.  E.  Kurzweg 
Associate  Technical  Director 
Naval  Ordnance  Laboratory 

Simulation  of  the  Atmospheric  Entry  of  Ballistic 
Missiles  (CONFIDENTIAL) 

S.  E.  Nelce  and  J.  A.  Carson 
Ames  Aeronautical  Laboratory,  NACA 

Some  Problems  Associated  with  the  Determination, 
from  Range  Firings,  of  Dynamic  Stability  of 
Ballistic  Missile  Re-Entry  Shapes  (CONFIDENTIAL) 

L.  C,  MacAlllster 

Ballistic  Research  Laboratories 

Lunch  At  Transonic  Range  and  Officers  Club 

Session  II  1:00  -  VAO  pm 

Chairman,  Mr.  H.  Julian  Allen 
Chief,  High-Speed  Research  Division 
Ames  Aeronautical  Laboratory,  NACA 

A  Correlation  of  Free -Flight  Transition  Measurements 
On  Various  Blunt  Nose  Shapes  by  Use  of  the  Momentum- 
Thickness  Reynolds  Number 

W.  R.  Witt,  Jr.  and  J.  Eersh 
Naval  Ordnance  Laboratory 

Advances  in  the  Dynamic  Analysis  of  Range  Data 
C.  5.  Murphy 

Ballistic  Research  Laboratories 


* 

Papers  were  presented  by  those  whose  names  are  underlined. 
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Survey,  Calibration,  and  Reduction  Techniques 
Used  at  the  Thompson  Aeroballistics  Laboratory 
W.  H.  Allan  and  E,  L,  Dunn 
Naval  Ordnance  Test  Station 

Wake  Visualization  Studies  in  the  Aeroballistics 
Range 

p.  V.  Bnl 1 ,  and  C.  B.  Jeffery 

Canadian  Armament  Research  and  Development 

Establishment 

The  Controlled-Temperature-Pressure  Range 
E,  D.  Bennett 

Ballistic  Research  Laboratories 

Tour  5:00  -  6:00  pm 

Aerodynamics  Range 

CTF  Range  and  Wind  Tunnels 

Social  Hour  and  Dinner 


Friday,  25  January  1957 


Session  III  8:50  -  11:20  am 

Chairman,  Brigadier ,iD.  A.  G.  Waldock 
Chief  Superintendent 

Canadian  Armament  Research  and  Development 
Establishment 

Two  Aeroballistic  Range  Topics: 

(l)  Mass  Asymmetry*]  (2)  Dynamic  Stability 
J.  D.  Nlcolaldes 
-Bureau  of  Ordnance 

Design  and  Initial  Tests  of  the  N0L  Shock  Gun  (CONFIDENTIAL) 

V.  C.  D.  Davson 

Naval  Ordnance  Laboratory 

An  Application  of  Aeroballistics  Range  Techniques 
to  Models  vith  DigonaJ  Rotational  Symmetry 
G.  H,  Tidy  and  M.  E.  Thomas 
Oanad  1  an  Armp-rnsnt  Research  and  Development 
Establishment 

Sabots  Used  at  the  Thompson  Aeroballistics  Laboratory 

W.  H.  Allan 

Naval  Ordnance  Test  Station 
Lunch  At  Transonic  Range  and  Officers  Club 


With  J.  E.  Long,  U.  S,  Naval  Ordnance  Laboratory  and  G.  Parrish, 
Bureau  of  Ordnance, 


Xk 


Session  IV  1:00  -  3:20  pm 

Chairman,  Dr.  W.  R.  Haseltine 
Head,  Ballistics  Division 
Naval  Ordnance  Test  Station 

Research  Investigations  in  the  Ames  Supersonic  Free- 
Flight  Facilities  (CONFIDENTIAL) 

T.  N»  Canning 

Ames  Aeronautical  Laboratory,  NACA 

Aeroballistic  Range  Measurements  of  the  Performance 
and  Stability  of  a  Supersonic  Fighter  Aircraft 
(CONFIDENTIAL) 

H.  R.  Warren*,  ft...  J.  Templin'**  and  B.  Cheers 
Canadian  Armament  Research  and  Development 
Establishment 

Model  Launching  Techniques  and  Other  Items  Related 
to  Range  Firings 
J.  E.  Long 

Naval  Ordnance  Laboratory 

Closing  Remarks  3:20  -  3:30  pm 


Mr.  C.  L.  Poor 


Saturday,  26  January  1957 


Tour 


10: 00  -  1:00  pm 


Naval  Ordnance  Laboratory 
Aeroballistic  Facilities 


* 

DeHavilland  Aircraft  of  Canada  Limited. 
National  Aeronautical  Establishment. 
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UNCLASSIFIED  ABSTRACTS 


SIMULATION  OF  THE  ATMOSPHERIC  ENTRY  OF  BALLISTIC  MISSILES  II:  15 

(CONFIDENTIAL) 

S.  E.  Neice  and  J.  A.  Carson 
Ames  Aeronautical  Laboratory,  NACA 


Abstract  is  Classified 


SOME  PROBLEMS  ASSOCIATED  WITH  THE  DETERMINATION,  FROM  II:  37 

RANGE  FIRINGS,  OF  DYNAMIC  STABILITY  OF 
BALLISTIC  MISSILE  RE-ENTRY  SHAPES 

(CONFIDENTIAL) 

L.  C.  MacAllister 
Ballistic  Research  Laboratories 

In  the  past,  free  flight  ranges  have  been  useful  in  the 
determination  of  the  dynamic  stability  of  bodies  of  revolution 
and  of  symmetric  missiles.  Recently  a  considerable  amount  of 
work  has  been  devoted  to  firings  of  models  of  var  heads.  The 
conditions  under  which  the  models  are  fired  and  the  aerodynamic 
properties  of  the  shapes  make  the  determination  of  the  dynamic 
stability  of  the  models  quite  difficult.  The  problems  that 
arise  and  some  possible  solutions  ar,e  discussed. 


A  CORRELATION  OF  FREE-FLIGHT  TRANSITION  MEASUREMENTS  ON  I:  25 

VARIOUS  BLUNT  NOSE  SHAPES  BY  USE  OF  THE 
MOMENTUM-THICKNESS  REYNOLDS  NUMBER 

W.  R.  Witt,  Jr.  and  J.  Pers^ 

U.  S.  Naval  Ordnance  Laboratory 

A  systematic  series  of  blunt  nose  shapes  has  been  fired  in 
the  Pressurized  Ballistics  Range  for  boundary- layer  transition 
studies.  The  transition  of  the  boundary- layer  flow  from  laminar 
to  turbulent  is  determined  directly  from  the  shadowgraph  plates . 

The  nose  shapes  have  all  been  fired  near  a  Mach  number  of  5  and 
the  Reynolds  number  per  foot  has  been  varied  by  changing  the 
pressure  (density  In  the  firing  range). 


IT 


The  Reynolds  number  based  on  momentum  thickness,  Re^,  at  the 
observed  transition  location,  was  calculated  using  the  laminar 
boundary- layer  calculation  method  given  by  Cohen  and  Reshotko  in 
"The  Compressible  Laminar  Boundary  Layer  with  Heat  Transfer  and 
Arbitrary  Pressure  Gradient"  (NACA  TN  3326).  In  general,  the 
results  indicate  that  transition  occurs  at  values  of  ReQ  which 
are  of  the  same  order  of  magnitude  as  the  values  of  minimum 
critical  Reynolds  number  usually  associated  with  incompressible 
flow. 


ADVANCES  IN  THE  DYNAMIC  ANALYSIS  OP  RANGE  DATA  I:  45 

C.  H.  Murphy 

Ballistic  Research  Laboratories 

The  range  technique  has  been  usually  restricted  to  dynamic 
analysis  of  the  motion  of  thrustless  symmetric  missiles  acted  on  by 
linear  aerodynamic  forces.  In  recent  years  all  three  of  these 
restrictions  have  been  relaxed. 

Pirst,  a  gun-boosted  burning  rocket  program  fired  on  the 
Transonic  Range  is  described  and  the  relatively  minor  alterations 
to  the  data  analysis  procedure  are  indicated.  Next  the  more 
difficult  problem  of  a  finned  missile  with  bent  fins  and  spin 
rate  varying  through  resonance  is  discussed.  Finally,  the 
successful  treatment  of  cubic  nonlinearities  in  static  and  Magnus 
moments  and  their  associated  forces  is  described. 


SURVEY,  CALIBRATION,  AND- REDUCTION  TECHNIQUES  USED  I:  85 

AT  THE  THOMPSON  AEROBALLISTICS  LABORATORY 

W.  H.  Allan  and  E,  L.  Dunn 
U.  S.  Naval  Ordnance  Test  Station 

The  cameras  of  the  Thompson  Aeroballistics  Laboratory  are 
calibrated  once  a  year  by  photographing  reference  markers  in 
the  field  of  view  of  each  camera.  Three  wires,  anchored  to 
towers  at  each  end  of  the  range,  are  suspended  in  space  near 
the  range  line.  The  reference  markers  are  l/8-inch  plastic 
beads' Tocated  every  two  feet  along  the  wire.  This  paper  describes 
the  techniques  used  to  measure  the  coordinates  of  the  cameras 
and  the  plastic  beads;  the  mathematical  treatment  of  these 
coordinates  to  furnish  calibration  equations  for  each  camera; 
the  use  of  the  equations  In  solving  for  orientation  and  location 
of  the  missile  in  space;  the  results  obtained  through  this 
system;  and  a  discussion  of  the  limitations  and  advantages  of  this 
system. 
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WAKE  VISUALIZATION  STUDIES  IN  THE  AEROBALLISTIC S  RANGE 


1:127 


G.  V.  Bull  and  C.  B.  Jeffery 
Canadian  Armament  Research  and 
Development  Establishment 

By  the  interaction  of  hydrochloric  acid  and  ammonium  hydroxide 
vapours,  a  plane  sheet  of  laminar  smoke  filaments  can  be  produced  in 
the  range  along  the  flight  trajectory.  Models  developing  lift  due 
to  incidence  were  fired  through  these  planar  sheets.  Spark  and  fastax 
photography  was  used  to  record  the  development  of  the  wake  profiles  in 
the  plane  of  the  smoke.  Wake  distortions  and  vortex  formations  have 
been  studied  for  several  types  of  bodies;  for  a  cruciform  arrangement 
of  rectangular  panels  of  aspect  ratio  1.9  on  a  cylindrical  body,  the 
wake  distortion  as  determined  from  these  tests  have  been  compared 
with  computations  based  on  the  assumptions  of  linear  theory. 


THE  CONTROLLED-TEMPERATURE-PRESSURE  RANGE  I;  1^9 

F.  D.  Bennett 

Ballistic  Research  Laboratories 

A  survey  is  given  of  the  research  and  development  program  which 
has  culminated  in  operation  of  the  Controlled-Temperature-Pressure 
Range  (CTFR)  for  production  of  flows  up  to  Mach  11.  Methods  of 
control  of  temperature  and  pressure  in  the  45*  working  section  are 
described.  The  instrumentation  necessary  for  (l)  measurement  of 
projectile  drag  coefficient  and  (2)  measurement  of  density  through¬ 
out  the  field  of  flow  is  discussed  in  some  detail. .  The  10"  Mach- 
Zehnder  interferometer  is  a  special  feature.  Various  research 
problems  encountered  in  the  development  of  light  sources,  projectile 
launchers  and  data  handling  schemes  are  briefly  touched  upon. 

At  low  supersonic  Mach  numbers  a  problem  requiring  the  full 
field  of  the  10"  Mach-Zehnder  interferometer  has  recently  been 
completed.  Here  a  study  of  phenomena  in  the  distant  N-wave  flow 
about  a  small  sphere  has  led  to  a  new  and  simple  experimental 
criterion  for  N-wave  flow  and  to  information  about  convergence 
to  N-wave  flow  with  radial  distance  from  the  projectile. 
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TWO  AEROBALLISTIC  RANGE  TOPICS 


I:  185 


J.  D.  Nicolaides 
Bureau  of  Ordnance 

(l)  MASS  ASYMMETRY 
with 

J.  E.  Long,  U.  S.  Naval  Ordnance  Laboratory- 
Gene  Parrish,  Bureau  of  Ordnance 

A  simple  approximate  theory  for  the  free  flight  motion  of 
ballistic  missiles  having  mass  asymmetry  is  given  and  proofed  by 
experimental  firings  in  the  NOL  Pressurized  Aeroballistic  Range. 

(2)  DYNAMIC  STABILITY 

The  Epicyclic  Theory  for  the  flight  dynamics  of  ballistic 
missiles  has  yielded  various  "Dynamic  Stability  Criteria"  which 
are  often  used  to  evaluate  missile  performance.  Recent  mislead¬ 
ing  uses  of  the  theory  and  criteria  in  appraising  missile  perform¬ 
ance  require  a  simple  restatement  of  the  theory,  its  assumptions 
and  its  use. 

The  parameters  of  Nutation  Half-Life,  Precession  Half Life 
and  Total  Motion  Half-Life  are  suggested  as  better  criteria  for 
missile  dynamics  than  those  classically  based  on  the  Linear  Theory. 

Also  a  summary  of  important  Nonlinear  Cases  of  ballistic 
missile  flight  performance  is  given. 


DESIGN  AND  INITIAL  TESTS  OF  TEE  NOL  SHOCK  GUN  II:  6l 

(CONFIDENTIAL) 

V.  C.  D.  Dawson 

U.  S.  Naval  Ordnance  Laboratory 

The  design  and  initial  tests  of  the  NOL  Shock  Gun  are  described. 
The  operation  of  this  gun,  which  was  conceived  by  Dr.  A.  E.  Seigel 
and  Dr.  Z.  I.  Slawsky.  Is  based  upon  a  new  principle  and  missiles 
weighing  two  grams  have  been  launched  from  a  0. 50-caliber  smooth¬ 
bore  gun  at  velocities  in  excess  of  15,000  feet  per  second. 
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M  APPLICATION  OP  AEROBALLISTICS  RANGE  TECHNIQUES 
TO  MODELS  WITH  DIGONAL  ROTATIONAL  SYMMETRY 


1:219 


G.  H.  Tidy  and  M.  E.  Thomas 
Canadian  Armament  Research  and 
Development  Establishment 

A  series  of  flat  plate  wings  of  triangular  plamforms  has 
been  fired  at  Mach  numbers  1.5  and  2  and  their  trajectories 
have  been  measured. 

Preliminary  manual  reduction  of  the  data  is  presented  and 
the  derived  values  of  some  aerodynamic  coefficients  are  compared 
with  NACA  wind  tunnel  measurements.  The  possibility  of  more 
complete  analysis  and  of  application  of  the  range  technique  to 
airplane  configurations  are  considered. 


SABOTS  USED  AT  THE  THOMPSON  AEROBALLISTICS  LABORATORY  I:  24l 

W.  H.  Allan 

U.  S.  Naval  Ordnance  Test  Station 

A  review  of  sabots  used  at  the  Thompson  Aeroballistics 
Laboratory  from  the  beginning  of  operations  to  the  present. 

The  discussion  will  cover  spinner  and  finner  sabots  used  in 
guns  ranging  from  LO-mm  to  8- inch  bore  diameter.  Follow- thru, 
breakapart,  slug  styrofoam,  and  slow-spin  smooth  bore  sabots 
will  be  discussed  along  with  the  use  of  the  sabot  retarder. 


RESEARCH  INVESTIGATIONS  IN  THE  AMES  II:  8l 

SUPERSONIC  FREE-FLIGHT  FACILITIES 

T.  N.  Canning 

Ames  Aeronautical  Laboratory,  NACA 


Abstract  is  classified. 
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AEROBALLISTIC  RANGE  MEASUREMENTS  OF  THE  II:  101 

PERFORMANCE  AND  STABILITY  OF  A  SUPERSONIC 
FIGHTER  AIRCRAFT 

(CONFIDENTIAL) 

H.  R.  Warr6rv*,  R.  J.  Templin**,  and  B.  Cheers 
Canadian  Armament  Research  and 
Development  Establishment 

This  paper  describes  a  method  being  developed  for  measuring 
aircraft  performance  and  stability  characteristics  in  free  flight. 

T^sts  have  been  made  firing  Into  the  Aeroballistlcs  Ranges  smalJ 
scale  models  of  a  current  delta  wing  fighter  at  a  supersonic  Mach 
number  and  approximately  1/10  its  combat  Reynolds  number.  Velocity 
screens,  schlieren  and  yaw  card  measurements  are  used  to  obtain 
histories  of  the  models  speed,  altitude  and  later  motion  during 
flight.  From  the  analysis  of  these  records  Information  la  Obtained 
about  the  aircraft  drag,  lift,  lateral  and  longitudinal  aerodynamic 
derivatives . 


* 

DeHavilland  of  Canada,  Limited 

-**■ 

National  Aeronautical  Establishment 


MODEL  LAUNCHING  TECHNIQUES  AND  OTHER  ITEMS  I:  275 

RELATED  TO  RANGE  FIRINGS 

J.  E.  Long 

U.  S.  Naval  Ordnance  Laboratory 

A  discussion  of  the  variety  of  methods  used  to  launch  scaled 
models  in  the  free-flight  precision  ranges  at  NOL  is  given.  These 
methods  include  such  items  as:  (l)  launching  finned  missiles  from 
rifled  guns;  (2)  launching  subcaliber  spinning  models  from  over¬ 
sized  sabots;  (3)  launching  spheres  as  small  as  l/32  inch  in  diameter 
for  drag;  and  (4)  launching  model  aircraft. 

Under  related  items  the  discussion  will  be  centered  about  the 
following  techniques:  (l)  firing  models  with  a  jet  exhausting  from 
the  model  base;  (2)  investigating  the  arming  of  fuzes  by  X-raying 
the  recovered  round;  (3)  firing  spinning  models  with  hot  and  cold 
plastics  rotating  hands;  (4)  development  of  the  spin  sonde;  and 
(5)  firing  models  from  powder  guns  at  10,000  ft /sec. 
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A  CORRELATION  OF  FREE-FLIGHT  TRANSITION 
MEASUREMENTS  ON  VARIOUS  BLUNT  NOSE  SHAPES 
BY  USE  OF  THE  MOMENTUM-THICKNESS  REYNOLDS  NUMBER 


W.  R.  Witt,  Jr. 

J.  Persh 

U.  S.  Naval  Ordnance  Laboratory 
■White  Oak,  Silver  Spring,  Maryland 


INTENTIONALLY  LEFT  BLANK, 


A  CQRREIATXON  <F  FREE -FLIGHT  TRANSITION 

measurements  on  various  blunt  nose  shapes 

BY  USE  OF  THE  MOMENTUM- THICKNESS  REYNOLDS  NUMBER* 

INTRODUCTION 

One  of  the  greatest  uncertainties  in  the  prediction  of  the  heat 
transfer  to  arbitrary  body  shapes  is  the  location  of  the  transition 
"point."  While  laminar  boundary- layer  theory  yields  indications  of 
the  probable  effects  of  several  variables  on  boundary- layer  transition, 
it  does  not  permit  prediction  of  the  transition  "point." 

For  some  time  ve  have  been  investigating  boundary- layer  transition 
in  the  Pressurized  Ballistics  Range  at  the  Naval  Ordnance  Laboratory. 
The  transition  "point"  is  determined  directly  from  the  shadowgraph 
plates  and  is  taken  as  the  position  at  which  small  eddies  first  appear 
in  the  boundary  layer. 

Our  earlier  investigations,  sponsored  by  the  U.  S.  Air  Force,  have 
previously  been  reported  (reference  a) .  These  firings  were  made  with 
10-degree  included  angle  cones  that  had  polished  or  ground-finished 


*  Published  in  report  form  at  NOL  as  NavOrd  4400.  Supported  by 
the  U.  S.  Air  Force 
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surfaces,  and  sharp  or  slightly  blunted  tips.  Results  of  these 
earlier  Investigations  were  as  foil  oust 

a.  For  a  highly-polished  (l  to  2  micro  inch)  sharp-pointed 
cone,  transition  uas  found  to  occur  at  free-Btream  Reynolds  numbers 
varying  from  1  million  to  21  million  at  small  angles  of  attack  under 
similar  firing  conditions 

b.  On  a  rougher  cone  (8  to  10  microinches)  vith  a  slightly 
blunted  tip  (0.0625  Dia.)  the  same  maximum  Reynolds  numbers  of 
transition  vere  obtained  as  on  the  polished  cone.  The  low  values 
of  Reynolds  number  of  transition,  however,  observed  on  the  smooth 
cone  were  not  obtained  with  the  slightly  blunted,  rougher  cone. 

The  recent  tests  which,  are  to  be  described  here  were  also 
fired  in  the  Pressurized  Ballistics  Range.  This  range  is  well-suited 
to  the  investigation  of  boundary- layer  behavior  for  the  following 
reasons: 

a.  The  turbulence  level  is  negligible 

b.  The  air  pressure  (density)  can  be  varied,  permitting  a 
wide  range  of  Reynolds  number  per  foot  (much  larger  than  can  usually 
be  obtained  in  wind  tunnels) 

c.  Numerous  spark  shadowgraph  stations  are  available  permitting 
both  horizontal  and  vertical  observations. 

The  quality  of  the  results  obtained  from  a  boundary- layer  study 
in  free  flight  will  depend  on  how  veil  the  following  parameters  can 
be  controlled  or  known  accurately: 

a.  Mach  number 

b.  Body  geometry  or  pressure  distribution 
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c.  Range  pressure  (Re/ft) 

d .  Surface  finish 

e.  Angle  of  attack 

f.  Heat  transfer  (Ty/T^) 

The  present  study  was  made  with  the  models  shown  on  Figure  1 
with  most  of  these  six  parameters  held  constant.  These  nose  shapes 
include  20-degree  and  UO-degree  cones  with  nose-to-base  radii  of 
0.75^  0.500,  0.250,  and  0.125  and  cone-like  noses  based  on  1/2  and 
l/3  power  curves.  For  these  firings  the  Mach  number  has  been  in  the 
neighborhood  of  3  and  the  range  pressure  was  varied  from  l/2  to  2-1/4 
atmospheres,  with  one  shot  at  5  atmospheres. 

For  the  firings  discussed  herein,  the  wall  temperature  of  the 
models  was  essentially  room  temperature.  However,  equipment  haB 
been  developed  which  will  permit  the  temperature  of  10-degree  cones 
to  be  raised  as  high  as  1000  degrees  Fahrenheit  immediately  before 
firing. 

Complete  control  of  the  angle  of  attack  is  impossible  in  range 
firings.  The  models  for  transition  studies  are  always  designed  so 
as  to  keep  these  angles  small  (0<tcn:4o).  Since  the  angle  of  attack 
cannot  be  controlled  accurately,  it  can  be  treated  as  a  correlation 
parameter  (see  reference  a). 

Surface  roughness  cannot  be  eliminated  but  polishing  techniques 
can  reproduce  a  smooth  surface  rather  well.  A  number  of  the  l/2  power 
and  l/3  power  noses  were  polished  to  a  1  to  2  microinch  finish.  The 
cones  were  unpolished  (8  -  10  microinches).  A  detailed  account  of  the 
polishing  technique  and  of  the  sui*face-measuring  method  can  be  found 
in  reference  (a). 
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Pressurized  Ballistics  Range 

The  Pressurized  Ballistics  Range  (Figure  2)  is  an  enclosed 
tube  three  feet  in  diameter  and  300  feet  long  in  which  the  pressure 
can  be  varied  from  0.01  to  6  atmospheres.  There  are  twenty- five 
stations  at  which  spark  shadowgraphs  can  be  taken  of  the  model  in 
flight.  As  the  model  approaches  each  station  (Figure  3)  a  spark 
is  initiated  by  the  pulse  from  a  photoelectric  pickup  unit.  At 
each  station  the  light  from  the  spark  fa]_ls  directly  on  a  vertical 
photographic  plate  and  is  reflected  by  a  mirror  onto  a  horizontal 
plate.  The  times  at  which  the  sparks  flash  at  any  13  of  the  25 
stations  are  recorded  on  electronic  counters. 
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Symbols  and  Definitions 


Mach  number 

free -stream  stagnation  pressure 
Reynolds  number  based  on  ramentum  thickness  (©) 
local  radius  of  model 
radius  at  base  of  model 
distance  measured  along  surface  of  model 
local  static  temperature 
longitudinal  velocity  component 
distance  perpendicular  to  surface 
boundary- layer  momentum  thickness 


<*y 


half -angle  of  cones 

viscosity  coefficient 
kinematic  viscosity  coefficient 


density 
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Subscripts : 

e  =  local  conditions  outside  boundary  layer 

v  =  vail  or  surface  value 

e>o  =  initial  or  free -stream  value 

Data  Recorded 

The  two  spark  shadowgraphs  obtained  from  each  station  give 
horizontal  and  vertical  projections  of  the  model  and  of  the  flow 
pattern.  The  transition  is  veil  defined  except  in  the  optically 
distorted  region  near  the  front  of  blunt-nosed  bodies.  As  mentioned 
before,  transition  is  detected  by  locating  the  formation  of  small 
eddies  or  vortices  in  the  boundary  layer.  From  each  plate  the 
position  of  transition  is  measured  both  on  the  windward  and  sheltered 
sides  of  the  cone.  Thus  there  are  four  transition  readings  at  each 
station.  Three  typical  spark  shadowgraph  prints  are  shown  as 
Figures  4,  5,  and  6.  The  true  angle  of  attack  is  calculated  from 
the  horizontal  and  vertical  components  of  the  angle  of  attack  that 
are  read  from  the  plates.  Since  neither  of  the  shadowgraph  plates 
is,  in  general,  parallel  to  the  plane  of  the  angle  of  attack,  the 
fully  windward  and  fully  sheltered  sides  are  not  seen.  The  fully 
windward  side  can  be  determined  and  the  four  transition  readings  at 
each  station  can  be  related  (in  azimuth)  to  the  true  wind  vector. 

Since  from  one  shot  sixty  transition  measurements  are  normally 

obtained,  one  representative  transition  length  had  to  be  selected 
for  each  shot.  This  was  done  by  taking  an  average  of  the  readings 

from  the  vertical  plates.  The  vertical  plate  readings  were  selected 
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since  these  readings  are  more  accurate  than  the  readings  from  the 
horizontal  plates  because  the  reflected  light  from  the  mirror 
(Figure  3)  does  not  give  as  good  a  shadowgraph. 

Results 

It  appears  reasonable  that  any  correlation  of  boundary- layer 
transition  on  body  shapes  vith  pressure  gradients  should  involve 
Reynolds  numbers  based  on  flow  properties  outside  the  boundary  layer 
and  some  characteristic  thickness  associated  with  the  boundary  layer. 

Of  the  boundary- layer  parameters  considered  (total  thickness, 
displacement  thickness,  and  momentum  thickness),  the  boundary-layer 
momentum  thickness  appears  most  reasonable.  The  momentum  thickness 
may  be  defined  as  the  thickness  of  a  layer  in  the  free-stream  flow 
having  a  momentum  equal  to  the  loss  of  momentum  in  the  boundary  layer 
due  to  viscous  effects.  Nov  for  a  look  at  the  results  and  the  success 
of  the  correlation. 

The  transition  data  discussed  here  were  obtained  from  thirty-one 
shots  including  sixteen  firings  with  the  l/2  and  l/3  power  noses. 

The  data  from  these  firings  are  listed  on  Table  I.  The  variation  of 
ReQ  with  distance  along  the  surface  was  calculated  tip  to  the  observed 
transition  point  using  the  method  given  by  Cohen  and  Reshotko  in 
reference  (b ) .  The  use  of  this  method  requires  a  knowledge  of  the 
surface  pressure  and  temperature  distributions  at  the  edge  of  the 
boundary  layer  in  addition  to  the  vail  temperature.  The  pressure 
distribution  was  determined  by  use  of  the  modified  Newtonian  flow 
concept  (reference  c),  and  the  temperature  distribution  was  calculated 
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using  the  perfect  gas  isentropic  flow  relations.  While  it  is 
realized  that  the  Nevtonian  flow  calculation  procedure  is  inadequate 
for  the  slender  cone  models,  it  was  felt  that  for  consistency  it 
.should  be  used  for  all  of  the  configurations  investigated.  Hie 
model  wall  temperature  was  assumed  to  be  room  temperature  and 
unchanged  during  the  short  time  of  flight  (duration  of  less  than  0.1 
seconds).  This  is  a  reasonable  assumption  since  the  calculation 
scheme  utilized  is  not  particularly  sensitive  to  small  changes  in 
wall  temperature.  It  is  important  to  mention  that  the  transition 
"point"  Reg  values  obtained  are  presumably  associated  with  the  end 
of  transition  even  though  they  are  based  on  the  average  of  numerous 
readings  per  shot.  Justification  for  the  statement  that  the  first 
appearance  of  turbulence  on  the  photographs  indicates  the  end  of 
transition  was  reported  by  Grunevald  (reference  d).  In  Grunevald ’ s 
wind-tunnel  investigation  it  was  shown  by  comparisons  between  surface 
temperature  data  and  schlieren  observations  that  the  end  of  transition 
is  observed.  Effectively,  if  this  is  the  actual  case,  this  means 
that  the  transition  region,  which  may  be  of  appreciable  length,  is 
neglected  and  it  is  assumed  that  laminar  flow  exists  in  this  region. 
Since  the  length  of  the  transition  region  depends  on  a  combination  of 
several  parameters,  it  is  not  possible  to  assign  a  specific  adjustment 
factor  to  the  data  which  are  presented.  However,  Mr.  Persh  of  NOL  has 
a  HavOrd  report  ready  for  publication  (reference  g)  in  which  he 
describes  a  technique  for  calculating  the  boundary- layer  development  in 
the  transition  region. 
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Figure  7  shows  how  the  values  of  Re^^,  determined  as  previously 

described,  vary  with  the  incompressible  Pohlhausen  (reference  f) 

9  due 

pressure  gradient  parameter,  Re^  ^  -j^r.  This  pressure  gradient 

parameter  was  chosen  as  the  abscissa  because  of  its  familiarity  and 
apparent  suitability  as  a  correlation  parameter  for  compressible  flows. 
Also  included  in  this  figure  are  the  subsonic  data  reported  by  Persh 
(reference  e).  It  is  important  to  point  out  that  the  subsonic  data 
are  for  the  start  of  transition  and,  therefore,  represent  an  earlier 
part  of  the  transition  regime  than  do  the  range  data.  It  is  seen  that 
all  of  the  experimental  values  of  Re^p  lie  between  300  and  800,  and 
are  reasonably  grouped  about  the  theoretical  incompressible  curve  of 
Tetervin  and  Levine  (reference  f).  The  favorable  influence  of  strong 
heat  transfer  and  strong  pressure  gradient  on  laminar  stability  is  not 
evident  from  these  data.  In  fact,  the  data  indicate  that  neither  a 
cool  surface  nor  favorable  pressure  gradient  appreciably  affect  the 
values  of  Reg  at  the  end  of  transition.  In  general,  it  appears  that 
transition  on  these  body  shapes  occurs  at  Reynolds  numbers  which  are 
of  the  same  order  of  magnitude  as  the  values  of  minimum  critical 
Reynolds  number  usually  associated  with  incompressible  flow. 

Conclusions 

From  a  systematic  investigation  of  boundary- layer  transition  on 
blunt-nose  body  shapes  in  the  NOL  Pressurized  Ballistics  Range,  the 
following  conclusions  may  be  drawn: 

a.  For  all  the  data  analyzed,  the  calculated  values  of  the 
momentum  thickness  Reynolds  number,  ReQ,  at  the  observed  transition 
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point  were  never  larger  than  800  nor  less  than  300.  It  appears 
therefore  that  transition  occurs  at  values  of  ReQ,  which  are  of 
the  same  order  of  magnitude  as  the  values  of  minimum  critical  Reynolds 
number  usually  associated  with  incompressible  flow 

t.  Although  the  data  have  a  "trend”  with  incompressible  pressure 

6  due 

gradient  parameter  Re^  Uq  ds  suggestive  of  that  described  by 
incompressible  stability  theory,  it  does  not  appear  that  a  favorable 
pressure  gradient  appreciably  affects  the  values  of  ReQ  at  which 
transition  occurs 
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Fig.  3  Typical  Spark  Station,  Pressurized  Ballistics  Range 
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FIG. 7  VARIATION  OF  CALCULATED  VALUES  OF  Reg 
WITH  THE  INCOMPRESSIBLE  PRESSURE  GRADIENT 
PARAMETER,  Refl-S-  ^ 
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ADVANCES  IN  THE  DYNAMIC  ANALYSIS  OP  RANGE  DATA 


Charles  H .  Murphy 
Ballistic  Research  Laboratories 

INTRODUCTION 

The  first  use  of  aerodynamic  ranges  was  the  accurate  measurement 
of  the  drag  coefficients  of  bullets  and  shells .  As  the  instrumentation 
improved  in  both  quantity  as  well  as  quality,  good  measurements  of  the 
pitching  and  yawing  motion  of  projectiles  became  available.  The  trans¬ 
verse  motion  of  their  centers  of  mass  were  also  measured  with  good 
accuracy.  From  these  data  values  of  various  linear  aerodynamic  coeffi¬ 
cients  could  be  computed.  Ranges  were  then  used  to  test  symmetric  finaed 
missiles.  Since  finned  missiles  with  angles  of  rotational  symmetry 
equal  to  or  less  than  120°  have  the  same  linear  aerodynamic  force  and 
moment  expansion  as  bodies  of  revolution,  this  introduced  no  new 
difficulties. 

Thus,  it  can  be  said  that  the  "classical"  use  of  aerodynamic 
ranges  has  been  for  the  dynamic  analysis  of  thrustless  symmetric  missiles 
acted  on  by  linear  aerodynamic  forces  and  moments.  All  three  of  these 
restrictions  arise  from  the  fact  that  the  range  method  infers  the 
aerodynamic  force  and  moment  from  the  observed  motion  and,  therefore, 
is  usually  restricted  to  problems  for  which  the  differential  equations 
have  solutions  in  closed  forms.  As  we  shall  see  these  restrictions 
can  be  and  have  been  relaxed.  Although  the  swerving  motion  of  the  center 
of  mass  can  be  considered,  we  will  limit’  ourselves  in  this  paper  to  the 
pitching  and  yawing  motion  and  its  associated  aerodynamic  coefficients. 

This  paper  will  be  divided  into  four  parts.  First  we  will  discuss 
the  "classical"  case.  Then  a  gun-boosted  burning  rocket  program  fired 
on  the  Transonic  Range  will  be  described  and  the  minor  alterations  to 


the  data  analysis  indicated.  Next  the  more  difficult  problem  of  an 
asymmetrical  finned  missile  with  roll  rate  varying  through  resonance 
is  considered.  Finally,  the  successful  treatment  of  cubic  non-linear¬ 
ities  in  static  and  Magnus  moments  and  their  associated  forces  is 
described. 
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THE  CLASSICAL  THEORY 


Although  the  linear  expansion  of  the  transverse  aerodynamic  force 

and  moment  in  terms  of  the  angles  a,  p,  the  angular  velocities  q,  r, 

and  their  time  derivatives  require  the  use  of  thirty- two  derivatives, 

the  assumed  rotational  symmetry  reduces  this  number  to  sixteen.  The 

2 

introduction  of  a  plane  of  mirror  symmetry  reveals  that  half  of  these 

derivatives  are  odd  functions  of  the  roll  rate  p  and,  therefore,  vanish 

for  zero  roll  rate.  These  will  be  identified  by  products  of  second 

order  derivatives  and  the  roll  rate  and  are  called  Magnus  coefficients. 

The  remaining  derivatives  can  be  independent  of  roll  rate.  This  total 

of  eight  non-Magnus  and  eight  Magnus  derivatives  can  be  reduced  to  five 

if  only  those  quantities  which  usually  have  a  measurable  effect  on 

* 

the  motion  are  retained  .  (Due  to  the  rotational  symmetry  the  use  of 
complex  variables  is  a  convenient  shorthand.) 


where 


.\Y  +  iZ  =  -  i  p  CN  (p  +  ia)  (1) 


1 

2  p 

j  r  (M) 

°M  '  1CM  ]  [P  +  10 

l L 

pa  a  J  u 

r 

“ 

n 

r  (q  +  ±r)l  1 

(p  +  i &)i 

ipi  (p  +  ia) 

M 

L  V  J. 

"  1  M.  L  v  + 

v  J 

q 

a 

(2) 

- 

p  is  air  density 

V  is  total  velocity 

S  is  reference  area 

i  is  reference  length 

p,q,r  are  components  of  the  angular  velocity  vector 
a  is  angle  of  attack 

p  is  angle  of  sideslip 
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1  „ 
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-  -  cz 

(Normal  force  derivative) 

a 

0 
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M  ^ 
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m 

=  -  c 

(Static  moment  derivative) 
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m  „ 
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4. 

i  _ 

C  =  -  c  1 

*  (Damping  moment  derivatives) 

M&  " 

md 

% 

(The  expression  following 
rolling  coordinate  system. 


Cj^  is  the  derivative  of  (3  +  ia  in  a  non- 
As  has  been  pointed  out  by  Sacks^  and  later 


by  Davis 
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this  choice  has  certain  advantages.) 


In  Eqs.  (l-2)  both  the  angles  and  the  angular  velocities  are  measured 
in  missile -fixed  coordinates.  Since  non-rolling  coordinates  give  rise  to 
slightly  less  complex  differential  equation,  we  will  make  use  of  a  non¬ 
rolling  coordinate  system.  Quantities  in  this  coordinate  system  will  be 
distinguished  by  tilde  superscripts.  (Fig.  l) 


If  the  effect  of  gravity  is  neglected,  the  aerodynamic  force  assumed 
to  be  linear,  and  the  mass  distribution  assumed  to  be  symmetric  (ly  =  Iz), 
the  following  differential  equation  for  the  complex  angle  of  attack,  Tj 
can  be  derived: 

T"  +  (H  -  iP)'|*  -  (M  +  iFT)  T  =  0  (5) 

where  T*  =  fT  +  ia  =  (p  +  ia)  P  ^ 
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m  =  mass 

CD  drag  coefficient 


transverse  radius  of  gyration 


axial  radius  of  gyration 


I  axial  moment  of  inertia 
x 


I  -I  transverse  moment  of  inertia  and  primed  denote  differentia- 
y  z  r  y 

tion  with  respect  to  dimensionless  arc  length  s  ~  For  constant 

spin  to  velocity  ratio  and  constant  moment  coefficients,  Eq.  (3)  can  "be 

easily  solved: 


£  =  I^e  *  rl  +  e1  02 


1  01  +  K2  eJ 


■(h) 


where  0^  =  0^o  +  0^  3 

KJ  =  KJo  -  * 


and 


Xj  +  i0j  =  i  [_■  H  +  iP  +  H2  -  P2  +  2iP  (2T  -  "h)J 


The  pitching  and  yawing  motion  is  described  in  a  complex  plane  by 
the  Siam  of  two  two-dimenBional  vectors  which  are  exponentially  damped 
in  magnitude  and  rotating  with  angular  velocities  0*.  (Fig.  2)  For  a 
statically  stable  missile  with  zero  spin  (M<0,  P  =  0)  the  motion  is 
described  by  a  damped  ellipse.  (Fig.  3)  If  the  missile  is  allowed  to 
spin,  the  ellipse  will  precessand  become  an  epicycle.  (Fig.  4)  A 

2 

statically  unstable  missile  which  is  spin-stabilized,  (M>0,  4m  -  P  <0), 
has  an  epicyclic  motion  with  the  nodes  on  the  inside.  (Fig.  5) 
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The  basic  procedure  in  range  tests  is  to  fit  observed  pitching 

and  yawing  notion  by  Eq.  (t)  and  to  infer  from  the  measured  damping 

factors  and  frequencies  the  actual  aerodynamic  moments  which  act  on 

the  missile*  In  a  similar  fashion  the  center  of  mass  motion  can  be 

used  for  measurements  of  the  normal  force  derivative  C.-  . 

N 
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THRUST 


The  consideration  of  a  burning  rocket  introduces  three  modifications 
of  Eq.  (5): 

(1)  The  drag  equation,  C^,  is  no  longer  valid. 

(2)  The  mass  and  radii  of  gyration  are  not  constants. 


(3)  A  jet  damping  term*,  — -  1  jj 


where  r  is  distance 
e 


from  the  c.m.  to  the  nozzle  exit,  must  be  added  to  the  right 
side  of  Eq.  (3)* 

A  further  practical  handicap  to  aerodynamic  range  tests  of  burning 
rockets  is  the  large  dispersion  inherent  is  conventional  rockets. 


* 

Live  gun-boosted  rockets  have  tolerable  dispersions  and  have  been 
launched  in  the  Transonic  Range^ .  Figure  6  is  a  direct  microflash  of  a 
burning  2.75  in.  rocket  while  Figure  7  is  a  shadowgraph  of  the  same 
rocket.  Although  the  rocket  motor  does  expose  part  of  the  shadowgraph, 
the  shadow  itself  is  of  good  quality. 


For  this  firing  the  modifications  listed  above  can  be  easily 
handled: 

(1)  -y-  was  approximately  constant  and  positive.  This  meawt  that 
a  negative  pseudo  drag  coefficient  may  be  used  in  Eq.  (3) 

(2)  The  variations  in  mass  and  radii  of  gyration  were  relatively 
mild  and  averages  can  be  estimated  to  within  10$. 

(3)  The  {jet  damping  term  for  the  high  initial  velocities  was  small 
and  so  could  be  neglected. 


The  conventional  data  analysis  was  then  employed  to  yield  values  of  C„ 

% 
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This  term  arises  from  the  fact  that  the  gases  in  the  jet  carry  angular 
momentum  from  the  rocket. 
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CjW  \ +  \  ,  and  CK  .  Figure  8  gives  a  comparison  of  for 

■burning  and  non-turning  rockets.  The  .Jet  appears  to  have  little 
effect  on  this  aerodynamic  coefficient. 


AERODYNAMIC  ASYMMETRY 


In  1952,  NicolaideJ  treated  the  problem  of  missiles  with  slight 
configurational  asymmetry  and  constant  roll  rates .  This  asymmetry 
can  be  caused  by  misalligned  fin  surfaces  or  deflected  control  surfaces 
and  can  be  described  by  a  constant  amplitude  aerodynamic  force  and 
moment  whose  direction's  are  constant,  in  missile  -  fixed  coordinates. 

For  this  case  our  basic  equation  (Eq.  (3))  is  modified  by  the  intro¬ 
duction  of  an  inhomogeneous  term. 

T"  +  (H  -  iP)  7*  -  (M  +  iPT)  7  =  e^9  +  9e^  (5) 

where  M  =  -  k~2  CM  e 

e  2m  t  M£ 

e  is  magnitude  of  asymmetry  angle 
0£  is  its  initial  orientation  and 


t  s 


For  constant  ^  Nicolaides  showed  that  the  solution  of  the  pitching  and 
yawing  motion  assumed  the  form: 


7=1^  e1^  +  K2  e  ^2  +  K  e  +  Ho* 


(6) 


where 


K,  e^50  = 

3 


-q  Q 

9T  =  and  the  other  symbols  are  defined  after  Eq.  (4). 


-Me  e 
e 


0fe2  (1-  Y~)  +  M  -  10* (H-  T) 
7  7 


The  third  term  in  Eq.  (6)  possesses  a  rather  important  resonance  property 

for  statically  stable  missiles.  This  occurs  when  the  roll  rate  is  near  the 

1 2  j 

pitching  frequency.  (More  precisely  it  is  when  0  (  1  -  Tx)  +  M  =  0.) 

Ty 
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For  this  roll  rate  can  become  ^_quite  large.  Its  size  is  determined 
by  the  aerodynamic  ^.sj^ing,  H  -  ^  T.  In  Figure  9,  the  amplification  of 
equilibrium  yaw,  ,  is  plotted  versus  the  ratio  of  spin  rate  to 

yawing  rate  for  different  aerodynamic  damping.  The  resonance  peaks  can 
be  quite  large  for  reasonable  aerodynamic  damping. 


Naturally,  the  curves  of  Figure  9  are  of  no  quantitative  value 
for  the  variable  roll  rate  problem.  For  this  more  difficult  data  analysis 
problem  of  free  flight  missiles  whose  roll  rates  vary  through  resonance 

Q 

considerable  use  was  made  of  an  analog  computer  .  (Figure  10) 

For  this  method  Eq.  (5)  together  with  the  measured  roll  history 
was  placed  on  the  computer,  and  a  series  of  solution  curves  were  drawn. 
These  were  compared  with  actual  measurements  and  the  parameters*  H, 

M,  Mp  and  ©  together  with  initial  conditions  i  ,  were  then  varied  hy 
trial  and  error  in  a  curve  fitting  process.  For  an  experienced  operator 
a  good  fit  could  he  obtained  in  two  to  four  hours .  Figure  11  shows  such 
a  fit  for  a  missile  whose  roll  rate  grows  to  1.25  times  its  resonant 
value . 


This  approach,  however,  suffers  from  a  crucial  dependence  on  an 
experienced  operator.  If  it  is  possible,  it  is  desirable  to  have  the 
machine  make  the  necessary  decisions.  With  this  in  mind  it  was  decided 
to  try  to  make  use  of  one  of  BEL's  large  digital  machines.  To  do  this 
it  was  necessary  to  he  able  to  state  the  solution  of  Eq.  (5)  for  non¬ 
constant  ~. 

At  this  point  it  is  important  to  note  that  although  the  particular 
solution  is  quite  sensitive  to  roll  rates  near  resonance,  the  complement¬ 
ary  solution  is  relatively  insensitive.  This  means  that  we  can  either 
neglect  on  the  left  side  of  Eq.(5)  or  replace  it  by  an  average  value. 

The  problem  is  now  to  find  the  solution  of  an  inhomogeneous  linear  equation 


* 

T  was  neglected  in  this  work. 
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with  constant  coefficients.  Since  we  know  that  the  complementary 
solution  is  given  by  Eq.  (4),  the  general  solution  can  be  found 
by  the  method  of  variation  of  parameters. 


According  to  this  method  the 

are  assumed  to  be  functions  of  s  . 

i  0 

the  quantities  A ^  e  *jo 


constants  K.  ,  in  Eq.  (h) 
jo’  pJo  ^  x  ' 

For  convenience  we  will  introduce 


where 


e  (X1  +  1  *  s  +  A2  e  (X2  +  1  3 

is  the  particular  solution  due  to  asymmetry  and 
are  complex  functions  of  a. 


Differentiating  we  see  that 


=  (\+  10[)  Alfi(Xl  +  3  +  (X2+  10*  )  A2e 

+  A*  e  ^Xl+  i0{^  “  +  A2  e  ^X2  +  ^2^  S 
If  the  last  two  terms  are  set  equal  to  zero, 

A>  a  <  V  *  +  A-  e  (X2  +  ^  8  .  0 

and 


(X2  +  i0»)  h 
(8) 


(9) 


r* 

V 


(X1  ♦  10')  Al  e(Xl  +  8 


(x2+ 


i  fh) 


A2e(X2+ 


i0p 

(10) 


fi 


Differentiating  again  we 

1; » +  is*!)2  a, 

+  (X1  +  10’) 


obtain  the  following: 

e(Xl+  i0|)  s  +  +  ±0*)2  ^  e(x2+ 

e(X1+  i01)  s+  ^  +  e(X2+ 


10^)  a 

102 )  s 
(11) 
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If  Eqs.  (7,  10,  ll)  are  substituted  in  Eq.  (5)  and  it  is  recalled 
that  terms  without  come  from  the  solution  of  the  homogeneous  equation 
and  must  cancel,  the  following  equation  can  be  obtained: 

^  +  10* )  A*  e  (X1  +  %  +  (X2+  10*)  A *  e^2+  *$2)  R=  M^e1^0  +  V 

(12) 

Substituting  Eq.  (9)  in  Eq.  (12), 

M  e‘(Xl+  ^\)S  +  i^9  +  6e^ 


*4 


Xx  -  x2  +  i(0*_  -  02) 

M  e-(x2  +  i0J>+  i(e  +  ee) 

€ _ 

X2  -  X1  +  i(02  -  0J.) 


(13) 


(14) 


If  Eqs.  (13-14)  are  integrated  and  the  integrals  substituted  in  Eq.  (7), 
the  following  general  solution  can  be  obtained: 


where 


*S  -  e1|Z)l  +  K2  e1^  +  XA 

M^ejj  [  e(V  i0i)(s-r)_e(X2+  i^)(s-r)]  ^(r)^ 


(15) 


(X1  -  X2)  +  1(0*  -  0») 


33q.  (15)  can  now  be  fitted  by  the  usual  least  squares  differential 
corrections  method.  This  process  has  been  coded  for  the  QRDVAC  by 

Q 

P.  Yonda  . 

The  rounds  which  had  been  fitted  by  the  analog  computer  in  Ref.  8 
have  now  been  fitted  by  the  OKDYAC  and  the  results  are  given  in  Table  1. 
The  standard  error  a  of  the  OKDYAC  fit  is  about  60$  of  the  standard  error 
of  the  analog  computer  fit.  Since  the  latter  work  was  a  trial  and  error 
process  and  the  OKDYAC  method  did  not  depend  on  the  skill  of  the  machine 
operator,  this  superiority  of  the  digital  method  was  not  surprising.  As 
can  be  seen  from  the  other  entries  in  the  table,  the  analog  computer 
operator  was,  however,  very  capable. 
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TABLE  I 


Rd. 

No. 

^ofl/^Reson 

/  10 

im 

M£ 

Ref.  8 

ORDVAC 

3099 

.72 

2.0 

1.3 

3.4 

-0.1 

-4.8 

3088 

•91 

1.5 

1.0 

1.4 

-6.5 

-6.4 

3100 

1.00 

1.8 

1.1 

6.1 

-6.5 

+30.3 

3098 

1.02 

2.0 

1.1 

2.9 

-7-1 

-3.1 

3105 

1.08 

1.9 

1.3 

2.6 

-1.9 

+8.5 

3103 

1.09 

1.3 

1.3 

3.9 

+9.8 

+33.3 

3090 

1.10 

1.8 

1.5 

l.l 

-13.3 

+14.7 

3104 

1.11 

1.6 

1.0 

1.2 

+3-1 

+3.8 

3102 

1.23 

2.0 

0.9 

0.3 

+3.1 

+4.8 

Where 


■p  is  steady  state  roll  rate 

p  is  roll  rate  at  resonance 

*reson 

a  standard  error  of  fit 


£dM  = 

= 


100  1 


' 

1 

L 


100  1 


] 


^ef  8 
M0RDVAC 

^Ref  8 

TT 

ORDVAC 


Kf 

€ 


100 


1  _  Mg  Ref  8 

***  ORDVAC , 
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NON- LINEAR  AERODYNAMIC  FORCE  AND  MOMENT 


Although  non-linearities  can  arise  from  geometric  terms  as  well 
as  all  the  aerodynamic  terms,  we  will  confine  ourselves  to  angles  of 
attack  for  which  the  geometric  linearizations  are  correct  and  only  the 
static  and  Magnus  moments  are  non-linear.  The  more  complete  problem 
is  treated  in  References  10-11.  We  will  also  limit  ourselves  to  those 
non-linearities  which  arise  from  the  angles,  a,  B,  and  can  occur  for 
bodies  of  revolution  .  (For  configurations  with  less  symmetry,  terms 
other  than  these  can  occur.)  Under  these  assumptions  our  basic  equation 
assumes  the  form: 

1"  +  (H  -  iP)  -  £m(52)  +  iPT  (62)J  T=  0  (l6) 

„2  i^i  2  ~2  ,~<2 

where  6  =  Jfj|  =  a  +  P 

Range  tests  of  missiles  which  were  known  to  have  non-linear 
aerodynamic  coefficients  did  not  reveal  pitching  and  yawing  motions  which 
were  qualitatively  different  from  those  associated  with  linear  forces 
and  moments.  In  other  words,  it  was  not  possible  to  identify  non- 
linearities  from  the  angular  motion  of  individual  rounds  over  the 
relatively  short  observed  trajectories.  On  comparing  fitted  damping 
exponents  and  frequencies  from  different  firings  of  the  same  configuration 
at  the  same  Mach  number  it  was  found  that  they  varied  with  the  amplitude 
of  the  motion.  Hence,  the  non-linear  motion  over  short  trajectories 
can  be  approximated  by  the  linear  epicycles  but  the  epicyclic  parameters 
are  functions  of  some" effective"  magnitude  of  the  motion. 


In  order  to  obtain  quantitative  relations  we  will  make  use  of  a 

perturbation  process  which  1b  a  generalization  of  the  Kryloff-Bogoliuboff 
12 

method  and  is  similar  to  the  method  of  variation  of  parameters.  Eq. 
(16)  can  be  written  in  the  form: 

T"  -  iPT*  -  M(0)T  =  -  hT*  +  Tm(B2)  -  M(o)  +  iFT  (52)| 

L  J  (17) 
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If  the  right  side  of  Eq.  (17)  is  neglected,  then  for  a  gyro- 
scopically  stable  missile  the  solution  is  an  epicycle  without  damping: 


7=  K  e1^  +  K  e1^2 
5  10  20 


(18) 


where  s* 

-  0jo  +  » 


0'  =  \  [p  +  ~\j  -  4M(0)  J  and 


P2  -  k  M(0)  >  0, 


We  now  assume  a  solution  of  the  form 


'%  =  K±  e  1^1  +  *7  +  K2  e  1^2  +  *2^ 


(19) 


where  K  (s)  and  (s)  are  unknown  functions  of  s.  Substituting  Eq.  (19) 

J  J 

in  the  left  side  of  Eq.  (17),  and  inserting  a  condition  similar  to  Eq.  (9), 
we  can  now  write  equations  which  are  formally  equivalent  to  Eqs.  ( 13-14-) 
for  Aj  =  Kj  ell|rJ  and  =  0. 

to  +  i«  K.  )  e1^  =  -E  ^  1  i  (20) 

i 

to  ♦  i«  K,)  ^*2  .  -Hi'  t[M(6A2)^M(0)  ±  1FT(S2)]  1  (2l) 

■  K0£  -  0{) 


We  now  substitute  Eq.  (l9)  on  the  right  side  of  Eq.  (20)  and  divide 


by  K±  e 


i+l. 


K, 


-E  0^ 


+  i^  = 


1  + 


0 


^  e^  +  PT(52)  -  i[M(52)-M(0)|  1+^ 

*  K  L  J  l 

11 

-*•  -*• _ •>  W  _ M  «► _ 


,*0 


0[  -  0g 


(22) 
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where 


0  =  P2  +  +2  -  01  -  ti  and 

&2=  III2  =  f^2  +  Kg2  +  2Kj_  Kg  cos  0. 


Since  effective  values  of  the  damping  exponents  and  frequencies 
as  desired,  ve  average  Eq.  (22)  and  separate  into  real  and  imaginary 


Symmetric  expressions  apply  for  the  other  mode  of  oscillation. 


As  an  application  of  this  technique  we  now  consider  the  case  of  a 
cubic  static  moment.  The  static  moment  coefficient  is  quadratic  and  has 
the  form: 

CM  "  ac  +  *2  ^  («) 

Of 

H(52)  =  *;2  [*0  +  a2  B2J  (26) 

Substituting  Eq.  (26)  in  Eqs.  (24-)  and  its  partner  for  the  other  mode, 
we  can  obtain  the  following  relations. 


Since  T  and  M  are  even  functions  of  0,  the  sine  integrals  vanish. 
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(28) 


^  *t2  a2  (K22  *  2KI2> 

%  -  A 

Multiplying  Eq.  (2j)  by  Eq.  (28),  neglecting  fourth  powers  of  the  K's, 

'  J 

and  simplifying,  we  can  write  our  final  relation. 


C  •  t'2  =  aQ  +  a2 


where 


o  1 

c  m  k2  _  y 

S5  1 7 

2  2  2  ^  01 

5e  =  ^  +  K2  +  - 


K22  ^ 


01  -  01 


(29) 


If  we  call  the  left  side  of  Eq.  (29)  the  "effective"  static  moment 
coefficient,  this  equation  states  that  the  pairs  of  effective  static 
moment  coefficients  and  effective  angles  of  attack  squared  as  obtained  from 
different  tests  of  the  same  configurations  should  fall  on  a  line.  The 
slope  of  this  line  is  a 2  and  the  intercept  is  aQ.  Figure  12  shows  such 
a  plot  for  three  center  of  mass  locations .  Since  we  sould  expect  the 
coefficients  aQ,  a 2  in  the  moment  expansion  to  be  linear  in  c .m.  position, 
these  values  are  plotted  against  c.m.  position  in  Figure  13.  Good  internal 
consistency  can  be  seen.  The  slopes  of  the  lines  in  Figure  13  are  the 
corresponding  normal  force  coefficient,  i.e. 


CN  =bo+b2  5 

a 

These  values  and  others  have  been  compared  with  wind  tunnel 
measurements  with  good  correlations'^. 


(30) 


Since  the  maximum  angles  in  Figure  12  are  about  10°,  we  consider 
a  large  angle  of  attack  program  as  our  final  example.  (Fig.  14)  Here 
the  maximum  angle  of  attack  is  35°  and  a  single  line  does  not  fit  the 
data*.  A  pair  of  lines  form  a  good  approximation i  Since  a  line  in  this 


The  geometric  non-linearities  could  not  be  neglected  for  this  program 
and  the  appropriate  formulas  from  Reference  10  ha^e  been  used. 
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plot  corresponds  to  a  cubic  in  the  moment  plane,  this  means  the  moment 
is  approximated  by  a  pair  of  cubics  (Fig.  15)*  These  curves  have  been 
found  to  be  in  good  agreement.  The  discontinuous  derivative  can  be 
explained  by  a  shadowgraph  of  the  projectile.  At  about  22°  flow  separates 
on  the  lee  side. (Fig.  l6 ) 

In  addition  to  a  cubic  static  moment  this  equivalent  linearization 
method  has  been  successfully  used  for  the  measurement  of  cubic  Magnus 
moments  as  well  as  cubic  normal  and  Magnus  forces^.  As  we  can  see  the 
limited  concept  of  the  aerodynamics  range  as  a  tool  for  the  measurement 
of  linear  forces  and  moments  is  no  longer  correct. 


6k 


SUMMARY 


In  summary  it  should  be  noted  that  the  classical  use  of 
aerodynamic  ranges  has  been  extended  in  three  directions: 

(1)  Gun  boosted  rockets  have  been  successfully  handled 
by  the  range  technique . 

(2)  The  pitching  and  yaving  motion  of  an  symmetric  missile 
with  roll  rate  varying  through  resonance  has  been 
properly  analysed  by  use  of  both  analog  and  digital 
computing  machines. 

(3)  Accurate  measurements  have  been  made  of  cubic  static 
and  Magnus  moments  together  with  their  associated 
forces . 
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PLANE  NORMAL  TO  FLIGHT  PATH 


NON-ROLLING  STATICALLY  STABLE  SYMMETRIC  MISSILE 


70 


GENERAL  MOTION  PLANAR  MOTION 
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Figure 
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Figure  10 
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Figure  l6 
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SURVEY,  CALIBRATION,  AND  REDUCTION  TECHNIQUES  USED 
AT  THE  THOMPSON  AEROBALLISTICS  LABORATORY 

Wallace  H.  Allan 
Eldon  L.  Dunn 

U.  S.  Naval  Ordnance  Test  Station 

INTRODUCTION 

An  essential  requirement  for  free-flight  aeroballistics  ranges  is  an 
accurate  means  of  determining  the  orientation  and  position  of  the  missile. 

This  is  usually  done  photographically;  it  has  been  done  by  yaw  cards  and 
ideally  we  would  like  to  do  it  electronically.  For  the  present,  photographic 
systems  are  the  most  popular.  This  paper  is  concerned  with  the  photogrammetric 
problem  of  reducing  data  from'  the  Thompson  Aeroballistics  Laboratory.  The 
paper  describes  the  precise  position  survey,  the  photographic  calibration  of 
the  cameras  to  furnish  transformation  equations,  the  use  of  the  transformation 
equations  to  yield  position  and  orientation  of  the  missile,  alternate  solu¬ 
tions  for  position  and  orientation,  solutions  for  roll  orientation,  and  the 
results  of  the  reduction  as  applied  to  several  typical  rounds.  The  basic 
reduction  methods  were  developed  by  N.  S.  Hall,  formerly  of  the  Naval 
Ordnance  Test  Station. 

Photographic  methods  fall  into  two  classes;  those  used  on  small  ranges 
where  the  image  is  full  size  or  larger,  and  those  used  on  the  larger  ranges 
that  work  with  greatly  reduced  image  sizes.  The  scale  reduction  used  in  the 
Thompson  Aeroballistics  Laboratory  is  l/hO,  Most  ranges,  small  and  large, 
use  external  reference  systems.  In  these  systems  the  position  of  the  missile 
is  determined  by  using  reference  markers  in  the  field  of  view  of  the  camera. 

We  use  an  internal  reference  system  in  which  two  fiducial  markers  inside 
the  camera  are  used  to  pinpoint  the  location  of  the  missile.  The  accuracy 
of  the  external  reference  system  depends  upon  the  long-term  stability 
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of  the  building  or  tbs  reference  cables.  Tbs  interna.!  reference  system  depends 
upon  well  built,  temperature-compensated  cameras  and  stable  camera  piers.  Once 
tie  calibration  is  completed,  the  cameras  are  expected  to  maintain  their  posi¬ 
tions  and  orientations  for  at  least  a  year.  Figure  la  illustrates  the  camera 
and  its  mount-  The  camera  is  mounted  on  axes  that  intersect  at  the  front  nodal 
or  perspective  point  of  the  lens.  The  camera  may  be  swung  in  azimuth  or  eleva¬ 
tion  angles  without  disturbing  the  location  of  the  front  nodal  point.  The 
camera  may  be  oriented  and  securely  locked  in  position.  The  pier  is  made  of 
reinforced  concrete  and  is  sunk  more  than  8  ft.  underground.  Figure  lb  is 
a  view  of  the  inside  of  the  camera;  the  two  fiducial  lamp  projectors  may  be 
seen  in  this  view-  The  six  projecting  pins  are  located  in  the  focal  plane 
of  the  lens  and  the  glass  photographic  plate  is  positioned  against  these  pins. 
The  lenses  and  fiducial  projectors  are  mounted  so  that  as  the  temperature  of 
the  camera  changes,  the  distance  between  the  two  fiducial  marks  and  the  foeal 
distance  change  as  if  the  camera  were  made  of  glass.  This  maintains  the 
same  scale  values,  provided  the  photographic  plates  are  left  in  the  cameras 
long  enough  to  assume  the  temperature  of  the  cameras. 

PRECISE  SURVEY 

Our  cameras  are  calibrated  once  a  year  by  photographing  some  30  reference 
beads  in  the  field  of  view  of  each  camera.  The  positions  of  these  beads  in 
space  are  accurately  determined  by  means  of  a  precise  survey,  and  the  positions 
of  the  bead  images  on  the  plate  are  used  in  a  least-squares  fit  yielding  the 
coefficients  of  the  transformation  equations  for  each  camera.  Figure  2  is  a 
transverse  section  of  the  laboratory  showing  a  typical  camera  station.  The 
beads  are  fixed  at  two-foot  intervals  along  each  of  3  calibration  wires. 

These  wires  hang  in  a  tight  catenary  carve  near  the  range  line  within  the 
circle  entitled  "circle  of  missile  model  dispersion,"  in  Fig.  2.  The  vires 
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are  suspended  from,  towers  located  at  each  end  of  the  range.  Figure  3  is  a 
view  of  the  laboratory  looking  downrange.  The  downrange  tower  may  be  seen  at 
the  extreme  end  of  the  range.  The  catenary  wires  are  made  of  .030- inch-diameter 
Ni-Span-C  alloy,  an  invar-type  material.  The  wires  are  suspended  so  that  their 
effective  coefficient  of  expansion  is  approximately  l/4o  of  that  of  invar. 

In  order  to  calibrate  the  cameras,  the  locations  of  the  beads  and  cameras 
must  be  known  to  an  accuracy  of  at  least  a  .001  of  a  foot.  (Approximately  a 
l/lOO  of  an  inch. )  Since  there  are  bb  cameras  and  more  than  TOO  beads  used  in 
this  system  the  task  is  not  a  simple  one .  The  positions  of  the  beads  and 
cameras  are  found  in  a  coordinate  system  in  which  the  X-axis  is  parallel  to 
the  range  line  and  located  1  foot  to  the  left  (or  west)  of  the  left  line  of 
cameras  as  seen  in  Fig.  3*  Hie  Y-axis  is  cross-range  and  located  40  feet  up- 
range  from  the  first  camera  station  and  the  Z-axis  is  vertical  with  its  origin 
1  foot  underground.  In  laboratory  coordinates,  X  is  positive  downrange,  Y  is 
positive  to  the  right,  and  Z  is  positive  upward.  Thus  it  is  seen  that  this 
coordinate  system  is  lefthanded  and,  barring  some  catastrophe,  all  coordinates 
of  each  object  of  interest  are  positive. 

Our  survey  is  accomplished  with  Swiss  surveying  instruments.  The  Wild 
N-III  level  is  used  to  determine  the  level  of  each  camera.  Figure  4  illustrates 
this  instrument.  All  of  our  camera  elevations  are  nearly  the  same,  so  the 
same  rod  reading  is  used  between  fore  and  back  sights.  The  elevation  differences 
are  measured,  using  the  optical  micrometer.  This  instrument  can  give  differences 
good  to  a  .001  of  an  inch  between  cameras  20  feet  apart.  The  elevation  sights 
are  taken  on  a  small  rod  or  scale  resting  on  a  bench  mark  mounted  on  the  camera. 
Figure  5  illustrates  the  rod  and  bench  mark.  The  reference  mark  on  the  bench 
mark  is  located  directly  above  the  front  nodal  point  of  the  camera  at  a  fixed 
distance  from  the  nodal  point.  There  are  four  such  bench  marks,  any  one  of  which 
may  be  placed  on  a  camera  and  when  properly  leveled  will  indicate  the  same 
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location  within  several  .001  of  an  inch.  Figure  6  illustrates  an  alternate 
method  of  taking  levels  using  a  hydraulic  micrometer  level. 

The  distance  between  cameras  is  determined  with  an  invar  tape.  We  use  a 
20-foot  tape  for  measuring  distances  between  cameras  and  a  100-foot  tape  as  a 
primary  standard  against  which  the  20-foot  tape  is  compared.  The  two  tapes 
agree  with  each  other  over  the  100-foot  distance  within  a  few  .001  of  an  inch. 
Both  were  standardized  at  the  Bureau  of  Standards.  Figure  7  illustrates  one 
end  of  our  taping  rig.  The  tape  is  tensioned  by  means  of  a  20  lb.  dead  weight 
and  the  tape  read  with  the  magnifying  glass.  Figure  8  shows  the  other  end  of 
the  tape.  The  spring  scale  is  used  as  a  shock  absorber  to  prevent  overstraining 
the  tape  while  setting  the  rig  up.  Only  the  west  or  left  cameras  are  taped  for 
distance.  The  east  cameras  and  all  the  beads  are  located  by  means  of  theodolite 
triangulation.  Figure  9  shows  the  Wild  T-3  theodolite  mounted  on  the  basic 
bench  mark.  The  bench  mark  contains  a  center  hole,  and  a  stud  on  the  base  of 
the  theodolite  plate  is  positioned  in  the  hole.  The  theodolite  is  centered  on 
the  base  plate  by  means  of  a  dial  indicator  referenced  against  a  diameter 
concentric  with  the  inner  stud.  The  theodolite  may  be  centered  to  within  a 
.001  of  an  inch.  Figure  10  illustrates  the  Wild  T-3  and  an  operator  observing 
a  head  on  the  wire.  The  cameras  at  station  2  and  station  23  are  defined  as 
being  parallel  to  the  range  line  and  they  define  the  basic  camera  line.  All 
other  cameras  in  this  line  are  measured  as  to  their  cross-range  coordinates  by 
the  use  of  the  theodolite. 

The  Wild  level  is  used  to  give  the  z  or  elevation  coordinate  of  the  cameras. 
The  tape  gives  the  x  or  downrange  coordinate.  The  theodolite  gives  the  y  or 
cross-range  coordinate  of  the  cameras.  Each  camera  in  the  west  line  may  be  used 
as  a  theodolite  station  to  observe  the  beads  or  the  east  line  of  cameras. 

Azimuth  angles  taken  to  the  east  cameras  serve  to  fix  the  xy  coordinates  of  the 
cameras.  The  basic  equations  used  in  determining  the  xy  coordinates  from  the 
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azimuth  angles  are  illustrated  in  Fig.  11.  It  my  be  seen  from  Fig.  11  that 

the  position  of  an  unknown  such  as  Bead  x  y  is  determined  by  measuring  the 

3  3 

angles  a  and  £  ,  using  the  basic  camera  line  as  a  reference  line.  The 
coordinates  of  the  theodolites  are  the  same  as  those  of  the  cameras  upon 
which  they  are  mounted.  The  two  equations  have  two  unknowns  x^y^  ^or  which 
they  can  be  solved.  The  basic  equations  for  determining  the  z  coordinate  are 
illustrated  in  Fig.  12.  The  elevations  of  the  beads  may  be  measured  by  means 
of  the  Wild  T-3  theodolite.  An  elevation  angle  0  to  the  bead  may  be  measured 
and  the  elevation  calculated  using  the  elevation  z ^  of  the  theodolite  and  the 
horizontal  distance  D  from  the  bead  to  the  theodolite.  The  elevation  of  the 
theodolite  is  determined  by  finding  its  height  L  over  the  elevation  of  the 
front  nodal  point  of  the  camera.  The  height  L  is  measured  by  finding  the 
elevation  z+;  then  L  =  z  -  z  .  z  is  found  by  observing,  with  the  theodolite, 
a  bench  mark  on  an  adjacent  station.  The  dip  angle  e  and  distance  D  to  the 
camera  may  be  used  to  compute  the  distance  h.  This  distance  plus  1,  the  height 


of  the  bench  mark  above  the  nodal  point  of  the  camera,  plus  the  elevation  of 

the  observed  camera  may  be  used  to  calculate  z,  .  Then  L  is  determined.  L 

t 

remains  constant  for  a  particular  theodolite  as  it  is  moved  from  camera  to 


camera. 

The  horizontal  distances  between  beads  on  the  wires  are  secured  by  a  process 
by  which  each  bead  is  fixed  by  triar.gulation  in  space.  Only  azimuth  angles  need 
to  be  taken  for  this  measurement.  The  distances  between  beads  are  determined 
by  taking  three  separate  measurements  several  days  apart  and  using  the  average 
of  the  three  measurements.  The  distances  are  believed  to  be  accurate  to  within 
a  few  .0001  of  a  foot.  As  these  distances  are  added  together  accumulated  errors 
can  become  large,  so  distances  are  checked  every  20  ft.  against  triangulated 
data  and  then  corrected  if  necessary.  It  is  believed  that  the  accumulated 
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distances  are  good  to  within  .001  of  a  foot. 

Die  towers  move  slightly  every  day  and  the  wire  changes  its  sag  slightly 
with  temperature  so  that  the  positions  of  all  "beads  need  to  "be  determined  by  a 
method  that  can  be  accomplished  very  quickly.  Ve  do  this  by  what  we  call  a 
speed  run.  In  the  speed  run  all  of  the  coordinates  of  the  beads  are  measured 
by  a  process  which  takes  about  an  hour  and  15  minutes.  The  speed  run  is 
accomplished  by  a  technique  that  requires  two  theodolites  and  two  three-man 
crews  plus  a  supervisor  who  operates  a  third  theodolite.  Advantage  is  taken 
of  the  fact  that  the  three  lines  are  catenaries.  If  the  end  beads  and  midpoint 
bead  position  of  each  wire  sure  known,  then  the  position  of  the  rest  of  the  beads 
may  be  calculated  provided  the  distance  between  beads  is  known.  Actually,  to 
increase  accuracy,  more  data  than  this  is  taken.  Since  the  xy  coordinates  of 
the  beads  change  slowly  only  the  end  beads  are  observed  in  azimuth  during  the 
speed  run  and  these  measurements  can  be  verified  against  determinations  taken 
hours  before  and  after  the  speed  run.  In  order  to  get  a  good  fix  on  elevation, 
twelve  beads  on  each  wire  are  observed.  The  elevation  readings  are  taken  from 
the  nearest  theodolite  station,  which  means  that  twelve  different  theodolite 
stations  are  occupied  during  the  speed  run.  The  third  theodolite  serves  only 
to  observe  sag  to  correct  the  elevation  readings  of  the  other  theodolites  for 
changes  in  sag  with  time.  This  correction  is  very  small  because  the  theodolite 
crews  Btart  work  at  the  end  points  and  finish  at  the  middle  of  the  wire  where 
the  change  in  sag  is  largest.  No  bead  changes  more  than  .002  of  a  foot  in  sag 
between  its  measurement  and  the  taking  of  the  calibration  plate  and  this  is 
accounted  for.  Now  the  positions  of  all  the  beads  may  be  calculated  knowing 
the  distances  between  beads,  the  coordinates  of  the  end  points,  and  the  elevations 
of  the  twelve  beads  on  each  wire.  As  soon  as  this  speed  run  is  completed,  a 
calibration  picture  is  taken  of  the  beads  by  the  precision  cameras.  The  plates 
are  used  to  find  the  transformation  equations  of  each  camera. 
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DETERMDf  ATI  ON  OF  TRANSFORMATION  EQUATIONS 

The  calibration  plates  are  measured  in  a  Mann  comparator  that  is  accurate 
to  a  few  microns.  The  positions  of  each,  of  the  two  fiducial  marks  and  of  the 
individual  bead  images  on  each  photographic  plate  are  determined  with  respect 
to  the  rectangular  coordinate  frame  of  the  comparator;  the  plates  do  not  need 
to  be  positioned  in  the  comparator  with  particular  care  as  far  as  the  angular 
orientation  is  concerned.  The  cooparator  measurements  are  made  by  either 
one  or  two  readers,  and  each  reader  makes  N  Independent  measurements  of  each 
point.  The  individual  "x”  and  ny"  readings  are  automatically  recorded  on  IBM 
cards,  along  with  the  necessary  identification  code.  These  cards  are  the  input 
data  to  an  IBM  701  or  704  electronic  calculator,  both  of  which  are  coded  to 
calculate  the  coefficients  of  the  transformation  equations  for  each  camera. 

From,  this  point  on,  any  arithmetic  mentioned  is  performed  by  these  high-speed 
computers. 

The  original.  N  readings  of  each  point  are  averaged  separately  for  each, 
reader.  The  orthogonal  axes  of  the  comparator  coordinate  system  are  then  ro¬ 
tated  and  translated  by  standard  procedures  so  that  the  new  x-axis  passes 
through  the  two  fiducial  marks,  with  the  new  origin  midway  between  the  fiducial 
marks.  Figure  13  shows  the  relationship  between  the  two  coordinate  systems. 

The  positive  direction  of  the  y-axls  is  upward  as  defined  by  the  image  of  the 
photographic  background  in  the  laboratory,  and  the  positive  direction  of  the 
x-axis  is  directed  toward  the  downrange  fiducial  mark.  It  is  seen  that  the 
plate  coordinate  system  is  "left-handed”  for  the  right  or  east  cameras  and 
"right-handed"  for  the  left  or  west  cameras.  This  rotation  and  translation 
of  the  comparator  coordinate  system  to  the  camera  plate  axis  is  performed 
separately  for  the  measurements  of  each  reader.  If  two  readings  have  been  made, 
the  two  sets  of  plate  coordinates  are  averaged. 
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The  bead  coordinates  are  used  to  determine  the  coefficients  of  transforma¬ 


tion  equations  relating  focal  plane  positions  with  positions  in  space.  The 
positions  in  space  are  all  located  in  two  space-planes  approximately  parallel 
to  the  focal  plane  of  the  camera  and  located  near  the  range  line  of  the 
laboratory.  Figure  l4  shows  the  space-planes,  which  are  inclined  at  45  degrees 
to  the  horizontal  and  are  exactly  parallel  to  the  x-axis  or  range  line.  The 
coordinates  of  the  beads  on  the  calibration  plates  are  used  to  give  a  least- 
squares  determination  of  the  coefficients  of  two  general  quadratic  functions 
to  fit  the  known  positions  of  points  in  the  space-plane.  The  points  in  the 
space-plane  are  determined  from  the  survey  by  projecting  the  bead  coordinates 
into  the  space-plane  as  seen  from  the  camera.  This  is  simply  done  by  forming 
the  equation  of  a  line  through  the  bead  and  nodal  point  of  the  camera  and  then 
solving  for  the  intersection  of  this  line  with  the  space-plane.  Figure  15  shows 
the  transformation  equations  for  one  camera  relating  focal-plane  positions 
with  space-plane  positions.  The  10  coefficients  are  determined  by  the  previously- 
mentioned  least-squares  fit.  These  equations  take  account  of  the  camera  lens 
distortion  and  remove  the  need  for  the  cameras  to  be  exactly  aligned  to  the 
assumed  geometry.  All  scale  factors  between  plate  coordinates  and  space  coor¬ 
dinates  are  likewise  accounted  for  by  the  coefficients.  If  the  readings  of  the 
beads  are  now  applied  to  these  equations  and  the  positions  of  the  beads  in  space 
calculated,  then  the  standard  deviation  of  the  bead  positions  as  computed  with 
respect  to  the  actual  positions  is  approximately  .001  of  a  foot. 

CENTER  OF  GRAVITY  REDUCTION 

When  a  missile  is  tested  in  the  range,  accurately-timed  light  flashes  of 
one-microsecond  duration  are  utilized  to  obtain  multiple- image  photographs  of 
the  missile  at  each  station.  Figure  16  shews  a  typical  data  record  for  one  of 
the  cameras.  The  fiducial  marks  visible  at  the  ends  of  the  plate,  in  addition 
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to  relating  the  missile  records  with  those  of  the  calibration,  also  provide  the 
means  for  correcting  for  the  difference  in  plate  temperatures  at  the  times  of 
reading  the  calibration  plate  and  the  data  plate.  The  correction  factor  is 
equal  to  the  distance  between  the  fiducial  marks  for  the  calibration  plate 
divided  by  the  same  distance  for  the  missile  record.  This  correction  factor 
is  computed  separately  for  each  data  plate,  and  is  applied  to  all  comparator 
readings  of  the  missile  images. 

As  in  the  case  of  the  calibration  plates,  each  photographic  record  is 
measured  on  the  Mann  comparator  by  one  or  two  readers,  each  reader  making  U 
independent  measurements  of  each  nose,  tail  and  fiducial  mark.  The  comparator 
coordinate  axes  are  translated  and  rotated  so  that  the  plate  coordinates  are 
defined  in  the  same  manner  as  for  the  calibration  plates. 

At  this  time  the  transformation  equations  are  introduced  to  convert  the 
plate  coordinates  of  each  point  to  its  x  and  z  coordinates  in  the  appropriate 
space-plane.  Figure  14  illustrates  the  space-planes  previously  mentioned; 
from  an  inspection  of  the  figure  it  is  apparent  that 

y^=  29.28k)  -  z^ 

y~2  =  z2  + 

where  the  prime  indicates  a  right  camera,  a  double  prime  indicates  a  left  camera, 
and  the  subscript  2  indicates  the  point  on  the  ray  from  the  camera  to  the  point 
under  observation.  Since  the  coordinates  of  the  front  nodal  point  of  each  camera 
lens  are  known  and  the  positions  of  the  points  in  the  space-plane  are  known  from 
the  transformation  equations,  then  two  points  on  each  ray  are  known,  and  the 
equations  of  the  rays  can  be  written.  These  equations  are  displayed  in  Fig.  17; 
here  the  subscript  1  indicates  the  camera  coordinates.  The  expressions  in  the 
denominators  are  the  direction  numbers  of  the  rays. 
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It  is  unlikely  that  the  two  rays  will  intersect  and,  hence,  the  most 
probable  point  must  be  chosen  to  represent  the  position  of  the  point  under 
observation.  It  is  felt  that  neither  ray  is  more  accurate;  hence  the  midpoint 
of  the  common  perpendicular  to  the  rays  is  chosen.  The  direction  numbers  of 
the  common  perpendicular  may  be  obtained  as  the  direction  numbers  of  the  vector 
product  of  the  two  rays.  These  numbers  are  given  by  the  expressions  in  Pig.  l8; 
this  figure  also  shows  two  equations  of  the  common  perpendicular.  In  these 
equations  the  x's,  y's,  and  z's  have  the  property  that  they  lie  both  on  the 
perpendicular  and  on  the  ray.  Since  no  points  on  the  common  perpendicular 
are  known,  an  arbitrary  point  X^Y^Z^  on  this  line  is  chosen  to  write  its  equation; 
this  equation  and  the  equations  of  the  rays  my  then  be  solved  simultaneously 
to  yield  the  points  of  intersection.  The  fixed  point  is  chosen  in  such  a  way 

that  the  numbers  appearing  in  the  equations  are  of  reasonable  magnitude .  To 

\ 

do  this  it  is  sufficient  to  fix  the  coordinate  corresponding  to  the  largest 

direction  number.  As  the  common  perpendicular  is  in  all  cases  nearly  parallel 

to  the  x-axis,  the  direction  number  a  will  in  all  cases  be  largest.  Therefore, 

we  set  Xq  equal  to  20  times  the  station  number,  and  YQ  and  ZQ  can  be  determined 

in  terms  of  y1  or  yM  as  shown  in  Fig,  19«  Here  the  coefficients  are  functions 

of  the  two  sets  of  direction  numbers,  the  camera  coordinates,  and  X  and  hence 

0 

can  be  calculated . 

The  solutions  of  the  equations  of  the  rays  and  the  equations  of  the  common 
perpendicular  are  given  in  Fig.  20,  where  again  the  coefficients  are  the  same 
as  those  used  in  the  equations  of  Yq  and  Zq,  and  have  been  determined.  The 
final  values  of  x,y,  and  z  are  taken  as  the  average  of  the  x'  and  x",  y*  and  y* 
and  z1  and  zT1,  respectively. 

The  differences  between  x*  and  x",  etc.  are  a  direct  measure  of  the  miss 
distance  between  the  rays  from  opposite  cameras.  This  distance  is  usually  very 
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small;  however,  if  a  camera  has  changed  its  orientation  since  its  calibration, 
then  the  miss  distance  becomes  significant.  Thus  we  have  a  measure  every 
round  of  how  well  the  cameras  are  maintaining  their  calibration. 

Once  the  positions  of  the  nose  and  tail  have  been  determined,  the  direction 
cosines  of  the  missile  axis  may  be  calculated  by  standard  geometric  procedures. 
Then  the  position  of  the  center-of- gravity  of  the  missile  is  given  by  the 
expressions  in  Fig.  21.  In  these  equations  1  is  the  distance  from  the  CG 
to  the  nose  of  the  round  and  1*  is  the  distance  from  the  CG  to  the  base.  It 
is  seen  that  the  expression  for  each  coordinate  of  the  CG  is  the  average  of  the 
values  for  the  coordinate  based  on  the  nose  and  tail  points. 

CENTER  OF  GRAVITY  REDUCTION  FOR  BLUNT  BODIES 

Recently,  interest  has  arisen  in  determining  the  aerodynamic  characteristics 
of  blunt-nosed  re-entry  shapes  for  Intercontinental  Ballistic  Missiles.  Missiles 
having  these  shapes  are  not  susceptible  to  reduction  by  the  method  just  described, 
because  they  have  no  definite  "nose"  point  that  can  be  seen  from  all  perspectives 
available  to  the  cameras  in  the  laboratory.  However,  a  method  of  determining 
their  position  and  orientation  has  been  devised  and  will  now  be  described.  In 
case  the  base  of  the  missile  is  too  large  for  accurate  determination  of  its 
center,  a  slender  pin  is  installed  in  the  center  of  the  base  and  the  tip  of 
this  pin  is  read  as  the  base  of  the  missile;  at  the  same  time,  another  point 
on  the  axis  of  the  round  near  the  nose  is  read.  Figure  22  illustrates  such  a 
missile  and  the  points  read.  By  the  procedure  just  described,  the  coordinates 
of  the  base  are  determined,  as  well  as  the  coordinates  x^,  y^,  and  z^  in  the 
space-plane  of  the  points  on  the  rays  from  each  camera  to  each  of  the  two  points 
On  the  axis.  These  two  points  together  with  the  coordinates  of  the  front  nodal 
point  of  the  lens  of  each  camera  determine  a  plane  through  the  axis  of  the  round. 
There  are  two  of  these  planes,  one  for  each  camera,  and  they  Intersect  in  the 
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axis  of  the  round.  One  can  vise  a  solution  of  the  equations  of  the  planes 
analogous  to  the  vector  product  and  obtain  the  direction  numbers  and  hence 
the  direction  cosines  of  the  axis  of  the  missile.  Now,  knowing  the  direction 
cosines  of  the  missile  axis  and  the  location  of. the  base,  one  can  find  the 
location  of  the  CG  by  use  of  the  equations  given  in  Fig.  23,  where  d  is  the 
distance  from  the  CG  to  the  tip  of  the  pin. 

ROLL  ORIENTATION 

The  roll  orientation  of  the  missile  is  measured  with  the  aid  of  spin  pins 
fastened  in  the  base  of  the  missile.  These  pins  are  coded  by  means  of  their= 
shape;  one  pin  is  pointed,  and  the  other  pin  is  blunt.  These  pins  are  read  in 
the  comparator  and  their  position  found  in  space.  The  roll  orientation  of  the 
missile  is  found  by  finding  the  line  that  goes  through  the  position  data  for 
the  two  pins.  The  change  of  orientation  of  this  line  gives  a  measure  of  the 
spin. 


EXAMPLES 

Figure  24  illustrates  the  yawing  and  swerving  motions  of  a  finned  missile 
plotted  against  downrange  distance.  The  orientation  angles  (derived  from  the 
direction  cosines  by  use  of  the  small-angle  approximations)  were  converted  to 
yaw  angles  by  correcting  for  gravity  drop  and  drift  due  to  gun  alignment.  The 
damping  and  frequency  terms  for  this  round  were  well-determined  by  a  digital 
reduction.  Figure  25  shows  the  yawing  motion  of  a  typical  spinning  missile 
based  on  data  obtained  from  the  reduction  applied  to  a  5- Inch  projectile.  The 
data  were  subjected  to  an  iterative  epicyclic  yaw  reduction  on  a  digital 
computer  using  the  method  of  differential  corrections,  and  the  standard 
deviation  of  the  fit  was  about  3  milliradians . 

We  believe  that  by  the  application  of  the  survey  and  reduction  techniques 
described  in  this  paper,  we  are  able  to  achieve  excellent  accuracy  for  a  range 
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of  this  size.  One  of  the  principal  advantages  of  this  calibration  method  over 
some  of  the  methods  using  external  reference  systems  is  that  the  camera  is 
calibrated  over  its  entire  field  of  view  and  the  data  can  be  assessed  with 
equal  accuracy  throughout  the  entire  range  regardless  of  where  the  missile 
is  located  along  the  range  line.  This  permits  us  to  take  multiple -image  photo¬ 
graphs  with  each  camera,  thus  increasing  the  frequency  of  data  points.  Al¬ 
though  our  cameras  are  spaced  20  ft.  apart,  we  are  able  to  secure  accurate 
data  every  five  feet  along  the  trajectory.  This  feature  is  very  important 
to  us.  As  an  example,  the  curve  presented  in  Fig.  25  would  yield  a  less  accurate 
determination  of  turning  and  damping  rates  if  we  used  every  fourth  point  in 
the  fit. 

Thus  accuracy  and  frequent  data  points  are  one  of  our  principal  advan¬ 
tages.  However,  when  we  speak  of  accuracy  it  should  be  borne  in  mind  that 
we  are  referring  to  the  larger  type  of  aerobal  l.istic  range.  The  smaller 
spark  and  pressure  ranges  are  capable  of  obtaining  more  accurate  position  data, 
because  they  work  with  full  scale  or  larger  images.  The  principal  disadvantage 
of  our  range  is  that  we  are  not  able  to  achieve  good  orientation  data  on  short 
missiles.  Since  we  find  orientation  by  locating  the  position  of  the  nose  and 
tail  of  a  missile,  short  base  lines  reduce  the  accuracy  of  the  orientation 
data.  We  do  not  like  to  work  with  missiles  shorter  than  .5  of  a  foot. 

In  conclusion,  it  is  felt  that  the  survey,  calibration,  and  position 
reduction  techniques  just  described  render  the  Thompson  Aeroballistics  Laboratory 
a  well-equipped  tool  to  assess  the  aeroballistic  characteristics  of  models  or 
full  scale  missiles  ranging  from  1  to  3  feet  in  length. 
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FIG  DCWKRANGE  VIEW  OF  THOMPSON  AEROBALLISTICS  LABORATORY 


FIC  .  4  ViTLD  N-lll  LKVLL  US4D  FOR  Li;A2UTIKG  41417  AT  I  OK  OF  CAIJERA  BLriCH 

MARK. 


KF/45-69F20,  USI«TCTC,  Inyokem,  California 
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August  1953  ulC  UNIFIED 


FIG.  5  ELEVATION  ROD  ON  CAMERA  BENCH  MARX,  U3ED  FOR  READING  ELEVATIONS 

wEEH  THE  WILD  N-lll  LEVEL. 


NF/45-69821,  US  HOTS,  Inyokern,  California. 
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August  1953 


UNCLASSIFIED 


FIG..  5  HYDRAULIC  MICRCMETHt  LEVEL  USED  FOR  MEASURING  ELEVATION  OF 

CAMERA  BENCH  MARKS. 

HP/45-69822,  USMOTS,  Ir^okorn,  California  August  1?53  UHlIALSIb’XcD 
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FIG.  7  TAPING  STAND  EQUIPPED  WITH  DEAD  HEIGHT  TENSION  CONTROL,  USED  FOR  TAPING  DISTANCES 

BETWEEN  CAMERAS. 

NP/^-69C23,  USNOTS,  I«yoI:ern,  California  August  1953  UIELASSQ'BD 


FIG.  8  TAPING  STAND  EQUIPPED  WITH  SPRING  SHOCK  ABSORBER,  USED  FCR  TAPING  DISTANCES  BETWEEN 

CAMERAS. 

NP/45-69824,  us  ROTS,  Inyokern,  California  August  1953  UNCLASSIFIED 


FIG.  9.  WILD  T-3  THEODOLITE 
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;0N  EQUATIONS 


ELEVATION  PLAN 


FIG.  12.  ELEVATION  EQUATIONS 


FIG.  13.  TRANSFORMATION  0X3  PLATE  COORD: 
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FIG.  lk.  PROFILE  OF  LABORATORY  SHOWING  PROJECTION  PLANES 
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15.  TRANSFORMATION  EQUATIONS 


116 


FIG.  16.  TYPICAL  DATA  PLATE 
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PIG.  17.  EQUATIONS  OP  RAYS 
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FIG.  18.  EQUATIONS  OF  COMMON  PERPENDICULAR 
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FIG.  19.  SOLUTIONS  FOR  Y 
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2RSECTI0NS  OF  RATS  AND  COMMON 
PERPENDICULAR 
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FIG.  21.  POSITION  OF  CENTER-GF-GRAVITY 
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WAKE  VISUALIZATION  STUDIES  IN  THE  AEROBALLISTICS  RANGE 

G.  V,  Bull  and  G.  B.  Jeffery* 

Canadian  Armament  Research,  aid  Development 
Establishment 

1.0  INTRODUCTION 

In  the  computation  of  stability  and  control  characteristics  of 
wing-body- tail  configurations  such  as  aeroplanes  and  guided  missiles, 
knowledge  of  the  nature  of  vortex  shedding  is  essential.  Various 
techniques  have  been  used  to  date  to  study  vortex  shedding  phenomena 
in  wind  tunnels;  for  low  speed  work  smoke  tunnels  and  tuft-grid 
techniques  have  been  used  extensively;  in  the  high  speed  range  vapour 
screen  techniques  have  been  developed  by  the  NACA.  This  present 
paper  will  review  a  technique  presently  being  used  and  under  develop¬ 
ment  in  the  aeroballistics  range  of  the  Canadian  Armament  Research 
and  Development  Establishment. 

2.0  EXPERIMENTAL  TECHNIQUE 

The  experimental  technique  consists  of  introducing  a  laminar 
sheet  of  smoke  filaments  normal  to  the  line  of  fire  and  obtaining  by 
photographic  methods  a  time  history  of  filament  distortion  during 
and  after  the  passage  of  the  model.  Photographic  observation  has 
been  made  both  by  reflected  and  transmitted  light  techniques.  The 
rather  simple  experimental  set-up  is  shown  in  Figures  l(a)and  1(b) 

For  reflected  light  observations  a  microflash  unit  and  simple  camera 
were  employed;  due  to  the  high  speed  shutter  problem  encountered  this 
technique  was  limited  to  a  single  photograph.  A  high  speed  drum 
camera  was  tried  for  reflected  light  photography  but  was  temporarily 
abandoned  due  to  the  amount  of  development  of  the  technique  which 
appeared  necessary  for  successful  results .  An  on-axis  parallel 
beam  shadowgraph  system  was  then  set  up  using  a  multiple  timed 
spark  unit  as  a  light  source  (Figure  1(c)).  Geometrical  and 
diffraction  problems  with  this  system  led  to  its  modification  to 
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a  focussed  shadowgraph  system  using  lenses  focussed  on  the  plane 
of  the  smoke  to  form  an  image.  This  system  is  similar  to  the 
Cranz-Schardin  multiple  beam  system  and  is  shown  diagrammatically 
in  Figure  l(b).  The  spark  units  are  located  downstream  at  a 
distance  of  twice  the  fiacal  length  from  a  1 6"  diam  f/8  paraboloidal 
mirror.  This  mirror  is  centered  one  diameter  off  the  flight  path. 
The  reflected  beams  pass  through  the  plane  of  the  smoke  and  focus  at 
lenses  opposite  the  sparks.  The  lenses  are  focussed  on  the  plane  of 
the  smoke.  The  whole  apparatus  is  enclosed  in  a  light  tight  room 
so  that  shutters  are  not  needed.  This  system  is  essentially  a 
backlighted  normal  photography  technique  rather  than  a  shadowgraph 
technique,  and  is  found  to  produce  rather  high  quality  pictures. 

The  multiple  spark  unit  is  started  by  a  trigger  pulse  obtained  from 
an  infra-red  screen  located  just  downstream  of  the  smoke  and  there¬ 
after  the  sparks  automatically  discharge  consecutively,  with  time 
intervals  variable  as  desired  between  10  and  100  microseconds 
approximately .  A  reference  grid  of  wires  is  located  as  shown  in 
Figure  l(b)  to  allow,  by  simple  geometrical  calculations  an  accurate 
determination  of  model  position  with  respect  to  the  smoke.  Model 
velocity  is  measured  by  light  screens  in  the  normal  manner  and  the 
flight  trajectory  and  angular  orientation  are  measured  using  yaw 
cards . 


The  production  of  a  laminar  sheet  of  smoke  filaments  presents 
a  rather  difficult  problem.  In  the  test  procedures  employed  to 
date  it  is  necessary  for  the  smoke  generator  to  operate  for  intervals 
up  to  five  minutes.  The  photographic  quality  of  the  smoke  employed 
had  to  he  matched  to  the  particular  method  of  observation;  for 
reflected  light  work  ammonium  chloride  vapour  was  found  to  be 
photographically  suitable;  at  the  same  time  it  could  be  generated 
easily  and  gave  the  most  laminar  results  of  any  smoke  trial. 

The  simple  smoke  generator  developed  for  this  purpose  is  shown 
in  Figure  2.  In  this  unit  hydrochloric  acid  and  ammonium  hydroxide 
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are  allowed  to  mix  by  dripping  into  a  common  dish  and  the  resulting 
ammonium- chloride  vapour  is  ejected  at  very  low  velocity  through  ft 
series  of  glass  nozzles  by  wHvtng  with  pressurized  air  from  an  air 
bottle.  A  chimney  is  placed  above  the  generator  to  provide  an  updraft 
and  assist  filament  laminar ity.  With  this  unit  1  aminar  filaments  of 
vapour  10  to  12  inches  in  length  could  be  obtained  quite  readily. 

The  ammonium  chloride  vapour  did  not  prove  satisfactory  for  the 
backlighting  scheme  of  Figure  1-b,  and  it  was  found  necessary  to  use 
the  smoke  from  oily-rags,  ejected  through  an  identical  nozzles  system. 

The  experimental  set-up  described  above  has  oeen  used  to  study  the 
vortex- shedding  from  a  pitching  and  yawing  W&r  combination.  The  problem 
will  first  be  described  along  with  the  theoretical  background  work  and 
then  the  experimental  results  will  be  presented. 

I 

3.0  DESCRIPTION  OF  THE  PROBLEM 

The  specific  problem  studied  to  date  by  this  method  Is  associated 
with  the  determination  of  the  induced  velocity  field  (downwash  and 
sldewash)  behind  rectangular  plan-form  wings  located  in  cruciform  on 
a  cylindrical  after-body  with  a  conical  nose.  The  particular  configur¬ 
ation  studied  is  shown  with  geometrical  details  in  Figure  3. 

Linearized  aerodynamic  theory  has  been  applied  by  many  authors 
to  the  theoretical  study  of  this  problem  employing  an  assumed  form 
of  the  trailing  vortex  sheet*  Lagerstrom  and  Graham  (re|terence  l) 
used  conical  superposition  techniques  assuming  a  planar  undisturbed 
sheet  of  vorticity  extending  from  the  wing  trailing  edge  downstream 
to  infinity  in  the  streamwise  direction.  Other  authors  have  chosen 
a  fully-rolled  vortex  sheet  extending  from  various  positions  along  the 
wing-chord  downstream  to  infinity  in  the  streamwise  direction. 

Application  of  these  methods  to  a  variety  of  configurations  showB 
some  rather  large  discrepancies  with  experiment  primarily  associated 
with  the  manner  of  rolling  up  of  the  trailing  vortex  sheet.  Spreiter  and 


Sacks  (reference  2)  discussed  the  mutual  Interaction  effects  of  shed 
vorticlty  and  used  water-tank  studies  to  shew  experimentally  the 
manner  of  vortex  shedding  which  could  occur.  Rogers  (reference  3) 
applied  slender  body  theory  to  compute  the  distortion  of  the  wake 
behind  low  aspect  ratio  delta  wing-body  combinations ,  in  both 
subsonic  and  supersonic  flows.  In  the  present  studies,  the  method 
employed  by  Rogers  was  applied  to  cruciform  rectangular  wings  in 
supersonic  flows  for  which  slender  body  theory  could  be  assumed  to 
hold. 


4.0  THEORETICAL  WAKE  ROLLING  UP 

The  manner  of  rolling  up  of  the  trailing  vortex  sheet  behind 
lifting  surfaces  is  determined  by  the  spanwise  distribution  of  lift 
on  the  surface.  For  the  configuration  considered  here  -the  wing  panel 
loading  was  computed  by  using  strip  theory  and  accounting  for  body 
upwash  lift  cm  the  wings  by  the  method  of  Beskin  (reference  4).  Within 
the  framework  of  slender-body  theory,  changes  in  the  downstream  direction 
are  neglected  and  the  governing  linearized  differential  equation  becomes 
two-dimensional  in  form  and  independent  of  Mach  number. 

Under  these  assumptions  standard  line  vortex  theory  may  be  applied. 
The  lifting  surface  is  then  replaced  by  a  distribution  of  horseshoe 
vortices  leaving  the  surface  at  the  trailing  edge  and  passing  downstream 
in  the  streamnise  direction.  The  body  effect  is  taken  into  account  in 
the  normal  fashion  by  distributing  Inverse  vortices  corresponding  to 
the  external  field.  For  the  computational  work,  the  panel  loading  was 
approximated  by  using  ten  equal  strength  vortices  adjusted  slightly 
from  their  normal  position  to  give  a  smooth  curve  of  their  divided 
differences.  For  small  angles,  it  was  assumed  that  pitch  and  yaw  effects 
do  not  couple  through  the  wake,  and  therefore  from  the  symmetry  of  the 
model  it  is  necessary  only  to  compute  the  case  of  pure  pitch,  simple 
linear  superposition  yielding  the  combined  case.  The  induced  velocity 
at  each  vortex  position  due  to  the  complete  vortex  field  was  computed 
and  the  local  streamwise  direction  established.  Step  by  step  construction 
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of  the  vortex  paths  downstream  was  then  made  using  increments  in 
downstream  direction  of  approximately  .05  chord  lengths.  The 
numerical  computations  were  performed  under  contract  by  the  Computing 
Devices  of  Canada  Ltd.  on  a  digital  computer.  For  a  pure  pitch  of 
5°,  the  distortion  of  the  trailing  vortex  sheet  in  the  downstream 
direction  as  computed  is  shown  in  Figure  4  for  the  case  of  a  free 
stream  Mach  number  of  1-5*  From  this  figure  it  is  seen  that  the  wake 
tends  to  roll  into  two  vortex  cores,  one  inboard  and  one  outboard 
of  roughly  equal  strength. 

5.0  EXPERIMENTALLY  DETERMINED  WAKE  DISTORTION 

(a)  Reflected-Light  Technique 

Because  only  a  single  photograph  could  be  obtained  per  trial, 
this  method  was  not  extensively  used. 

A  typical  result  obtained  by  this  method  Is  presented  in  Figures 
5  (a)  and  (b)  which  show  two  views  of  the  model  of  Figure  3  fired  at  a 
Mach  number  of  1.5  and  having  a  pitch  angle  of  attack  of  l/2°  and 
zero  yaw.  The  side  view  (Figure  5  (a))  shows  the  position  of  the  model 
with  respect  to  the  smoke;  the  smoke  plane  Is  1.25  chord  lengths  back  of 
the  model  trailing  edge.  The  back  view  shows  the  vorticity  pattern  in 
the  smoke;  careful  inspection  will  show  an  outboard  translation  of  the 
filaments  above  the  wing  and  an  inboard  translation  of  those  below  the 
wing.  The  outboard  vortex  core  is  in  the  position  to  be  expected  from 
the  calculations;  the  inboard  vortex  cannot  be  detected.  Another 
interesting  result  obtained  from  a  test  round  similar  but  not  Identical 
to  that  of  Figure  3  is  shown  in  Figure  6.  This  projectile  had  a  pitch 
angle  of  attack  of  1°,  a  yaw  angle  of  attack  of  l/2°  and  was  fired  at 
a  nominal  Mach  number  of  1.5-  The  photograph  Is  taken  about  2  chord 
lengths  aft  of  the  trailing  edge.  The  vortex  formation  at  the  tips 
of  the  lifting  surfaces  Is  plainly  visible  with  no  visible  inboard 
vortex  formation  being  detectable.  There  also  appears  in  this  photo¬ 
graph  traces  of  the  intersection  of  the  shock  waves  from  the  leading 
edges  of  the  lifting  surfaces. 
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(b)  Back  Lighting  Results 


Because  of  the  simplicity  of  equipment  involved  and  the  time- 
sequence  of  wake  distortion  obtainable,  the  backlighting  scheme  has 
been  used  extensively  in  the  present  tests.  Figure  7  shows  the 
development  of  the  vorticity  pattern  in  the  smoke  field,  x/c  is  the 
distance  downstream  of  the  trailing  edge  in  wing  chords.  The  round 
shown  here  has  identical  geometry  to  that  shown  in  Figure  3  and  was 
fired  at  a  Mach  number  of  1.1*1  ;  the  model  developed  an  angle  of 
attack  in  pitch  of  1.1*°  and  in  yaw  of  0.25°.  The  reference  wire 
grid  is  visible  in  this  series.  In  interpreting  these  photographs, 
the  geometry  outlined  in  Figure  l(b)  should  be  studied  since  the 
model  appears  distorted  in  the  perspective  view  obtained. 

Figure  8  is  a  single  photograph  of  another  series  taken  with 
a  round  fired  at  a  nominal  Mach  number  of  1.1*1,  and  again  with  the 
geometry  depicted  in  Figure  3.  This  model  has  an  angle  of  attack  in 
pitch  of  2.25°  and  in  yaw  of  .6°,  and  is  rolled  20°  from  the  vertical 
(the  smoke  fl laments  are  flowing  vertically).  The  smoke  plane  is  3*3 
wing  chords  aft  of  the  wing  trailing  edge.  The  tip  vortex  cores  are 
visible  from  three  of  the  panels,  and  lost  from  the  fourth  in  the 
turbulent  region  of  the  wing  cut  in  the  smoke.  In  this  photograph 
refraction  effects  enable  visualization  of  the  vortex  cores;  it  can 
be  seen  readily  that  the  vortex  cores  in  the  tip  regions  leave  the 
wing  ahead  of  the  mid-chord  point.  Again  no  inner  vortex  cores  can 
he  detected. 


.0  ANALYSIS  AND  DISCUSSION  OF  RESULTS 

The  studies  described  in  this  report  are  still  in  progress;  to 
illustrate  the  kind  of  analysis  performed  on  the  experimental  data 
and  the  results  obtained  the  case  of  Figure  8  will  be  considered  here. 
Since  within  the  framework  of  the  slender  body  theory,  conditions  in 
any  plane  downstream  of  the  trailing  edge  may  be  computed  on  the  basis 


of  two-dimensional  Incompressible  theory,  the  theoretical  distortion 
of  the  smoke  filaments  can  be  computed  readily. 

From  Figure  8,  It  would  appear  that  most  of  the  vorticity  Is 
concentrated  in  the  tip  core;  two  sets  of  computations  were  performed, 
the  first  assuming  that  l/2  of  the  total  panel  vorticity  was 
concentrated  in  the  tip  core,  and  the  second  assuming  the  total  panel 
vorticity  to  be  concentrated  in  this  core.  Commencing  with  a  set  of 
planar  filaments  (stream  lines  of  negligible  velocity)  the  distortion 
of  these  stream  lines  due  to  the  two  different  strength  vortices  was 
computed  assuming  that  the  vortices  had  been  flowing  for  a  time  equal 
to  the  time  taken  for  the  model  to  travel  3.3  wing  chord  lengths  at  its 
velocity  of  15OO  feet  per  second.  The  resultant  streamline  patterns 
were  then  compared  to  that  obtained  by  measurement  from  Figure  8  correct¬ 
ing  for  linear  and  angular  distortions.  These  results  are  shown  in 
Figure  9  where  the  l/2  strength  and  full  strength  results  are  plotted 
to  the  left  and  right  respectively  of  the  experimental  results. 

From  these  computations  as  shown  in  Figure  9,  the  experimental 
field  corresponds  closely  to  the  resultant  field  assuming  a  total 
concentration  of  vorticity  in  the  tip  core.  The  strength  of  the 
inner  core  (See  Figure  4)  is  evidently  over-estimated  considerably 
by  the  theory.  Further  iterations  of  the  results  can  be  carried  out 
to  obtain  the  exact  vortex  strength  of  the  tip  core;  in  this  case 
the  petition  of  vortex  shedding  should  be  corrected  from  the  trailing 
edge  to  possible  the  mid-chord  point.  However,  for  conditions  far 
aft  of  the  trailing  edge,  such  as  represented  in  this  result,  the 
correction  is  small  (of  the  order  10$) .  From  the  analysis  of  the 
smoke  distortions  at  different  stations  downstream  of  the  trailing  edge, 
the  point  of  vortex  shedding  can  be  determined  accurately  as  well  as 
the  core  strength.  These  computations  are  being  conducted  for  the 
results  shown  in  Figure  7* 
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Some  interesting  observations  regarding  this  technique  can  be 
made  here.  By  theoretical  fitting  methods  the  observed  smoke  patterns 
can  be  used  to  determine  the  vorticity  distribution  in  the  vake  of 
lifting  surfaces  even  under  more  complex  cases  than  the  single  core 
problem  considered  here.  The  complexity  of  the  fitting  problem  increases 
with  the  number  of  cores  it  is  necessary  to  assume.  Once  having 
determined  the  vorticity  distribution  in  the  vake ,  the  inverse 
calculation  of  circulation  (lift)  on  the  surface  may  he  performed. 

This  technique  would  appear  to  he  peculiar  to  the  aerohallistics 
range,  since  it  requires  the  introduction  of  stream  filaments  in  a 
plane  perpendicular  to  the  air  flow. 

CONCLUDING  DISCUSSION 

The  analysis  of  experimental  results  already  obtained  is  being 
continued  at  CAKDE.  Other  experimental  techniques  are  being 
investigated:  among  these  two  may  he  mentioned.  The  first  is  the 
development  of  a  multiple  exposure  reflected  light  system,  since 
with  the  backlighting  method  blanketing  of  a  large  portion  of  the 
field  by  the  model  occurs.  This  system  will  be  inherently  more 
direct.  The  second  development  is  being  carried  on  by  A.  Easson 
to  develop  a  burning  wire  smoke  generator;  preliminary  results  from 
this  method  have  been  most  encouraging. 

ACKNOWLEDGEMENTS 

The  active  investigation  of  visualization  techniques  was  commenced 
after  a  discussion  with  Prof.  S.  M.  Bogdonoff  of  Princeton  University, 
who  suggested  that  such  a  line  of  attack  might  be  worthwhile . 

In  regard  to  the  experimental  work  the  authors  are  indebted  to 
Mr.  R.  Coderre  of  CARDE  for  his  enthusiastic  interest  in  the  problem 
and  in  particular  for  development  of  the  multiple  spark  system,  to  Lt. 

W.  Ready  for  conduct  of  the  numerous  test  firing  under  the  most  adverse 
conditions,  to  the  CARDE  photographic  section  in  general,  LAC  C.  MacKenzie 
in  particular. 


REFERENCES 


1.  Lagerstrom,  P.  and  Graham,  M. ,  Downwash  and  Sidewash  Induced 
by  Three-Dimensional  Lifting  Wings  in  Supersonic  Flow. 

2.  Spreiter,  John  R.  and  Sacks,  Alvin,  H. ,  The  Rolling  Up  of  the 
Trailing  Vortex  Sheet  and  its  Effect  on  the  Downwash  Behind 
Wings.  Jour.  Aero.  Sci.  Vol  18,  No.  1  -  Jan.  1951. 

Sacks,  A.  H. ,  Behaviour  of  Vortex  System  Behind  Cruciform 
Wings,  -  Motions  of  Fully  Rolled-up  Vortices  -  NACA  TN  2605, 
1952. 

3.  Rogers,  Arthur,  Application  of  Two-Dimensional  Vortex  Theory 
to  The  Prediction  of  Flow  Fields  Behind  Wings  of  Wing-Body 
Combinations  at  Subsonic  And  Supersonic  Speeds  NACA  TN  3227- 

4.  Beskin,  L.,  Determination  of  Upwash  Around  A  Body  of  Revolution 
at  Supersonic  Velocity  APL/jhU/CM  -  251. 


137 


SMOKE 

GENEPATOE 


(a.)  reflected  L/GHT  system 


(b)  Sack- i/<s nr/MG  system 

F/GUEE  /  /  E /PEE /MENTAL  SETUP  USED  F6Z 

Y&erEr  shedd/ng  studies 


138 


SMOKE  STATION 


159 


FIG  1(G) 


EXHAUST 


140 


COMPUTED  ACCORDING  TO  SLENDER  BODY 
ASSUMPTIONS  FOR  RECTANGULAR  WING-BODY 
COMBINATION 
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DISTORTION  OF  TRAILING  VORTEX  SHEET 


FIG.  5(a) 

REFLECTED  LIGHT  SIDE  VIEW  OF 
MODEL  8  SMOKE 
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FIG.  5  (b) 

REFLECTED  LIGHT  VIEW  OF 
SMOKE  DISTORTION 
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FIG.  6 

REFLECTED  LIGHT  PHOTO 
OF  SMOKE  SHOWING  SHOCK 
WAVE  INTERSECTION. 
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FIG.  7 

DEVELOPMENT  OF  VORTEX 
PATTERN  IN  SMOKE 
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FIGURE  8  :  SMOKE  PATERN 
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CIRCULATION  DUE  TO  S  VORTICES  PHOTOGRAPHIC  RESULT  CIRCULATION  DUE  TO  10  VORTICES 

(1/2  CIRCULATION)  (TOTAL  CIRCULATION) 


FIG.  9  COMPUTED  AND  EXPERIMENTAL  SMOKE  FILAMENT  DISTORTION 
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THE  GCEraiROLTM-TEMFERATUKE-PBESSUHE  RANGE: 

SPECIFICATION,  PROBLEMS  OF  MEASUREMENT, 

AND  RECENT  PROGRAMS 

by 
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Ballistic  Research  laboratories 
Aberdeen  Proving  Ground,  Maryland 

INTRODUCTION 

1.1  Historical  Survey 

The  development  program  which  has  culminated  in  the  high  Mach  number 
instrument  called  the  C ont rolled-Temperature -Pres sure  Range,  abbreviated 
CTPR,  had  its  first  beginnings  in  the  period  frcm  1944-1946.  During  this 
interval  Col.  Paul  Idbessart,  of  the  French  Army,  Drs.  R.  F.  Clippinger, 

I.  E.  Segal  and  R.  N.  Thomas  in  various  conferences  and  memoranda  evolved 
the  idea  of  an  enclosed  firing  range  in  which  the  working  fluid  could  be 
maintained  at  closely  controlled  pressure  and  temperature.  While  the 
earliest  record  is  in  the  form  of  a  typewritten  memo  dated  June  1944 
from  Col.  Libessart  to  the  late  Dr.  Edwin  P.  Hubble  who  was  then  director 
of  the  Exterior  Ballistics  laboratory,  the  final  design  is  the  conception 
of  its  principal  architect  Dr.  R.  F.  Clippinger.  Subsequently  many 
people  have  participated  in  the  construction,  modification  and  operation 
of  the  CTP  Range.  Some  of  their  contributions  appear  as  items  in  the 
bibliography  to  this  paper  and  will  be  indicated  by  bracketted  £  numbers 
in  the  text. 

From  the  beginning  it  was  clearly  recognized  that  control  of  temperature, 
particularly  through  use  of  refrigerants  such  as  liquefied  gases,  offered 
the  possibility  of  choosing  the  ambient  velocity  of  sound  so  as  to  obtain 
high  Mach  number  flows  about  freely  moving  supersonic  projectiles.  Through 
use  of  refined  optical  techniques  offered  by  the  Mach-Zehhder  interferometer 
it  was  hoped  that  comparisons  could  be  drawn  between  experimentally  produced 
and  theoretically  calculated  fields  of  flow. 
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The  seemingly  easier  objective  of  high  Mach  numbers  has  been  realized 
only  recently  while  the  apparently  more  complex  comparisons  of  experimental 
and  calculated  flows  [1,2]  could  be  carried  out  earlier.  Certain  simple 
cases  were  studied  in  a  prototype  Model  Range  with  pressure  control  only. 
This  prototype  which  is  now  in  its  Mark  II  design,  was  constricted  expressly 
to  provide  experience  with  an  enclosed,  pressurized  range.  Many  valuable 
results  have  been  obtained  from  it  and  lessons  have  been  learned  which  could 
be  applied  to  problems  of  instrumentation  for  the  CTPR .  The  Model  Range 
remains  a  useful  tool  for  studies  requiring  variable  density  only,  but  has 
been  superseded  by  the  CTER  where  high  Mach  numbers  are  desired. 


II  RANGE  DEVELOPMENT 


2.1  Control  of  Pressure  and  Temperature 

In  Figure  1  IS  given  an  artist’s  conception  of  the  CTFR  in  schematized 
form.  The  main  outlines  of  the  three  principal  parts,  viz.,  the  gun  chamber, 
central  working  section,  and  stopping  section  can  be  seen.  The  oil-diffusion 
vacuum  pumps,  and  the  pressure  and  nitrogen  lines  are  indicated  by  V,  P  and 
N  respectively.  Figure  2  shows  ah  enlarged  view  of  the  working  section  and 
stopping  tube.  The  dark  frame  in  the  center  represents  the  external  housing 
of  the  10"  Mach-Zehnder  Interferometer.  At  the  cutaway  section  may  be 
discerned  both  the  annular  evacuated  space  between  inner  and  outer  cylinders, 
and  the  cooling  and  heating  coils  wrapped  around  the  Inner  chamber. 

Both  the  type  and  density  of  the  working  fluid  can  be  controlled  In 
the  CTFR;  thus  7  and  may  be  chosen  as  desired  throughout  the  gamut  of 
available  gases  and  continuously  in  the  domain  of  pressures  0^  p  £  10  atm. 

For  most  room-temperature  experiments  dry  sir,  dew  point  t  k,  -70°F,  from 

,jS 

the  Wind  Tunnel  sphere  is  used  hoth  as  working  fluid  and  as  a  rinse  for  the 
gun  section  after  firings.  Since  the  combined  volumes  of  the  CTFR  sections 
do  not  exceed  600  cu.  ft.  and  ordinarily  only  about  one-third  of  this  volume 
is  changed  per  firing,  the  demands  on  the  sphere  supply  are  not  high.  At  the 
extreme  low  temperature  T  =  80°K,  only  permanent  gases  may  be  used.  The 
choice  so  far  has  been  dry  nitrogen,  mainly  because  of  its  similarity  to  air, 
secondarily  because  it  can  be  conveniently  introduced  from  the  piping 
system  carrying  the  refrigerant . 

Control  of  temperature  is  achieved  by  circulating  liquids  or  gases 
in  two  sets  of  codling  coils,  with  top  and  bottom  headers,  placed  as 
uniformly  as  possible  around  the  45’ >  stainless-steel,  inner  cylinder  of 
the  working  section. 

For  low  temperature  runs,  liquid  nitrogen  supplied  by  tank  truck  and 
centrifugally  pumped  through  one  set  of  coils  at  pressures  up  to  30  psi 
provides  a  minimum  temperature  near  85°K  over  the  entire  observation 
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Figure  3*  Schematic  Diagram  of  the  10"  Mach-Zehnder  interferometer.  Spark  source 
and  liquid- prism  monochromator  appear  at  the  upper  extreme  right.  One 
beam  passes  through  the  CTFR,  the  other  through  a  cylindrical  compensating 
chamber.  Components  may  be  identified  as  follows;  x  -  light  source; 

1,  9  -  lenses;  3,  6  -  divider  plates;  2,  k,  7 >  8-  plane  mirrors; 

10  -  rotating  mirror;  11  -  film  holder. 
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chamber.  Heat  gains  from  the  outside  shell  axe  minimized  by  maintaining 
the  intervening  annular  space  at  high  vacuum,  p  £  .05  micron,  and  insu¬ 
lating  the  inner  evacuated  wall  with  several  layers  of  aluminum  foil 
radiation-shield . 

Temperatures  higher  than  ambient  are  obtained  by  circulating  electri¬ 
cally  heated  oil  from  an  insulated  tank  through  the  other  set  of  coils. 

The  upper  temperature  limit  T  =  400°K  is  set  by  the  efficiency  of  the  oil 
heating  system  and  the  chemical  stability  of  the  oil. 

2.2  Instrumentation 

The  working  section  of  the  CTFR  has  five  sets  of  paired  crossed  windows 
with  6"  diameter  apertures  separated  by  9'  intervals.  Each  is  instrumented 
for  spark  shadowgraph  records  except  at  the  center  station.  Here  a  hori¬ 
zontal  window  pair  for  shadowgraph  and  triggering  beam  apparatus  precedes 
by  about  two  feet  a  vertical  station  with  10M  diameter,  optical  glass 
windows  accurately  aligned  for  the  10"  Mach-Zehnder  interferometer. 

Reference  marks  in  the  shadowgraph  stations  are  surveyed  to  an  accuracy 
of  .001'  on  the  centerline  and  the  stations  are  rendered  perpendicular  to 
the  centerline  to  within  a  few  seconds  of  arc  by  means  of  a  precision 
telescope  and  optical  square. 

Time  of  flight  intervals  between  shadowgraph  stations  are  measured 
to  within  +  0.1  microsecond  by  means  of  Hewlitt -Packard  10  Me  counters. 

The  10"  Mach-Zehnder  interferometer,  one  of  the  largest  of  its  kind 
used  for  full-field  flow  surveys,  is  equipped  with  an  exploding-wire 
light-source  [5]  ,  a  liquid-prism  monochromator  and  rotating - 

mirror  camera.  The  instrument  itself  is  arranged  as  shown  in  Figure  5 
with  its  10"  elements  at  the  vertices  of  a  60°  -  120°  parallelogram 
oriented  perpendicular  to  the  range  centerline .  Fringes  parallel  to  the 
centerline  fulfill  both  the  condition  for  optimum  light  source  £  and 
the  requirements  for  precision  measurement  of  fringe  shift  in  the  field 
of  flow  [2,  Sec  nr]  . 
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2.J  Projectile  Launchers 

The  absolute  velocity  attainable  in  a  free  flight  range  is  set  by  the 
types  of  projectile  launchers  available.  Until  recently  cal  .3$  and  cal 
.50  smoothbore  and  rifled  guns  were  the  principal  means  of  launching 
projectiles.  Full  size  spin-stabilized  missiles  with  driving  bands  and 
sub-caliber  missiles  with  separating  sabots  have  been  extensively  used  for 
experimental  firings.  With  cal  .JO  fixed  ammunition  and  special  powder 
loads  velocities  up  to  5000  fps  can  be  reliably  obtained.  With  the  cal  .50 
smoothbore  velocities  up  to  6000  fps  have  been  achieved  in  a  number  of 
programs . 

Quite  recently  a  cal  .75  smoothbore  gun  £6]]  with  detachable  chamber 
and  detached  ammunition  has  been  used  to  fire  sub-caliber,  non-spinning 
static ally-stable  rounds  in  separating  sabots .  Velocities  up  to  7000  fps 
at  the  muzzle  have  been  attained. 

2.4  Mach  and  Reynolds  Number  Bounds 

Since  a  principal  feature  of  the  CTPR  is  its  potential  to  utilize 
cooled  gases  in  which  the  velocity  of  sound  is  low,  one  expects  the  span 
of  accessible  Mach  numbers  to  be  considerably  enlarged.  Assuming  a 
perfect  gas  law  dependence  of  sound  velocity  on  absolute  temperature 
we  find  that  at  the  lowest  temperature  of  the  CTPR  Mach  number  can  be 
increased  by  a  factor  of  approximately  [3OO/85J  =1.9*  Thus  with 

the  cal  .75  projectile  launcher  the  highest  Mach  number  attainable  is 
slightly  less  than  12.  On  the  lover  side,  firings  below  M  =  1  have  been 
made  on  occasion,  but  routine  studies  much  below  M  =  1.2  are  complicated 
by  interference  from  reflected  shock  waves.  We  may  then  say  that  with 
its  present  instrumentation  the  Mach  number  scale  of  the  CTPR  satisfies 
the  inequality  1 . 2  £  M  £ 12 . 

Since  density  may  be  controlled  directly  in  the  CTPR,  it  is  possible 
to  make  Reynolds  number  Re  vary  over  a  considerable  range  of  values . 

With  standard  values  of  the  physical  constants  for  the  density,  viscosity 
and  molecular  weight  of  air  we  can  find  for  Re  =  (fi  VL/| 1)  the  value 
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Re  =  .55  x  10°  per  inch,  per  Mach  number  per  atmosphere  at  T  =  JOO°K. 

Since  at  constant  pressure  density  varies  inversely  with  temperature, 
this  specific  Reynolds  number  Re  increases  with  decreasing  temperature 
and  at  our  low  temperature  limit  would  be  larger  by  nearly  four  times . 

To  illustrate  the  possible  span  of  Reynolds  number  values,  we  note 
that  since  the  CTPR  can  be  evacuated  to  pressures  below  1mm  Hg,  the  lower 
bound  is  zero  and  may  be  approached  as  closely  as  we  like.  For  a  repre¬ 
sentative  projectile,  say  0.5"  in  diameter,  a  room  temperature  study 
could  be  carried  out  at  M  =  6  over  the  full  pressure  range  0£pf;10  atm. 
and  the  corresponding  Reynolds  numbers  would  be  O^Refrl-5  x  10^. 

Since  the  boiling  point  of  oxygen  is  90 °K,  air  cannot  be  used  at 
the  extreme  low  temperature.  Nitrogen  which  boils  about  77°K  has  been 
used  successfully  at  pressures  less  than  1  atm.  Here  the  maximum  Re  for 
our  example  would  be  smaller  by  a  factor  of  about  three . 

At  the  maximum  Mach  number  and  pressure,  a  room  temperature  study 

7 

could  be  made  at  about  Re  =  3  x  10  -  at  least  in  principle.  Actually 

at  densities  above  that  characteristic  of  air  at  room  temperature  i.e. 
pj  =  1.2  gm/i,  a  number  of  difficult  problems  of  measurement  appear  which 
severely  limit  the  accuracy  possible.  To  cite  a  few  of  these:  1)  a  high 
velocity  missile  decelerates  an  appreciable  fraction  of  its  velocity  in 
50'  of  travel  when  air  densities  in  excess  of  3Pi  are  used,  2)  shadowgraph 
and  interferometer  records  are  complicated  by  very  large  refraction  effects 
and  3)  refraction  may  cause  important  portions  of  the  optical  record  i.  e. 
large  bundles  of  light  rays,  to  miss  the  rotating  mirror  entirely  and 
fail  to  show  on  the  photographic  plate,  h)  stability  of  the  projectile 
if  spin  stabilized  may  be  adversely  affected. 

While  the  potentiality  tbr  specification  of  Re  through  a  considerable 
range  of  values  is  a  feature  of  the  CTHR,  practical  limitations,  which  can 
eventually  be  removed  or  minimized,  have  confined  variable  Reynolds  number 
studies  so  far  to  Mach  numbers  lees  than  six  (Mf:6)  and  densities  less 
than  3.2  gms/i  (p£3P1)- 
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It  is  clear  that  Reynolds  number  cannot  be  controlled  independently 
of  Mach  number  unless  compensating  changes  are  made  in  one  of  the  other 
variables  such  as  density,,  projectile  size  or  type  of  gas-  It  is  found 
experimentally  that  quantities  such  as  drag  coefficient  depend  only 
weakly  on  Re  over  considerable  ranges  of  M.  Thus  in  presenting  vs  M 
data  the  variation  of  Re  from  point  to  point  is  frequently  ignored. 


Ill  PROBLEMS  OF  MEASUREMENT 


3»1  Space  and  Time 

Measurements  of  projectile  center-of-mass  positions  relative  to  some 
reference  starting  point  and  of  the  corresponding  time  intervals  are  among 
the  simplest  and  most  fundamental  that  cnn  be  made  in  a  firing  range. 

From  these  may  be  derived  values  of  ballistic  drag  coefficient  K^. 

The  effect  on  drag  coefficient  of  random  errors  in  time  has  been 
studied,  for  ballistic  ranges,  by  Karpov  £7^  and,  for  range  windtunnel 
combinations,  by  Charters  £ 8 J  vho  applies  Karpov* s  results.  These  are 
based  on  the  following  assumptions:  l)  time  t  can  be  accurately  repre- 

"x'-v  i 

sented  as  a  function  of  distance  Z  by  a  cubic  polynomial  t  =  /  a^Z  , 

which  is  fitted  by  the  method  of  least  squares  to  the  set  w  $  of 

station  data,  2)  the  times  t,  are  subject  to  random  errors  of  measurement 

satisfying  independent  normal  distributions  each  of  which  has  the  Bame 

variance  cr^.  and  3)  the  range  has  an  odd  number  of  Stations.  The  effect 

of  random  errors  in  the  Zj  is  accounted  for  heuristically  by  supposing 

the  existence  of  an  equivalent  timing  error  a  =  -where  cr  is  the 

standard  deviation  of  Z,  assumed  the  same  for  each  station,  and  V  is  the 

average  velocity  of  the  projectile.  The  total  variance  in  time  is  then 
2  2  2 

taken  to  be  0^,  =  ( cr^  +  cre ) . 

Karpov  then  finds  that  the  principal  source  of  error  in  K^,  aside 
from  constants  of  the  projectile  and  atmosphere,  arises  from  the  quadratic 
coefficient  a^-  With  the  number  of  stations  fixed  this  error  varies 
inversely  as  the  square  of  range  length. 

The  significance  of  this  result  for  the  CTPR  may  be  seen  from  the 
following  consideration.  In  a  study  of  the  bounds  on  errors  In 
encountered  in  the  five  station,  20*  Model  Range,  the  author  and  J.  M. 
Bartos  PJ  established  under  representative  assumptions  the  maximum 
limits  j4.6$  and  +0.9$  for  firings  of  a  blunt  cone-cylinder  at  pressures 
of  1  and  5  atmospheres  respectively.  Since  the  base  line  of  the  CTFR  Is 
long  and  the  number  of  stations  5,  for  the  sSme  projectile  will  be 
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more  accurately  determined  by  a  factor  of  more  than  4.  Since  that  writing, 
timing  accuracy  has  "been  increased  by  a  factor  of  6  through  adoption  of 
10  me.  counters.  Thus  errors  of  the  order  of  1  $  or  less  in  may  be 
regarded  as  optimum  for  the  CTFR.  In  point  of  fact  unavoidable  uncertainties 
in  the  gas  density  p1  which  enters  the  calculation  are  sometimes  severe 
enough  to  mask  al 1  other  sources  of  error.  This  type  of  error  is  important 
only  at  temperatures  near  the  extreme  lower  limit,  85°K.  Temperature 
gradients  in  the  Range  and  uncertainties  in  pressure  and  temperature  de¬ 
terminations  combine  to  contribute  errors  in  excess  of  1$  but  probably 
less  than  3$  for  representative  high-drag  projectiles. 

Two  points  regarding  the  error  study  remain  to  be  mentioned.  A 
recent  re-examination  of  the  problem  by  G.  D.  Kahl  and  the  author  dhows 
with  range  length  constant  the  error  in  a^  diminishes  as  the  inverse 
square  root  of  number  of  stations  when  this  number  is  large.  Furthermore 
the  equivalent  timing  error  og  attributed  to  errors  in  Z  is  Just  that 
obtained  by  assuming  that  the  normal  equations  define  the  coefficients 
as  functions  of  t  and  Z.  With  this  assumption  the  variances  of  the  set 

may  be  now  calculated  directly  in  terms  of  the  rats  errors  cr^  and 
Thus  while  the  prospect  of  minimizing  squared  error  with  respect  to  both 
t  and  Z  is  explicitly  given  up,  the  introduction  of  the  equivalent  timing 
error  to  account  for  distance  error,  increases  the  variances  In  a  manner 
which  is  both  strictly  in  accord  with  the  most  convenient  assumption  of 
functional  relation,  between  t  and  Z,  and  in  agreement  with  intuition  as 
to  the  likely  ratio  of  errors  in  t  and  Z. 

3,2  Optical  Effects 

A  minute  examination  of  errors  in  time  measurement  for  the  CTPR  is 
relatively  unprofitable  for  the  reason  that  with  the  Introduction  of 
counters  capable  of  discriminating  0.1  Usee  Intervals,  the  rms  error 
is  below  the  equivalent  timing  error  cr^/v  by  a  factor  of  five  for 
V  =  10,000  ft  per  sec.  Thud  improvement  in  distance  measurement  has 
a  clear  priority  at  the  present  time.  One  of  the  major  Sources  of 
error  in  distance  measurement  is  uncertainty  inherent  in  the  optical 
methods  employed.  This  topic  will  be  discussed  under  three  headings,  viz.. 
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l)  Ideal  geometrical  shadowgraphs,  2)  refraction  effects,  and  3)  dif¬ 
fraction  effects. 

In  the  CTER,  measurement  of  distance  reduces  in  essence  to  determination 
of  the  interval  on  a  photographic  film  or  plate  between  the  shadowgraph 
Jjfljtges  of  a  fiducial  mark  and  the  projectile. 

An  ideal  shadowgraph  consists  of  rectilinear,  conical  projection  of 
the  projectile  Image  Onto  a  plane.  Under  these  conditions  a  perfectly 
sharp,  enlarged  image  of  the  body  is  obtained. 

If  geometrical  optics  is  assumed  and  the  refracting  medium  is  uniform, 
then  only  the  size  of  the  source  would  affect  image  sharpness.  If  we  call 
the  source  diameter  S,  the  distance  along  the  normal  to  the  film  frm  source 
to  projectile  d^  and  from  projectile  to  film  d^,  then  the  angle  subtended  by 
source  at  projectile  is  a  =  S/d^  and  the  width  of  penumbra  around  the 
projectile  is  W  =  adg  =  Sd^/d^-  The  magnification  of  the  shadowgraph  is 
m  =  (cLj  +  d^ ) / d^ .  por  the  CTFR,  Sparks  with  S  =  .05"  are  in  routine  use 
with  cLj  =  5 6"  and  =  32" .  These  dimensions  are  fixed  by  the  size  of  the 
outer  shell  of  the  range.  Thus  the  penumbra!  region  of  partial  shadow  and 
consequent  edge  uncertainty  will  be  .027"  on  an  image  about  twice  normal 
size  (m  =  1.9).  Estimation  of  the  shadow  center  can  reduce  the  edge  un¬ 
certainty  perhaps  by  a  factor  of  2,  but  cannot  eliminate  it  entirely. 

Smaller  magnifications  and  sharper  shadowgraphs  can  be  obtained  by  using 
a  camera  lens  to  focus  a  plane  nearer  the  projectile  onto  the  photographic 
plate.  These  "focussed  shadowgraphs"  are  preferred  whenever  possible  for 
these  reasons  and  others  to  be  explained  presently.  Nevertheless  Instal¬ 
lation  at  each  station  of  a  large  diameter  lens  of  good  quality  is  a 
considerable  complication  and  introduces  uncertainties  of  its  own. 

From  the  standpoint  of  flow  studies  the  most  valuable  portion  of  a 
shadowgraph  arises  not  from  its  approximation  to  an  accurate  conical 
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projection,  but  from  the  deviations  from  this  ideal  caused  by  inhomogeneities 
in  the  refracting  medium. 

Conventional  treatments  of  these  refraction  effects  [n,  12]  start 
from  the  assertion  that  the  curvature  of  a  light  ray  is  equal  to  the 
magnitude  of  the  vector  product  of  the  ray  tangent  with  the  gradient  of 
refractive  index  j"bf .  lj].  In  a  one -dimensional  situation  with  a  weakly 
refracting  medium  like  air,  |Vn  |  =  dn/dx  and  rays  crossing  transversely 
a  region  of  constant  gradient  will  become  and.  remain  circular  arcs  so  long 
as  the  angle  between  rays  and  Vn  does  not  appreciably  change.  Under  these 
conditions  the  beam  will  be  uniformly  deflected  by  the  disturbance  and  no 
changes  in  light  intensity  are  to  be  expected  on  the  photographic  plane. 

The  conclusion  is  then  drawn  that  shadowgraph  effects  are  connected  with 
variable  density  gradients  -  in  the  one -dimensional  case  with  non- zero 
values  of  d  n/bx  .  With  a  rounded- step  or  S-curve  density  function  such 
as  would  correspond  to  a  one-dimensional  shock  wave,  the  initial  region 
where  5  n/bx  >0  as  Sn/Sx  increases  will  correspond  to  increasingly  de¬ 
flected  rays  and  give  rise  to  lowered  illumination,  that  is  a  "shadow", 

2  2 

on  the  photographic  record;  while  /the  final  portion  where  c)  n/bx  <  0  with 
dn/dx  decreasing  Will  give  rise  to  decreasingly  deflected  rays .  Rays 
from  the  front  and  rear  portions  intersect  to  cause  a  double  mapping  or 
region  of  brightness  on  the  plate. 

By  this  reasoning  we  thus  expect  a  shock  wave  to  be  represented  by 
a  shadow  followed  by  a  bright  line  just  inside.  The  outer  boundary  of 
the  shock  where  normal  intensity  ends  is  an  accurate  projection  since  the 
bounding  rays  pass  through  the  uniform  medium  preceding  the  shock  and  are 
not  deflected.  Observe  that  we  are  still  neglecting  diffraction  effects. 

The  foregoing  discussion  does  not  exhaust  the  conditions  under  which 
shadows  and  bright  lines  may  be  obtained.  If  a  sharp  boundary  separates 
a  region  of  constant  index  n  from  one  of  constant  gradient  Vn,  the 
boundary  should  be  marked  by  a  shadow  separating  the  undeflected  linear 
rays  from  the  deflected  circular  arcs .  The  shock  shadow  followed  by  a 
bright  region  has  been  found  by  Sedney  and  Gerber  [lhjin  a  raytracing 
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Figure  4.  Shadowgraph  of  blunt  (35°  half-angle)  cone-cylinder  at  M  =  2.9  and 
p  =  5  atm.  The  horns  of  light  penetrating  the  shock  wave  represent  rays 
refracted  by  the  expansion  fans  at  the  corners  where  the  cone  and  cylinder  join 


investigation  of  refraction  by  a  cone  of  constant  index.  Here  the  conical 
flow  about  a  supersonic  conical  body  is  represented  by  an  equivalent  cone 
of  constant  refractive  index  chosen  to  match  the  jump  at  the  shook  wave 
but  to  ignore  the  slight  increase  in  density  as  the  projectile  cone  is 
approached.  The  cotte  is  assumed  to  he  illuminated  by  parallel  light  which 
travels  perpendicular  to  the  axis  of  the  cone  --a  cane  which  corresponds 
to  a  shadowgraph  source  at  infinity.  Rays  lying  in  planes  normal  to  the 
cane  axis  are  plotted.  Traces  in  the  image  space  of  these  planes  instead 
of  being  vertical  lines  are  hooked  curves,  like  inverted  Skew  parabolas, 
which  bend  into  the  geometrical  shadow  of  the  refracting  cone  leaving  the 
region  Just  behind  the  shock  in  shadow  and  doubly  mapping  a  flow  region 
further  inside  the  shadow.  The  envelope  of  these  curves  is  a  curve  of 
theoretically  infinite  intensity  which  corresponds  to  the  bright  inside 
boundary  of  the  shadowgraph  shock.  The  major  center  portion  of  each  trace 
is  obscured  by  the  solid  projectile  body  leaving  only  the  more  refracted 
rays  to  illuminate  the  flow.  For  cases  in  which  the  projectile-to-film 
distance  d^  is  large  compared  to  cone  length  and  refraction  displacements 
correspondingly  large  as  seen  in  Figure  4,  the  entire  flow  region  may 
appear  as  shadow  with  illumination  from  the  most  refracted  rays  appearing 
only  inside  the  body  shadow  and  tracing  out  two  nappes  of  a  luminous  cone. 

The  small  interior  shocks  caused  by  body  roughness  have  been  traced 
also  and  found  to  correspond  to  curves  lying  within  the  body  shadow.  In 
a  Shadowgraph  taken  at  5  norma]  densities  (Figure  4)  these  curves  are 
plainly  visible  as  hairy,  luminous  processes  of  approximately  the  shape 
expected  from  calculations . 

These  results,  while  Interesting  for  flow  interpretation,  do  not  hear 
significantly  on  problems  of  position  measurement.  As  already  noted,  the 
Shock  outline  is  accurate  if  diffraction  effects  can  he  ignored.  Thus 
the  forward  tip  of  the  usual  projectile  shadowgraph  may  he  measured 
with  some  confidence  providing  spark  duration  Is  small.  On  the  other  hand, 
base  position  would  appear  to  he  more  in  doubt;  for  a  qualitative  ray- 
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tracing  (no  quantitative  results  have  been  obtained)  just  behind  the  base 
of  a  cone-cylinder  projectile  indicates  that  the  grazing  rays  should  be 
displaced  by  a  Rwal  1  but  undetermined  amount.  Whether  forward  or  backward 
will  depend  oh  the  relative  strengths  of  the  shock  wave  and  the  expansion 
at  the  base,  and,  upon  the  angles  made  by  the  rays  with  the  density  gra¬ 
dients  at  these  places.  Since  the  two  regions  tend  to  compensate  one  might 
assume  the  effect  to  be  negligible.  Measurements  of  shadowgraph  projectile 
lengths  at  pressures  up  to  5  a±m  when  corrected  for  magnification  agree 
with  direct  length  measurements  to  within  .02",  or  about  2$. 

Diffraction  effects,  while  ordinarily  not  noticeable  in  shadowgraphs, 
are  nevertheless  present  and  should  be  considered  when  position  or  angle 
measurements  of  highest  precision  are  desired.  Diffraction  effects  are 
usually  not  noticeable  in  range  shadowgraphs  because  l)  spark  source 
duration  tends  to  blur  out  systems  of  fringes  not  parallel  to  the  line 
of  flight,  2)  source  size  tends  to  blur  any  system  of  fringes  through 
superposition  of  overlapping  systems,  and  3)  spectral  breadth  of  source 
tends  to  blur  all  but  the  zeroth  order  interference. 

When  rotating  mirror  techniques  are  used  to  minimize  blur  in  shadow¬ 
graphs  of  moving  projectile  taken  with  spark  sources  of  the  type  new  in 
use  on  CTPR,  Systems  of  fringes  both  inside  and  outside  the  Shock  wave  may 
be  seen  and  elsewhere  around  the  projectile.  The  persistence  of  a  large, 
bright  fringe  clearly  seen  in  Figure  k  around  the  body  outline  into  the 
turbulent  mixing  region  as  a  kind  of  encasement  may  conceivably  be  in¬ 
terpreted  as  evidence  of  a  laminar  superlayer  anteriorly,  between  the 
body  and  the  potential  flow;  and  posteriorly,  between  the  mixing  region 
and  the  outside  potential  flow. 

Because  the  ray  tracing  calculation  03  indicates  a  double  mapping 
inside  the  shock,  here  then  is  another  place  where  interference  may  occur 
if  the  phase  difference  between  coherent  rays  is  appropriate.  A  calcu¬ 
lation  of  fringe  positions  inside  the  shock  has  been  made  Ql^J  on  the 
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basis  of  the  geometrical  ray  tracing  alone  and  a  system  found  whose 
maxima  and  minima  lie  within  about  20  per  cent  of  the  measured  positions. 
One  would  not  expect  a  calculation  based  on  geometrical  optics  to 
represent  very  accurately  a  situation  where  diffraction  is  important; 
for  it  is  well  known  that  geometrical  optics  fails  to  predict  phenomena 
accurately  in  the  neighborhood  of  focal  points  or  caustic  curves  -  that 
is  to  say  near  ray  envelopes  such  as  we  encounter  behind  a  shock  wave. 

A  thoroughgoing  analysis  of  diffraction  by  shock  wave  boundaries 
from  the  point  of  view  of  wave  optics  has  not  been  made  to  the  author’ s 
knowledge  probably  because  such  a  calculation  would  he  very  difficult 
and  because  diffraction  phenomena  are  not  ordinarily  enough  in  evidence 
to  demand  attention.  Nevertheless,  when  high  precision  measurements  of 
unfocussed  shadowgraphs  are  required  for  the  detailed  examination  of  flow 
phenomena  mapped  by  large  refraction  effects,  account  of  diffraction 
effects  will  have  to  be  rendered  and  new  definitions  of  apparent  shock 
wave  and  projectile  boundaries  will  naturally  arise  from  these  studies . 

3-3  Measurement  of  Density  Fields 

The  interferometric  technique  of  measuring  density  throughout  an 
entire  oxisymmetric  field  Df  flow  is  relatively  well  understood  [>]and 
can  he  carried  out  in  representative  flows  produced  at  room  temperature 
to  an  accuracy  of  about  +  2$  save  near  the  axis  of  flow.  Various  studies 
of  the  effects  of  refraction  on  interferometric  accuracy  have  been  made 
[16,  2J  and  some  are  still  in  progress  [u]  •  These  indicate  that  largest 
errors  will  occur  near  shock  wave  boundaries.  Experimental  studies  with 
the  10”  Mach-Zehnder  instrument  of  the  CTPR  already  show  that  at  densities 
corresponding  to  one  normal  atmosphere  and  above,  appreciable  regions  of 
an  interferogram  near  the  shock  wave  will  be  lost  because  the  rays  mapping 
these  regions  are  refracted  off  the  rotating  mirror  face.  Thus  the 
rotating  mirror  acts  like  a  symmetrical  schlieren  stop  and  schlieren 
effects  on  high  density  interferograms  are  very  noticeable  indeed-  This 
particular  difficulty  is  easily  remedied  by  use  of  larger  mirror  faces; 
yet  it  is  clear  that  a  limitation  exists  here  for  mirrors  with  8  faces 
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2"  x  2”  in  size  are  rather  large  objects  and  to  achieve  rotational  speeds 
higher  than  the  6-8000  rpm  required  by  high  velocity  rounds  necessitates 
turbine  motors,  or  other  drives,  capable  of  producing  several  horsepower. 

I  do  not  wish  to  enter  the  complexities  of  the  light  source,  rotating- 
mirror,  camera  problem  for  interferometric  work,  but  rather  will  focus 
attention  briefly  on  another  problem  encountered  because  of  the  low  temper¬ 
atures  used  in  the  CTFR. 

As  is  well  known,  air  or  other  working  fluid  must  be  exceedingly 
stable  and  free  of  variable  density,  convection  currents  in  order  for 
interferometric  techniques  to  be  successful.  This  requirement  has  been 
one  of  the  most  difficult  to  satisfy  in  the  design  of  the  CTFR.  In  the 
original  plans,  the  interferometer  window  cells  of  the  inside  stainless 
steel  cylinder  were  joined  to  those  of  the  outside  mild  steel  chamber  by 
expansible,  stainless  steel  bellows.  On  these  were  attached  thermocouples, 
calrod  heating  elements  and  liquid  nitrogen  cooling  coils,  the  idea  being 
that  precise  metering  of  electrical  heat  and  of  liquid  coplapt  would  keep 
the  outer  optical  glass  windows  at  room  temperature  while  minimizing  heat 
flow  from  outer  to  inner  window  cells.  It  was  planned  to  reduce  radiated 
heat  losses  from  outer  to  inner  windows  by  coating  the  inner  facing  sur¬ 
faces  with  a  thin  layer  of  evaporated  gold. 

Subsequent  tests  have  shown  the  need  to  modify  the  original  designs 
in  almost  every  detail.  Heat  conduction  through  the  bellows  could  not  be 
controlled  either  rapidly  or  accurately.  They  have  been  removed  entirely. 
Gold  coatings  thick  enough  to  reflect  infra-red  radiation  were  found  to  be 
too  absorptive  of  all  frequencies  to  permit  their  use.  They  have  been 
replaced  by  two  pairs  of  solenoid-operated,  polished-aluminum  shutters  in 
the  vacuum  space  between  inner  and  outer  windows.  Annular  baffles  made  of 
cardboard  or  thin  plastic  materials  such  as  teflon  have  been  placed  between 
stations  inside  the  working  section  to  stabilize  the  working  fluid  and 
minimize  longitudinal  temperature  gradients.  Finally  a  cylindrical  en¬ 
closure  of  cardboard  or  plastic  sheet  perforated  with  12"  horizontal  holes 
to  permit  passage  of  the  projectile  has  been  used  between  upper  and  lower 
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Figure  5.  Interferogran.  of  l/2"  sphere  at  M  =  9-3  in  nitrogen  at  T  =  97°K  and  p  =  3.63"  Hg. 


interferometer  windows  inside  the  working  section. 

All  these  measures  together,  plus  improved  techniques  of  coolant 
circulation  have  resulted  in  sufficient  stability  of  the  cooled  gas  so 
that  reducible  interferograms,  of  which  Figure  5  is  an  example, can  be 
obtained.  These  are  not  of  as  high  quality  as  those  taken  at  room  temper¬ 
ature  because  slight,  local  density  fluctuations  are  still  present  and 
drift  through  the  interferometric  field;  nevertheless,  these  local 
variations  cause  only  fractional  fringe  shifts  which  are  small  compared 
to  the  shifts  caused  by  the  flow  itself.  Achievement  of  reasonably 
uniform  undisturbed  fields  at  the  extreme  low  temperatures  is  so  recent 
an  accomplishment  that  complete  flow  reductions  have  not  yet  been  made; 
nevertheless  interferograms  already  in  hand  show  all  the  essential 
features  of  those  obtained  at  room  temperature  and  should  yield  density 
fields  less  accurate  by  only  a  few  percent. 

Since  the  discovery  that  the  zone  method  of  interferogram  reduction 
deteriorates  in  accuracy  toward  the  axis  of  flow  [2]  ,  a  continuing  study 
[l8,  I9J  has  been  made  of  various  reduction  methods  which  involve  leant 
squares  procedures  of  fitting  consecutive  segments  of  the  fringe  shift 
data  with  polynomials.  The  advantages  of  least  squares  methods  are  twofold: 
1)  hand  plotting  and  smoothing  can  be  avoided  and  2)  random  errors  in 
density  values  on  the  axis  of  flow  appear  to  be  no  greater  than  elsewhere. 
The  principal  disadvantages  encountered  so  far  include  l)  spurious 
discontinuities  in  density  at  tie-points  of  adjacent  polynomials  and  2) 
propagation  of  errors  in  the  fit  through  the  tie  points  of  successively 
fitted  polynomial  segments  until  all  accuracy  is  lost.  These  difficulties 
arise  because  the  usual  fringe  shift  curve  is  not  readily  fitted  by  poly¬ 
nomials  of  degree  low  enough  to  be  convenient  i.e.  degree  three  or  four, 
and  so  must  be . approximated  over  short  sections .  Electronic  computing 
machine  programs  for  the  zone  method  (ENIAC)  and  least  squares  schemes 
(QRDVjftC)  have  been  devised  and  tested,  but  the  reduction  problem  is  not 
yet  solved  in  a  completely  satisfactory  manner ■  A  combination  of  the 
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two  methods ,  i.e.  approximating  the  fringe  shift  hy  the  zone  method  at 
intermediate  and  large  distances  and  hy  a  tied,  least-squares  cubic  near 
the  axis,  may  turn  out  to  he  the  optimum  choice  where  highest  accuracy 
is  desired. 
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IV  SOME  RECENT  FLOW  INVESTIGATIONS 


4.1  Kp  and  Wake  Angle  Measurements 

An  account  is  given,  in  a  companion  paper  ^20^]  of  this  symposium, 
of  application  of  the  CTPR  to  investigate  observable  properties  of  the 
flovs  about  proposed  nose  models  for  the  ICBM  program.  The  studies  made 
so  far  represent  attempts  to  simulate  normal  density  flovs  for  the  Mach 
number  span  4  <■  M  ^  10.  The  demand  of  this  program  for  data  at  higher  Mach 
numbers  has  greatly  stimulated  development  at  EBL  of  high-velocity  missile 
launchers.  This  laboratory  has  under  procurement  at  the  present  time  two 
guns  of  the  type  developed  by  Drs.  Slavsky  and  Seigel  of  NOL.  When  these 
have  been  obtained  and  proven,  extension  of  the  present  Mach  number  range 
by  a  factor  of  two  is  to  be  expected.  This  should  make  accessible  to  study 
flows  with  an  appreciable  degree  of  ionization.  Whether  the  interferometric 
techniques  discussed  here  can  be  used  to  explore  the  properties  of  a  self- 
luminous,  highly  excited  medium  remains  open  to  question. 

4.2  Kp  vs  M  for  a  Blunt  Cone-Cylinder 

In  any  discussion  of  an  inftrument  such  as  the  CTFR  the  question 
naturally  arises  as  to  how  well  flows  produced  by  models  at  low  temperatures 
and  high  Mach  numbers  simulate  those  produced  by  models  or  full  scale 
objects  at  the  same  Mach  numbers  arid  "ordinary"  temperatures.  To  provide 
a  partial  answer  to  this  complex  question,  studies  of  ^  vs  M  for  a  blunt 
cone-cylinder  i.e.  a  35°  half -angle,  cone-cylinder  have  been  made  at  300° 
105°  and  85°  Kelvin.  Before  discussing  some  of  these  results  a  few 
theoretical  remarks  will  be  interpolated  to  clarify  the  situation. 

If  the  equations  governing  the  conservation  of  mass,  momentum  and 
energy  for  a  viscous,  heat-conducting,  compressible  fluid  be  rendered 
dimensionless  through  the  introduction  of  suitable  transformations  on  the 
dependent  and  independent  variables,  then  a  set  of  partial  differential 
equations  emerges  which  is  not  only  free  of  dimensions  but  also  exhibits 
certain  dimensionless  functional  coefficients  which  appear  as  multipliers 
of  the  various  derivatives. 
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These  coefficients  are:  7  =  C ^/C^,  M  =  V/a,  Fr  =  V^/gL,  Re^  «  pYLfv- 
Reg  =  pVl/X  and  W  =  k 8/pV^.  Here  7,  M  and  Fr  represent  specific  heat 
ratio,  Mach  number  and  Froude  number  respectively;  while  p  and  X  stand  for 
the  first  and  second  coefficients  of  viscosity  and  Re^,  Reg  the  correspond¬ 
ing  Reynolds  numbers  -  The  quantity  W  in  which  heat  conductivity  k  and  a 
characteristic  temperature  difference  8  occur,  has  no  conventional  name  but 
is  discussed  by  Dryden  £30^  and.  is  recognizable  as  the  ratio  of  heat  flow 
to  energy  dissipated  by  viscous  friction.  Depending  upon  assumptions 
defining  relative  orders  of  magnitude  of  heat  flow,  viscous  dissipation, 
kinetic  energy,  other  dimensionless  quantities  such  as  Husselt  number  or 
Pr&ndtl  number  may  enter  through  the  boundary  conditions  of  a  given  problem. 

If  flows  about  two  similar  objects  of  different  size,  say  two  spheres 
in  an  infinite  medium,  are  to  have  the  same  dimensionless  solution  to  the 
flow  equations,  then  the  six  coefficients  Just  listed  must  be  presented 
without  change  in  the  two  situations.  In  addition,  initial  conditions 
and  conditions  at  the  boundaries  must  be  rendered  in  dimensionless  form 
and  dimensionless  coefficients  preserved. 

In  supersonic  flow  problems  encountered  in  firing  ranges  it  is 
customary  to  regard  gravitational  effects  as  negligible.  Thus  Froude 
number  may  be  disregarded.  If  heat  flow  were  small  compared  to  viscous 
dissipation  W  might  be  neglected.  On  the  other  hand,  it  has  long  been 
understood  that  conditions  on  the  heat  transfer  at  the  flow  boundaries 
determine  whether  W,  Nusselt  and/or  Prandtl  numbers  will  be  of  Importance. 

In  some  range  experiments  missiles  are  heated  or  cooled  in  order  to  control 
these  modelling  parameters.  In  experiments  to  date  with  the  CTFR  this 
has  not  been  done.  Comparisons  among  various  drag  data  are  made  simply 
Ignoring  the  fact  that  W,  Pr  and  Hu  have  not  been  controlled  in  the 
modelling.  As  for  second  coefficient  of  viscosity  the  view  is  taken  that 
either  the  Stokes'  relation  X  +  2/3 u  =  0  holds  or  the  effects  of  X  can 
be  ignored  or  If  not  ignored  are  strictly  proportional  to  those  of  p. 

Under  any  of  these  conditions  Reg  can  either  be  ignored  or  lumped  with 
Re^.  This  leaves  then  three  coefficients  to  preserve  in  the  preliminary 
modelling  experiments  viz.  7,  M,  and  Re^. 


The  firdt  of  these,  7,  is  usually  regarded  as  constant.  Tables  of 
this  ratio  for  air  or  nitrogen  [2lJ  show  deviations  from  7  =  1.40  of 
slightly  more  than  1$  in  the  temperature  range  100  ^  T  300°K.  As  the 
boiling  point  of  air  or  nitrogen  is  approached  the  deviations  become 
larger.  The  data  for  nitrogen  at  one  normal  density  i.e.  approximately 
l/j5  atm  suggest  that  at  T  =  85°K  the  1$  limit  may  not  be  exceeded.  For 
these  reasons  7  is  assumed  to  be  the  same  constant  in  both  the  room  and 
low  temperature  experiments. 

Mach  number  is  chosen  to  be  the  independent  variable  in  both  experiments. 
While  clearly  the  room  temperature  tests  cannot  reach  the  high  Mach  numbers 
accessible  at  low  temperatures,  a  considerable  overlap  in  the  range  of  values 
4  £  M  £  6  is  possible  where  values  ought  to  be  directly  comparable. 

Re^  =  pVL/p  is  the  one  remaining  coefficient  whose  value  must  he  preserved 
in  the  two  experiments.  Because  velocity  of  sound  variation  with  absolute 
temperature  is  the  essential  factor  in  attainment  of  the  high  Mach  number 
data,  the  dependence  of  Re^  on  M  is  shown  by  writing  Re^  =  pIMa/p  where 
V  =  sM. 

If  it  were  true  that  the  temperature  dependence  of  a  and  p  were  the 
same,  then  the  problem  would  simplify  to  that  of  maintaining  pL  =  const, 
for  a  given  Mach  number.  Simple  kinetic  theory  predicts  that  both  a  and 
p  should  be  proportional  to  T1/2.  In  the  present  study  of  the  blunt 
cone-cylinder  it  is  assumed  that  the  temperature  effects  of  a  and  p  da 
cancel.  For  constant  missile  size  p  is  merely  adjusted  to  have  the  same 
value  at  low  temperature  as  at  room  temperature.  This  means  working  at 
something  less  than  l/j  atm  when  T~85°K.  The  further  assumption  is  made 
that  in  the  non-overlapping  region  of  M,  i.e.  for  M  ^  6,  the  increased 
Re^  has  little  effect  on  or  other  flow  parameters  being  studied. 

The  assumption  that  the  temperature  effects  of  a  and  p  cancel  is 
not  strictly  valid.  For  velocity  of  sound  in  nitrogen  or  air,  the  relation 
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a  ^  CT"^  holds  ■within  about  3$  for  temperatures  100  — T^:  300°K  at  1  atm. 
pressure .  Yiscosity  of  air  or  nitrogen  varies  closely  as  T  raised  to  the 
8/9  power  for  the  same  range  of  temperature.  Thus  the  simplifying 
assumption  is  only  a  crude  representation  of  the  facts.  For  a  fixed 
projectile  size  L,  to  render  Re1  =  const,  precisely,  one  would  have  to 
adjust  density  so  that  p  =  const,  x  p/Ma.  The  measurements  of  vs  M 
made  for  comparison  purposes!  at  T  =  300°,  (105°)  and  85°K  for  the  blunt 
cone-cylinder  [lcfj  show  values  at  all  temperatures  falling  quite  well 
on  a  continuous  curve.  As  is  relatively  insensitive  to  Re^,  minor 
variations  in  this  quantity  appear  to  have  no  significance  even  up  to  a 
five-fold  variation  in  Re^  [l<5J.  This  is  in  contradiction  to  preliminary 
measurements  made  before  precise  determinations  of  ambient  density  could 
be  made  C2# 

The  relatively  small  failures  indicated  above  to  model  7  and  Re^ 
accurately  are  often  completely  masked  by  experimental  difficulties  in¬ 
dicative  of  the  crude  approach  followed  during  these  early  usee  of  the 
CTFR.  Until  recently,  temperature  gradients  in  the  range  introduced 
uncertainty  to  the  interpretation  of  thermocouple  measurements  and,  as  a 
consequence,  into  the  calculation  of  the  ambient  density  .  Some  of  the 
discrepancy  between  early  groups  of  measurements  can  be  traced  to  erroneous 
values  used  to  calculate  K^,  Improvements  in  temperature  control  and 
uniformity  have  brought  better  agreement  between  recent  series  of 
values . 

Among  other  sources  of  error  at  room  temperatures  are  the  unaccounted 
for  variations  in  7  and  mol.  wt.  of  the  working  fluid  caused  by  contaminants 
from  the  burnt  powder  gases.  In  spite  of  multiple  baffle  systems  some 
contamination  can  hardly  be  avoided  —  and  cannot  at  present  be  detected 
until  a  system  of  chemical  analysis  of  the  medium  or  precise  measurement 
of  its  characteristic  sound  speed  allows  continuous  monitoring  of  the 
composition  of  the  medium.  At  the  low  temperatures  contamination  from 
powder  gases  constitutes  a  two-fold  problem.  If  the  gases  are  condensible 
they  immediately  appear  either  on  walls  and  window  surfaces  as  "snow"  or 
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"frost"  of  various  kinds  or  as  a  fog  presumably  made  up  of  fine  crystals 
of  unknown  compos it on.  The  fog  does  not  settle  out  during  moderate 
lengths  of  time;  So  the  result  is  a  vorking  medium  with  poorly  known 
physical  constants  which  is  also  optically  of  poor  quality.  During  recent 
firings  at  low  temperatures,  the  baffle  systems  plus  improved  techniques 
of  flushing  the  gun  chamber  after  firing  have  minimized  problems  caused 
by  frost  or  fog;  nevertheless  after  10  rounds  have  been  fired  from  the 
cal  .75  launcher  seme  frost  can  be  observed  at  the  station  windows  and  fog 
though  slight,  is  easily  detectable. 

An  interesting  result  in  the  vs  M  curves  for  the  blunt  cone-cylinder 
has  emerged  as  improved  techniques  have  provided  smoother  data.  In  the 
region  near  M  =  4  a  slight  bump  appears  which  extends  over  about  half  a 
Mach  number.  This  phenomenon  was  first  noted  in  the  data  obtained  from 
firings  through  nitrogen  at  low  temperatures  £l0  J  .  More  recently  firings 
at  room  temperature  through  argon-  sulfurhemf luoride  mixtures  have 

revealed  the  same  hump  in  the  curve.  So  far  no  explanation  for  its 
presence  can  be  given. 

4.5  interferometric  Properties  of  the  N-Wave 

A  shadowgraph  Or  interferogram  of  an  axisyrametric  flow  presents  an 
enormous  amount  Of  data  in  visual  form.  Attempts  to  reduce  this  data  to 
quantitative  results  emphasize  the  desirability  of  principles  which  corre¬ 
late  whole  areas  of  flow  under  a  single  unifying  lav. 

Scafte  years  ago  the  discovery  was  made  [  2k]]  that  conical  flows  demon¬ 
strate  an  unusual  and  simplifying  symmetry  for  interferometric  purposes. 

The  gist  of  this  principle  is  that  fringe  Shift  _B,  divided  by  axial  distance 
z  is  a  unique  function  of  r/z  throughout  the  entire  conical  flow.  Here  r 
denotes  distance  in  the  radial  direction.  Plots  of  &/z  vs  r/z  yield  a 
single  curve  id  conical  flows  and  may  be  used  as  "test  gauges"  to  determine 
whether  a  flow  is  conical  and  where  it  departs  from  coni  city.  This  conical 
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flow  principle  has  been  used  successfully  as  a  tool  in  studying  other 
related  flows  ^25,  2^j  . 

Recently  a  similar  principle  has  been  found  to  hold  in  the  distant 
flow  between  the  front  and  rear  shocks  caused  by  an  axisymmetric,  super¬ 
sonic  projectile.  Along  a  trace  parallel  to  the  axis,  this  region  is 
characterized  by  a  vertical  pressure  rise  at  the  front  shock  followed  by 
a  linear  decrease  to  below  ambient  levels  and  finally  another  vertical 
pressure  rise  at  the  rear  shock.  Because  of  the  N-shaped  pressure  curve 
the  region  has  been  called  "the  N-wave" . 

Theoretical  studies  by  G.  B.  Whitham  [27] have  given  approximate 
expressions  for  the  pressure  in  the  N-wave  region.  Working  to  the  same 
approximation  D.  H.  Steininger  £28]  and  the  author  [29J  have  derived 
expressions  for  the  interferometric  fringe  shift  in  the  N-wave  and  compared 
these  to  data  from  cone -cylinders  and  spheres.  When  suitable  approximations 
are  applied,  a  simple  mapping  function  is  obtained  which  characterizes  N-wave 
flow.  If  £  represents  the  fractional  distance  between  front  and  rear  shocks 
along  a  trace  at  constant  radial  distance  r,  then  it  can  be  shown  that 

=  |1^2(1  -  45/3) j  thus  by  plotting  Sr^^G'  versus  5,  experimental- 
data  can  be  compared  with  the  function  £1//2(l  -  k|/3) .  G1  is  a  function 
which  depends  upon  Mach  number,  body  shape,  the  wavelength  of  light  used 
and  physical  constants  of  the  medium.  Convergence  to  N-wave  flow  has  been 
examined  for  a  cone-cylinder  of  7-5°  half -angle  and  for  small  spheres.  At 
70  diameters  from  the  axis  of  the  sphere  the  agreement  between  plotted  fringe 
shifts  and  the  mapping  function  is  excellent  and  one  concludes  that  N-wave 
flow  is  obtained. 

Not  only  may  fringe  shifts  within  a  given  region  of  flow  be  compared 
by  means  of  this  principle,  but  also  fringe  shifts  from  flows  about  similar 
bodies  at  different  Mach  numbers  and  under  different  conditions  of  ambient 
density.  This  latter  is  a  direct  consequence  of  the  fact  that  the  function 
G*  accounts  for  differences  in  Mach  number,  density,  etc.  It  is  perhaps 
apparent  that  interesting  studies  of  N-wave  flow  as  a  function  of  Reynolds 
number  might  now  be  undertaken. 
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The  N-vave  mapping  lav  yields  a  principle  vhich  permits  study  and 
comparison  of  axisymmetrlc  flow  data  over  a  large  region  of  the  flew.  It 
holds  this  property  in  common  with  the  conical  flew  principle.  Reasoning 
by  analogy  one  may  expect  that  the  N-vave  property  v±H  eventually  beccme 
useful  as  a  diagnostic  criterion  and  standard  of  comparison  in  the  same 
way  that  the  conical  flow  test  already  has  done. 
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ABSTRACT 

A  simple  approximate  theory  for  the  free  flight  motion 
of  ballistic  missiles  having  mass  asymmetry  is  given  and 
proofed  by  experimental  firings  in  the  NOL  Pressurized  Aero- 
ballistic  Range. 


INTRODUCTION 


A  simple  theory  for  the  free  flight  motion  of  "ballistic  missiles 

having  mass  asymmetry  appears  to  he  unavailable  in  spite  of  classical 

1/2-  3 

treatments  in  Dynamics  and  the  excellent  vork  of  Rankin.  In  Ref.  4  a 
theory  for  pure  angular  mass  asymmetries  vas  given,  and  in  Ref.  5  an 
experimental  program  to  proof  this  theory  vas  reported.  Difficulties 
arising  during  this  experimental  program  stimulated  Ref.  6,  where  the 
theory  is  extended  to  include  linear  mass  asymmetry  (C.g.  laterally  dis¬ 
placed  from  missile  axis  of  configurational  symmetry),  and  where  the 
importance  and  contributions  of  various  dynamical  and  aerodynamical  terms 
are  evaluated. 

Id  this  note  the  treatments  of  Refs.  4,  5,  and  6  are  combined  to 
yield  a  simple  tested  theory  for  the  free  flight  motion  of  ballistic 
missiles  having  slight  mass  asymmetry. 
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THEORY 


The  general  equations  for  the  free  flight  motion  of  a  missile 
having  mass  asymmetry  may  he  written  as 

m  -  vr+vq  -  x(q2 +r2)  *y(pq  -  f )  +  z(pr  +  4)]  -  XB 

m  [  v  -  vp  +ur  -  y(r%  p2)+  z(qr  -  p)+  x(qp  +  r)]  -  Yft 

2  2  7 

m[$  -  uq  +  vp  -  z(p  +  q  )  +  x(rp  -  q)+  y(rq  +  p)J  -  &B 

(') 

2  2 

n»p  t(l7-  IY)qr  -  IjmCf  +pq)+  Iy^r  -  q  )  +  IvyCpr  -  4) 
m  l  y(v  *  uq  +  vp)  -  z(v  -  vp  +  ur)  ]  : 

IYq  +(lx-  It )rp  -  Ixy(p+  qr)4  I^Cp2-  r2)  +  I^qp  -  r) 

m  l  z(u  -  vr  +  vq)  -  x(£  -  uq  +  vp)  ]  -  MB 

2  2 

l^r  +(IY-  Ix)pq  -  I^q  +  rp)  +Ixy(q  -  P  )  *Ixz(r<l  -  P> 
m  [x(v  -  vp  +ur)  -  y(u  -  vr  +  wq)  ]  ■  - 

where  Xt  y,  Z  are  the  coordinates  of  the  center  of  gravity  location  in 
tbe(XYZ)^yBtem  .* 

By  the  usual  linear  assumptions  of  Bmall  normal  velocities,  linear 
and  angular,  and  hy  assuming  constant  axial  velocity  and  spin^these 
equations  reduce  to  ** 

o 

m(V  -  vp  +  ur  -  yp  )  = 

m(*  -  uq+vp  -  zp2)  =  Z 

2  (3) 

Iq  -  (I  -  Ix)pr+  p  z 

^  (I  -  Ix)p<l  “  P^-^XY  r  N0 

Written  in  "Body  Axes"  coordinate  system  which  is  fixed  in  the  missile  (Fig.l). 
**We  have  also  assumed  that  the  origin  and  orientation  of  the  XYZ  Body  Axes 
are  so  selected  that  X  -  0  IY  =  1^ 
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For  the  Aeroballistic  Axes,  which  are  tied  to  the  missile  hut 

4 

do  not  roll  with  it,  these  equations  may  he  written  as* 
m(v  +  ru)  -  mp2  cos  0m  =  Y 
m(w  -  qu)  -  mp  Smq  sin  0m 
14  +  pr  Ix+  P2^o  cos  0m  (Ix-  I) 

p 

If  -  pq  lx+  p  S  n]o  sin  0m  (lx-  I) 

where 

|  SnJ  =  I  XXY~  1  XXaI 
I  -  IX 

<^>m-^o  +  tan-1  /  XXY  ,] 

'  _IXZ  / 


The  physical  significance  of  these  equations  may  he  noted  from  Fig.  1, 
where  Sj^  is  seen  to  he  the  lateral  displacement  of  the  center  of  gravity 
of  the  missile  and  where  &  no  is  Been  to  he  the  angle  between  the  X  -  axes 
and  the  X'  -  axis  which  is  one  of  the  three  Principal  Axes  located  at  the 
center  of  gravity  of  the  missile,  (it  may  he  shown  that  if  the  equations 


for  the  motion  of  the  missile  are  first  written  in  this  Principal  Axis 
System  and  then  translated  by  and  rotated  by  both  <5^  and  §  , 

the  aeroballistic  equations  above  would  be  obtained.**  ) 


*  We  have  assumed  that  31^  i  for  all  roll  orientations  of  the  Aero¬ 
ballistic  Axes  and  the  Boay  Axes. 


**  It  is  assumed  that  the  mass  asymmetry  is  slight  and  thus  that 
c$mo  are  small.  It  follows  that  Iyiss  Iz  and  Iy^  ^  0. 


and 
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The  aerodynamic  force  and  moment  system  is  given  "by 


Y  n  YyV  f  Yrr  +  YyV  +  Yfr  +  Ypypw  +  Ypqpq  +  Yp^-pv  ■+  Ypqpq 
t  Y^  coafi  -  Z^  £e  sin  <£ 


•  • 


Z  e  Z \fw  t  Zqq  +  Z$w  +  Zqq  4  ^pvpv  *  Zprpr  ■»  Zpypv  *■  Zpf-pr 
§E  cos  0 


+  Z 


se  °E 


+  Yg  sin  £ 


(S) 


M  =  MyW  +  M-q  +  M^w  t  M*q  +  M™pv  +  M_rpr  +  M^pv  +  M^pr 

-tM^COS  §  -NSj?<SBBiD0 

R  r  Nvv  +  Hrr  +  Hy^  ♦  4  IW)w  +  Wpa1^  4  MpwP*  +  ^6^4 

-t-Mj^cos  <$  +  ^  sin  § 

Following  the  development  of  Ref.  8  Eqs .  (3)  and  (S)  may  he 
combined  to  yield* 

ur  +  nx  dr  +  n2  *  rc3eipt 

K1  551  ‘  ^5a_  -  (Mg  4  «  M&)  _  ip  /  Z^_  +  Iy  \ 

m  I  \tb  I  ) 


fa) 


H2  *=s  -  u  M 


N3  sk  jpe  * 
mV 


-iL-Wc  +  V  (zJLMJv--  -  .1  u  Mpy-  p^I  ZpV  \ 
I  '  Im  I  Im  / 

a  - -£)  _£ 


where 


and  where 


e*  ^  e  ^  +  p^m  SM  uf  =  Of^iuf 

E  S£  -  P2  m  (lx-  I) 

e  ;  Z , 

S  =  -i  Mr 

H  <J4 

Equalities  of  the  Rotary  Derivatives  "are  obtained  by  the  assumption 
that  the  missile  configuration  has  basic  trigonal  or  greater  rotational 
symmetry  and  mirror  symmetry. 


(7) 


(*) 
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The  general  solution  for  the  complex  angle  of  attack  is  found,  to 


be  Tricyclic  Motion* 


4>t? 


where 


=  77  =  ^  ^ 


kt,£ 


not 


CAT  r'pw  ^ 

XK 


*>>,%.  *  (' 1 1 ) 

K  _  <^o  -  (\i2j^X}X  tXX_ ~_±p) 
(7,*-  ~  (a>,L  ~ 


**» 


Utr-  *i )  -  \Ja<p  -a)-> 7] 


5-  liA 

4 

r 


r=r 


*  The  general  solution  for  the  transverse  displacement  is  also  of  the 
Tricyclic  form  and  is  given  by 

S  fvC^  ■=  <X  ^  +  kf 


190 


Tbe  characteristics  of  Tricyclic  type  motion  are  discussed  in 
Ref.  8. 

These  equations  are  alleged  to  enable  the  computation  of  the  detailed 
free  flight  motion  of  a  missile  having  slight  mass  asymmetry,  or,  con¬ 
versely,  to  enable  the  reduction  of  observed  free  flight  motions  for  tbe 
determination  of  the  aerodynamic  coefficients.  In  order  to  proof  their 
accuracy  an  experimental  program  was  undertaken  in  the  Pressurized  Aero- 
ballistic  Range  at  the  Naval  Ordnance  Laboratory. 
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EXPERIMENT 


The  configuration  and  the  three  model  designs  used  in  the  exper¬ 
iment  program  are  illustrated  in  Fig.  2.  The  configuration  is  a  simple 
cone-cylinder  of  a  fineness  ratio  of  five.  The  three  model  designs, 

Basic,  One  Hole,  and  Two  Hole,  represent  three  values  for  Smo(.0,.30°, 
and  ) ,  and  andicJnoJO .  Table  I  gives  additional 

physical  data  on  the  models  tested  in  the  program.  (Basic,  Rounds,  1912- 
1915;  One  Hole,  Rds.  1372-1375;  Two  Hole,  Rds.  1586-1592) . 

The  models  were  gun  launched  in  the  NOL  Pressurized  Aeroballistic 
Range  (Photo  l),  and  their  free  flight  motion  was  observed  (Photo  2).  A 
typical  history  of  the  observed  complex  angle  of  attack  is  illustrated  in 

Fig*  3- 

The  complex  angle  of  attack  and  the  lateral  displacement  were  reduced, 
using  both  an  Epicyclic  Reduction  and  a  Tricyclic  Reduction.  The  values  of 
the  "fit"  (i.e.  P.E.--  and  P.E.-—  )  are  given  in  Table  II.  It  should  be 

oC  o 

noted  that  the  use  of  the  Tricyclic  reduction  improves  the  complex  angle 
of  attack  "fit"  by  30*7$  (One  Hole)  and  42. 6jo  (Two  Hole). 

The  aerodynamic  coeff icients*obtained  from  the  reductions  are  given 
with  their  Probable  Errors  in  Table  III. 

Values  for  the  Third  Arm,  which  exists  only  because  of  the  Mass 
Asymmetry,  were  obtained  from  the  experimental  tests  and  were  also  com¬ 
puted  from  the  measured  physical  parameters  using  the  Tricyclic  Theory  of 
this  note.  These  theoretical  and  experimental  values  for  are  given  in 
Table  IV. 

*  The  variation  in  the  values  of  the  coefficients  from  the  three  model 

designs  should  be  noted,  but  no  analysis  is  given  herein. 
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It  should  be  noted  that  although  the  size  of  the  Third  Aims  are 
of  the  same  order  of  magnitude  as  the  values  for  the  Tricyclic  "fit" 
(Table  II ),  the  probable  errors  for  the  Third  Arms  are  small  and  thus 
indicate  an  accurate  determination.  The  consistency  of  the  determin¬ 
ation  of  K3  round  to  round  appear  good  except  for  Rd.  1375  and  Rd.  1507 • 
In  Table  V  the  agreement  between  theory  and  experiment  is  given. 
The  values  in  parentheses  are  obtained  by  omitting  Rds.  1375  and  1587. 


**  It  should  be  noted  that  the  Probable  Errors  of  the  aero¬ 
dynamic  coefficients  are  improved  by  the  order  of  $0%  by 
using  the  Tricyclic  Theory, 
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CONCLUDING  REMARKS 


The  data  from  the  experimental  program  indicates  that  the  Tricyclic 
Theory  for  mass  asymmetry  accurately  represents  the  free  flight  motion 
of  the  models  tested  and  that  the  various  static  and  dynamic  aerodynamic 
coefficients  are  veil  determined  from  the  motion  using  this  theory. 

The  theoretically  predicted  values  for  the  Third  Ana  arising  from 
mass  asymmetry  are  in  good  agreement  with  the  values  obtained  from  the 
experimental  firings.  For  the  One  Hole  models^ this  difference  of  ex¬ 
periment  from  theory  is  11.1^ and  for  the  Tvo  Hole  models^  the  difference 
is  I.55&  (Table  V). 

An  alternate  evaluation  of  the  Tricyclic  Theory  for  mass  asyimnetry 
might  be  obtained  by  comparison  of  this  approximate  theory  with  numer¬ 
ical  integrations  of  the  exact  equations,  Eq.  (l),  on  an  advanced  com¬ 
puting  machine  using  representative  values  of  Sjj  and  8^* 
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X'Y'Z'  PRINCIPLE  AXES 
X  Y  2  AEROBALLISTIC  AXES 


I  Bm  sin  <£ 
0  cos  <f> 

SM  -  sin  <f> 


-Bm  cos 
sin  <f> 

cos  <f> 
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BASIC  MODEL 


TWO  HOLE  MODEL 


DESIGN  OF  MISSILES  USED  IN  EXPERIMENTAL  TESTS 
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MODEL  WITH  2  HOLES 
*/4  CAL.  IN  DIA.  I  CAL.  LONG 
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P.E.  of  yaw  and  swerve  reductions,  angular  velocities 
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Table  III  (cont»d) 
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1.6 

3.16 

3.15 

3.13 

3.09 

. 

H 

°v 

P.B. 

% 

P.E. 

% 

3372^ 

1373 

1374 

1375 

-3  .TO' 

-3.02 

-3.11 

-3.16 

2.2“ 

2.3 

1.7 

2.9 

-3700 

-2.88 

-3.00 

-3.05 

1 

CM 

q  J 

SUE. 

H 

1372 

1373 

1374 

1375 

-2d.  7 

-30.3 

-30.3 

-32.4 

2T.7 

25.3 
16.6 

29.4 

-25.6 

-26.7 

-26.9 

-29.1 

7.1 

7.8 

11.9 

16.9 

Cjy 

^poC 

P.E. 

$ 

°v 

P.E. 

A 

1372 

1373 

1374 
,  1375 

1.15 

1.45 

1.35 

1.58 

47.8 

42.2 

43.0 

45.4 

0.64 

0.80 

O.78 

0.98 

18.0 

15.6 

23.6 
29.9 

Aerodynartiic  coefficients 


Table  III  (cont'd) 

T(A)o  /-foes 


Table  IV 


round 

ft 

K3 

rad. 

- - - 

K3 

rad. 

m 

theoretical 

experimental 

■  ■ 

basic  model 

1912 

■BRM/kHi 

.00055 

.00052 

.00043 

1914 

.00068 

.00051 

1915 

.00022 

.00033 

one  hole  model 

1372 

.00599 

.00034 ' 

1373 

.00597 

.00036 

1374 

.00594 

.00501 

.00092 

1375 

.00596 

.00085 

two  hole  model 

1536 

. 00733 

.00686 

.00096 

1587 

.00782 

.00970 

.00102 

1589 

.00779 

.00758 

.00072 

1590 

.00776 

.00738 

.00120 

1591 

.00769 

.00768 

.00082 

1592 

.00788 

.00684 

.00144 

Third  Ann, 


theoretical  and  experimental 
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T  fteZ-E  TT 


MODEL 

FIT 

RDS. 

K3 

COMPARISON 

EXR 

TH. 

BASIC 

.  107° 

4 

.0263° 

0° 

ONE  HOLE 

.  206° 

4 

.303° 

.341° 

ii.i  % 

(11.6%) 

(14.9) 

TWO  HOLE 

.263° 

6 

.439° 

.446° 

1.5%. 

( 13.6%) 

(6.7) 

COMPARISON  OF  THEORY  8  EXPERIMENT 
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Part  Two 


DYNAMIC  STABILITY* 

John  D.  Nicolaides 
Bureau  of  Ordnance 


ABSTRACT 

The  Epicyclic  Theory  for  the  flight  dynamics  of  ballistic 
missiles  has  yielded  various  ^Dynamic  Stability  Criteria*  which 
are  often  used  to  evaluate  missile  performance.  Recent  mis¬ 
leading  uses  of  the  theory  and  criteria  in  appraising  missile 
performance  require  a  simple  restatement  of  the  theory,  its 
assumptions,  and  its  use. 

The  parameters  of  Nutation  Half-Life,  Precession  Half- 
Life,  and  Total  Motion  Half-Life  are  suggested  as  better  cri¬ 
teria  for  missile  dynamics  than  those  classically  based  on  the 
Linear  Theory. 

Also  a  summary  of  Important  Non-Linear  Cases  of  ballistic 
missile  flight  performance  is  given. 


•K-The  study  in  this  Technical  Note  was  suggested  by  A.  Wertheimer. 


LINEAR  AEROBALLI S TIC  THEORY 


The  complex  angle  of  attack  is  given  by: 

OC  -  +  K'p  & 

and  represents  two  rotating  arms,  the  nutation  and  the  pre- 
cessidn,  which  can  increase  or  decrease  in  size.  The  arms 
rotate  at  the  rates  <0^  and  cop,  and  damp  or  expand  at  the 
rates  and  Xp  • 

This  solution  for  the  complex  angle  of  attack  was  ob¬ 
tained  by  making  the  following  assumptions: 

1*  constant  missile  velocity 

2.  constant  missile  rolling  velocity 

3.  constant  missile  mass  and  moments  of  inertia 
(symmetrical  transverse  mass  distribution  and 
symmetrical  axis  mass  distribution) 

J4..  no  gravity  attraction 

5»  linear  aerodynamic  force  and  moment  system 

6.  trigonal  or  greater  missile  rotational  symmetry 

7.  mirror  symmetry  of  missile 

8.  small  angular  displacements  and  velocities,  except 
roll 

These  are  the  assumptions  of  the  Basic  Linear  Theory. 
The  Theory  can  be  extended  by  relaxing  some  of  these  assump¬ 
tions,  however,  for  our  purposes  here  we  need  only  consider 
the  Basic  Theory. 

The  requirement  for  “Dynamic  Stability”  is  that  both 
the  Nutation  and  the  Precession  Arms  damp  or*  ~  ' 

X^«p  <  0  (2) 


This  condition  is  the  basic  definition  of. Dynamic 
Stability  as  generally  used  and  it  insures  that  eventually 
all  initial  motion  will  fade  away.  The  dynamic  stability 
criteria  of  Refs.  3  and  5  are  other  more  complicated  ways 
of  stating  this  same  simple  condition.  These  criteria  and 
the  simple  expression  above  all  fail  to  indicate  ”how  fast” 
the  motion  dies  out.  In  order  to  overcome  this  shortcoming, 
the  concepts  of  Nutation  Half-Life  and  Precession  Half-Life 
are  introduced  as 


*See  Ref.  6  for  definition  of  all  symbols  used  and  a  deriva- 
tion  of  Eqs.  (1),  (2),... 
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Nutation  Half-Life  = 


(calibers)  (3) 


Precession  Half-Life  = 


(calibers)  (I4.) 


These  relations  reveal  how  far  the  missile  flies  in  missile 
diameters  until  the  initial  Nutation  size  is  reduced  by  one- 
half  and  how  far  the  missile  flies  until  the  Precession  is 
reduced  by  one-half.  Except  for  a  non-rolling  fin  stabilized 
missile,  these  distances  are  never  the  same,  and  as  a  result 
the  reduction  of  the  total  initial  motion  to  one-half  size 
must  lie  somewhere  in  between. 

The  distance  between  these  two  distances  that  Total 
Motion  is  reduced  to  one-half  its  original  size  will  depend 
on  the  relative  size  of  the  Nutation  and  Precession  arms 
initially  and  their  damping  rates.  This  critical  distance 
that  the  missiles  travels  in  order  for  the  Total  Motion  to 
damp  one-half  its  Initial  value  may  be  calculated  from  Eq.  (5) 
using  Eq.  (6).  (See  Fig.  1.) 


TOTAL  MOTION  HALF-LIFE 


V* 


K*. 


(calibers)  (5) 


The  time  of  one-half  life,  ,  Is  obtained  from 

^2. 


1 

2 


where 


-  K|e 


%  I4le 


Xoi 


P  >  Vz. 


(6) 


<ye  |  =  I  ^nI  + 1  ^ p 


(7) 
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Some  discussion  of  this  Total  Motion  is  in  order.  The 
damping  rates  ,Xp  are  independent  of  the  initial  con¬ 


ditions  of  motion  and  are  determined  from  the  aerodynamic 
coefficients  and  the  physical  parameters  of  the  missile. 


*  2 1 V  Cyr,|+ UC 1 1  fc ' ^  \  C*w2^( < 8 > 


The  values  of  and  K?  depend  on  the  initial  conditions 
of  motion  given  by 


Thus  it  is  clear  that  the  distance  of  Total  Motion 
Half-Life  will  depend  quite  critically  on  the  initial  com¬ 
plex  angle  of  attack,  ,  and  the  initial  complex  rate 

of  change  of  angle  of  attack,  . 

For  the  "Dynamically  Stable"  missile,  ,  the  cal¬ 

culation  and  interpretation  of  the  Total  Motion  Half-Life 
may  be  carried  out  in  a  routine  manner  using  Eq.  (5)* 
However,  for  the  case  where  rt  \p  )  O  the  picture  is 

not  so  simple  and  the  indiscriminate  use  of  the  usual  term 
"Dynamically  Unstable"  may  be  downright  misleading. 

To  point  up  this  often  misunderstood  and  misapplied 
flight  condition,  it  is  instructive  to  consider  the  special 
case  where  KK  is  small,  Kp  is  large,  >  o  ,  and  \p<?0 

Here  we  have  a  case  where  >  O  ;  but  is  the  missile 

“ actually  dynamically  unstable?  Consider  the  pictorial  rep¬ 
resentation  in  Fig.  1  (Dashed  lines). 

Here  we  note  that  the  Total  Motion  may  ££  first  damp  down 
in  an  excellent  manner  yielding  a  good  valve  for  Half-Life; 


211 


ultimately,  however,  the  Total  Motion  will  undamp.  If  the 
missile  has  reached  its  destination  before  the  Total  Motion 
begins  to  undamp  to  a  significant  degree,  then  "all  is  well" 
in  spite  of  the  fact  that  the  classical  criterion  for  "dynamic 
stability"  is  not  satisfied. 

Although  we  have  considered  a  special  case,  the  argument 
and  basic  idea  is  sound  in  general  and  thus  should  be  con* 
tained  in  any  evaluation  criteria  for  missile  performance. 

This  concept  can  be  easily  included  by  calculating  the  two 
values  for  Total  Motion  Half-Life  which  such  a  missile  may 
Eave”7rom  Eq.  (5).  or  by  computing  the  Total  Motion  Double- 
Life  from  Eq.  (11). 


TOTAL  MOTION  DOUBLE-LIFE  = 


where  the  time  of  double-life  , 

2/°4/  =  }^N{  & 

It  should  be  clear  from  the  above  discussion  that  the 
indiscriminate  use  of  the  terms  "Dynamic  Instability"  with¬ 
out  consideration  of  the  missile  performance  requirements 
should  be  carefully  avoided.  Constructive  evaluation  of  the 
performance  of  "Dynamically  Unstable"  missiles  may  be  pre¬ 
sented  by  using  the  parameters  of  Total  Motion  Half-Life  and 
Double -Life. 


Viz 

a 

A 


i 


j  ]/p  I  & 


(caliber)  (11) 

is  obtained  from 

Ap  jf'-L 


(12) 
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SOME  NON-LINEAR  CASES 


The  preceding  discussion  is  based  entirely  on  the  Linear 
Theory  and  its  assumptions.  There  are,  however,  a  number  of 
significant  flight  cases  where  the  assumptions  of  the  Linear 
Theory  fail  and  thus  the  theory  cannot  be  used  or  must  be 
used  with  great  care.  Since  these  cases  may  not  be  generally 
known  and  since  they  seem  to  be  occurring  with  greater  fre¬ 
quency,  it /seems  reasonable  to  enumerate  them  below. 

Case  1.  -  Changing  Velocity 

Ballistic  missiles  in  flight  may  not  have  constant  ve¬ 
locity.  When  the  dynamic  stability  criterion  is  evaluated 
as  a  function  of  Mach  number,  it  is  not  unusual  to  find  that 
many  missiles  have  \N  <y\  \p  >  O  over  some  range  of 

Mach  numbers.  For  example,  many  fin  stabilized  missiles 
and  bodies  of  revolution  are  unstable  at  transonic  speeds* 
Also  some  fin  stabili2ed  missiles  are  unstable  at  very  high 
Mach  numbers.  Going  strictly  by  the  criterion,  these  mis¬ 
siles  may  be  adjudged  to  be  dynamically  unstable  and  thus 
unacceptable.  However,  the  time  during  which  the  missile 
flies  in  this  critical  region  may  be  small,  and  thus  the 
increase  in  yaw  during  that  increment  of  time  may  be  small, 
and  may  have  negligible  effect  on  the  retardation  and  ac¬ 
curacy  of  the  missile.  In  other  words,  local  flight  in¬ 
stabilities  may  mot  mean  a  bad  missile.  What  is  important 
is  the  integrated  effect  of  missile  dynamics  over  the  en¬ 
tire  flight  of  the  missile.  Dr.  Kent,  in  Reference  7,  has 
suggested  such  an  approach.  The  seriousness  of  the  local 
Mach  number  instability  will  clearly  depend  upon  the  mag¬ 
nitude  of  the  instability,  X^  and  Xp  t  and  the  size  of  the 
nutation  and  precession  arms,  Kjj  and  Kp  at  the  point  whore 

the  instability  begins.  Qualitative  estimates  of  the  ef¬ 
fect  can  be  obtained  by  using  the  linear  aeroballistlc 
theory* 

Case  2.  -  Non-Linear  Moments 


Non-linearities  in  the  aerodynamic  coefficients  can 
cause  serious  failures  in  prediction  by  the  Linear  Aero- 
ballistic  Theory.  It  has  occurred  in  a  number  of  cases 
in  the  recent  past  that  missiles  which  were  predicted  to 
be  dynamically  stable  were  later  found  to  be  unstable  due 
to  a  non-linear  increase  in  the  magnus  moment  due  to  angle 
of  attack.  This  instability  was  such  that  at  small  angles 
of  attack  the  missiles  were  stable;  however,  at  large  angles 
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of  attack  the  missiles  went  unstable*  The  achievement  of 
stable  or  unstable  flight  simply  depended  on  the  initial 
launching  conditions.  For  large  cross-wind  launchings  and/or 
large  muzzle  blasts  and  bore  clearance,  unstable  flights 
were  obtained;  whereas  for  small  cross-wind  effects  and 
small  initial  launching  conditions,  good  flights  were  ob¬ 
tained*  It  is  clearly  necessary,  then,  that  the  aerodynamic 
stability  derivatives  be  obtained  over  the  entire  range  of 
angle  of  attack  in  which  the  missile  may  be  flying  under 
realistic  launching  conditions.  Reference  8  suggests  that 
the  linear  aeroballistic  theory  might  be  modified  to  employ 
these  non-linear  data  and  thus  yield  a  criterion  for  perform¬ 
ance  of  a  missile  being  acted  upon  by  non-linear  aerodynamic 
forces  and  moments  systems*  In  the  cases  tested  the  Qjiasi- 
Linear  Aeroballistic  Theory  has  yielded  results  consistent 
with  exact  numerical  integrations  and  actual  performance. 

Case  3«  -  Non-Linear  Moments 

Non-linearities  in  the  Aerodynamic  coefficients  and 
particularly  in  the  Magnus  moment  may  go  in  the  opposite 
direction  as  in  the  previous  case.  In  other  words,  the  lin¬ 
ear  criterion  has  predicted  Xw  XP  >  O  for  many  bodies 
of  revolution,  bullets,  shells,  etc.,  at  small  angles  of  at¬ 
tack;  however,  these  missiles  have  become  dynamically  stable 
at  values  of  the  angle  of  attack*  The  non-linearities, 

here,  were  such  as  to  make  the  missiles  dynamically  stable 
at  angles  of  attack  other  than  zero,  yet  dynamically  un¬ 
stable  at  angles  of  attack  near  gero.  {The  air  drag  re¬ 
sulting  from  this  small  yawing  motion  may  not  be  important 
compared  to  the  other  drag  sources  on  the  projectile.) 

Case  t;.*  -  Non-Linearities  and  Asymmetries 

As  Indicated  in  Case  1,  the  seriousness  of  the  dynamic 
stability  depends  upon  the  integrated  effect  over  the  entire 
flight  of  a  missile,  and  thus  upon  the  magnitude  of  the  in¬ 
stability  and  the  relative  size  of  the  nutation  and  preces¬ 
sion  arms.  The  random  introduction  of  conf igurations  asym¬ 
metries  can  lead  to  larger  nutation  and  precession  arms  and 
thus  larger  complex  angles  of  attack  than  normally  encoun¬ 
tered  due  to  Initial  launching  conditions,  and  thus  the 
flight  instabilities  of  Case  2  and  Case  3  can  be  seriously 
affected  and  worsened  by  any  random  configurational  asym¬ 
metries  which  might  exist.  One  obvious  example  is  the  Case 
where  the  rolling  velocity  of  a  missile  and  the  nutation 
rate  are  equal  and  thus  we  have  maximum  amplification  of 
the  "TRIM*  due  to  any  configurational  asymmetries  which  ex¬ 
ist.  This  leads  to  large  angles  of  attack  for  which  insta¬ 
bility  of  the  type  mentioned  in  Case  2  would  lead  to  trouble, 
where  otherwise  the  flight  of  the  missile  would  be  completely 
satisfactory. 
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Case  5 »  -  Roll  Lock-In 


In  addition  to  the  non-linearities  in  the  forces  and 
moments  considered  in  the  linear  theory,  there  may  be  forces 
and  moments  entirely  left  out  of  the  theory  which  must  be 
considered.  One  of  these  neglected  moments  is  the  roll 
moment  due  to  complex  angle  of  attack  and  missile  roll  ori¬ 
entation.  This  roll  moment  is  such  as  to  prohibit  the  canted 
fin  missile  from  picking  up  its  intended  rolling  velocity. 
This  induced  roll  moment  causes  Roll  Lock-In,  which  means 
that  the  rolling  velocity  and  the  yaw  rate  are  the  same. 

One  effect  of  this  phenomenon  is  to  force  the  missile  to 
remain  in  a  condition  of  resonance  instability.  (Other  sit¬ 
uations  are  discussed  in  Reference  8* )  This  failure  of  the 
missile  to  roll  in  the  designed  manner  clearly  negates  any 
evaluation  of  its  performance  baaed  on  assumed  satisfactory 
roll  and  the  linear  theory.  Missiles  prone  to  Roll  Lock-In 
may  be  evaluated  using  the  methods  of  Reference  8. 

Case  6. -Roll  Speed-Up 

It  has  been  observed  that  additional  non-linear  roll 
moments  act  on  finned  missiles  in  flight.  One  of  the  effects 
of  these  neglected  moments  Is  to  cause  a  marked  increase  in 
the  steady  state  rolling  velocity  of  the  missile  as  a  func¬ 
tion  of  angle  of  attack.  This  Roll  Speed-up  phenomenon  can 
lead  to  poor  flight  performance.  For  example:  for  the  de¬ 
signed  roll  rate,  the  missile  might  be  dynamically  stable; 
however,  due  to  the  roll  speed-up  phenomenon,  the  increase 
in  rolling  velocity  with  angle  of  attach  may  yield  a  jUtagnus 
type  instability.  A  more  critical  case  would  be  a  combin¬ 
ation  of  this  roll  speed-up  with  Case  2  (non-linear  magnus 
instability).  This  type  of  instability  is  believed  to  be 
the  key  problem  in  the  short  ranges  of  mortar  shells. 

Case  7.  -  Catastrophic  Yaw 

In  addition  to  the  roll  moments  neglected  In  Case3  £ 
and  6,  side  force  and  side  moments  are  also  left  out  of  the 
linear  aeroballistic  theory.  This  neglected  side  moment, 
when  coupled  with  the  roll  lock-in  of  Case  5>,  can  yield  a 
magnus  type  non-linear  dynamic  instability.  In  Reference  8 
a  Modified  Quasi -Linear  Theory  is  suggested  for  approximately 
evaluating  this  phenomenon.  The  application  of  this  theory 
to  a  few  cases  appears  encouraging. 
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An  Application  of  Aeroballistics  Range  Techniques 
to  Models  with  Digonal  Rotational  Symmetry 

By 

G.E.  Tidy  &  Kiss  M.E.  Thomas 
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Establishment 

SUMMARY 

A  series  of  flat  plate  wings  of  triangular  planforra  have  been  fired 
at  Mach  Numbers  1.5  and  2.  A  few  trajectories  have  been  measured. 
Preliminary  manual  reduction  of  the  data  is  presented  and  the  derived 
values  of  some  aerodynamic  stability  coefficients  are  conqjared  with  linear 
theory.  The  possibility  of  more  complete  analysis  and  of  application  to 
airplane  configurations  are  considered  briefly. 

INTRODUCTION 

The  history  of  the  free  flight  range  as  a  technique  for  aeroballistics 
measurements  has  been  one  of  steady  increase  in  range  of  applicability  to 
aerodynamic  problems.  The  measurement  of  characteristics  of  spinning 
shells  and  of  winged  missiles  with  trigonal  or  greater  symmetry  has  been 
highly  developed  on  the  basis  of  linear  theory.  Recent  extensions  of 
reduction  methods  to  handle  certain  non  linear  aerodynamic  coefficients, 
within  the  same  limits  of  symmetry,  have  greatly  increased  the  scope  of  the 
technique . 

The  possibility  of  applying  aeroballistics  range  methods  to 
measurements  on  aerodynamic  configurations  falling  outside  this  category  of 
symmetry  is  of  great  interest.  This  particularly  applies  to  airplane  like 
shapes,  as  the  development  of  modern  supersonic  aircraft  have  increased  the 
need  for  experimentally  determined  dynamic  stability  derivatives. 
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Though  such  measurements  can  be  obtained  from  rocket  propelled  telemetered 
vehicles  and  though  wind  tunnel  dynamic  balance  methods  have  become 
available,  there  is  no  doubt  that  all  possible  techniques  should  be 
exploited.  The  aeroballistics  range  offers  certain  advantages  as  well,  in 
that  the  atmosphere  is  uniform  and  non  turbulent,  and  a  fairly  high  Reynold's 
Number  is  obtainable. 

With  the  expectation  that  extension  of  range  methods  to  configurations 
of  one  and  two  plane  symmetry  would  thus  be  of  considerable  value,  we  have 
begun  at  CAEDE  programmes  of  firings  of  aircraft  models  and  of  simple  flat 
plate  triangular  wings.  The  latter  program  forms  the  subject  of  this  report. 
This  study  was  instituted  to  explore  possible  methods  of  data  reduction  for 
mirror  symmetry  and  digonal  rotational  symmetry  configurations,  and  in  the 
event  of  a  successful  solution  of  the  problem,  to  obtain  experimental 
verification  of  some  theoretical  dynamic  derivatives  of  wings  alone,  free  of 
body  interference. 

This  programme  is  by  no  means  complete  and  important  problems  remain  to 
be  solved  before  a  successful  data  analysis  method  can  be  assured.  This 
paper  should  be  considered  as  a  progress  report  only,  on  a  continuing  programme 
of  study. 

Seventeen  wings  have  been  fired  to  the  present  date,  of  which  13  were 
of  aspect  ratio  one,  and  the  remainder  aspect  ratio  2.  The  majority  have  been 
fired  at  a  Mach  number  of  about  1.5,  although  two  have  been  launched  at 
Mach  2 . 

The  models  are  of  simple  construction,  as  shown  in  figure  1.  Centre  of 
gravity  adjustments  were  made  by  cutting  out  part  of  the  rear,  filling  the 
slots  with  balsa  wood,  and  attaching  aluminum  coverplates.  The  leading 
edges  were  wedge  shaped,  with  a  total  wedge  angle  of  about  50  degrees. 
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The  thickness  to  root  chord  rati©  w as  approximately  6  percent.  Little 
effort  was  expended  on  surface  finish  since  the  main  purpose  of  these 
preliminary  rounds  was  to  investigate  analysis  procedures ,  rather  than  to 
obtain  accurate  data.  Schlieren  photographs  of  a  model  in  flight  appear 
in  figure  3, 

Due  to  a  launching  problem  which  developed  a  large  number  of  these  rounds 
have  been  failures  in  that  only  very  small  deviations  from  a  straight  line 
trajectory  developed.  Insufficient  measurement  stations  were  employed  for 
some  of  the  earlier  models,  so  that  of  the  total  number  fired,  trajectory 
measurements  in  sufficient  detail  have  been  obtained  for  three  models  of 
aspect  ratio  one  at  Mach  1.5,  and  for  one  at  Mach  2.  Analysis  has  been 
concentrated  on  two  of  these  at  Mach  1.5. 

RANGE  TBCHNIQOE 

Before  further  discussion  of  the  present  problem,  a  brief  description 
of  the  instrumentation  employed  at  CARDE  is  in  order.  The  measurement 
apparatus  differs  from  photographic  ranges,  and  analysis  methods  must  be 
tailored  in  some  degree  to  suit  our  instrumentation. 

The  models  are  sabot  launched  into  the  range  from  a  5-9  inch  calibre 
smooth  bore  gun  fitted  with  a  hi-lo  device  for  pressure  equalization. 

Light  screens  at  50*  intervals  measure  velocity  and  deceleration.  Plane 
sheets  of  paper,  called  "Taw  Cards",  are  surveyed  into  the  range  at  5  foot 
intervals  along  the  firing  line.  Their  normals  are  set  parallel  to  the 
fore  and  aft  axis  of  the  range,  very  nearly  that  of  the  line  of  fire. 
Horizontal  and  vertical  reference  lines  are  surveyed  on  each  sheet. 


The  model  leaves  a  clear  impression  of  the  apex  and  the  base  on 

each  card.  (See  Figure  4) •  The  angle  of  roll  is  »easured  directly  as  the 

angle  between  the  base  line  and  the  vertical-  The  distance  from  the  centre 

of  the  base  line  to  the  apex  imprint,  divided  by  the  root  chord,  is  equal  to 

the  tangent  of  the  total  angle  of  attack.  The  total  angle  of  attack  can  be 

resolved  into  angles  of  pitch  and  sideslip  in  body  fixed  coordinates,  as 

used  for  this  report.  (See  Figure  5).  Since  v  ,  w  are  the  tangents  of 

V  V  ' 

the  card  imprints  measure  these  velocity  ratios  directly.  At  angles  of 
attack  greater  than  ten  degrees  or  so  small  projection  corrections  must  be 
applied.  Further  corrections  are  required  for  rotation  of  the  body  axps 
during  the  interval  between  passage  of  the  nose  and  of  the  base  through  the 
card.  Independent  measurements  of  angle  of  attack  agree  normally  to  wlttiin 
l/lO  degree  and  roll  to  about  l/4  degree. 
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The  position  of  the  centre  of  gravity  may  also  be  measured  with  respect 
to  the  survey  references.  It  has  not  been  possible  as  yet  to  set  the 
references  to  better  than  1/4  in.  on  a  majority  of  the  measurement  stations, 
so  that  data  on  CG  displacement  are  not  reliable  for  accurate  work. 
DESCRIPTION  OF  TRAJECTORIES 

The  motion  of  the  Wings  is  interesting  because  of  its  unsual  character. 

A 

Figure  6  shows  the  history  of  roll  rate  and  angles  of  attack  and  sideslip 
in  body  fixed  axes. 


This  is  a  typical  history  for  rounds  launched  at  zero  pitch.  It  is  seen  that 
a  gradual  increase  in  ^  occurs,  until  a  small  angle  oC  is  sufficient 
to  produce  an  appreciable  rolling  moment.  The  roll  rate  induced  is  such  as 
to  increase  angle  of  attack  and  produce  a  high  acceleration  in  roll  rate. 

The  roll  rate  oscillates  about  a  quite  well  defined  constant  value. 
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o(  and  ^  also  oscillate  about  an  equilibrium  value,  though  in 
general  the  steady  c(0  and  (30  cannot  be  so  accurately  determined  as  the 
steady  roll  rate,  because  of  the  complex  form  of  the  oscillation.  The 
equilibrium  motion  can  be  most  simply  described  in  terms  of  the  total  angular 
velocity  vector  which  tends  to  line  up  with  the  displacement  velocity  . 

The  steady  initial  increase  of  (Z  under  equilibrium  launch  conditions 
is  due  to  lack  of  strong  restoring  moments  in  sideslip,  in  fact  it  was 
because  of  the  unknown  effects  of  such  low  restoring  moments  that  the  models 
were  originally  launched  with  a  low  initial  angular  rates.  Though  our  first 
half  dozen  models  quite  readily  entered  the  rolling  condition,  it  was 
gradually  found  that  the  initial  condition  of  low  angle  of  attack,  yaw  and 
roll  rate  could  persist  through  most  of  the  instrumented  range,  the  resultant 
histories  being  useless  for  analysis.  The  launching  problem  which  has  thus 
arisen  is  being  investigated  in  two  ways  of  which  one  la  to  set  the  model 
at  an  initial  angle  of  attack  in  the  sabot.  This  method  appears  to  stabilize 
the  wing  in  a  dutch  roll  like  trajectory  and  appears  promising  for  obtaining 
aerodynamic  data  near  zero  angles  of  attack.  Initial  angles  of  attack  of  up 
to  10  degrees,  combined  with  yaw  are  to  be  tried  in  the  near  future.  The  other 
consists  in  attempting  to  induce  the  rolling  trajectory  by  causing  an  initial 
yaw  rate.  The  rolling  trajectory  appears  of  interest  because  measurements 
can  then  be  made  at  high  angles  of  attack. 

THE  EQUATIONS  OF  MOTION  FOR  A  WING  IN  FREE  FLIGHT 

We  may  now  pass  on  to  a  discussion  of  the  equations  of  motion,  in  order 
to  explain  the  observed  trajectories  and  to  form  a  mathematical  basis  for 
deriving  aerodynamic  coefficients. 
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normally  negligable,  can  be  as  great  as  U%  at  peak  angles  of  attack  reached 
in  these  experiments.  Only  those  aerodynamic  coefficients  were  retained 
which  were  thought  likely  to  be  of  measurable  size  in  the  preliminary  work. 
All  non  linearities  in  coefficients  except  in  rolling  moment  have  been 
ignored.  The  quantities  Cy^  ,  Qv  p  were  retained  in  the  hope  they 
could  be  measured  though  they  were  expected  to  be  very  small.  They  are 
caused  by  differential  friction  and  leading  edge  effects.  The  equations  can 
be  generalised  if  necessary  to  account  for  observable  non  linearities  in  the 
aerodynamics. 

A  simple  steady  state  solution  to  the  equations  is  possible  by  setting 
the  time  derivatives  equal  to  zero  and  ignoring  the  small  moments  and  forces 
due  to  sideslip. 
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The  equations  are  shown  here  in  a  form  suitable  for  deriving  aerodynamic 
coefficients.  The  prime  sign  refers  to  differentiation  with  distance  along 
the  trajectory.  The  solution  for  <j>a  shows  two  possible  values,  zero  and  one 
determined  by  as  written.  The  zero  roll  rate  condition  is  only 

neutrally  stable,  due  to  low  restoring  moments  in  sideslip. 

It  appears  that  If  the  rolling  trajectory  can  be  induced,  useful  data 
may  be  obtained  from  the  steady  state.  The  values  of  obtained  to  date 

appear  reasonably  accurate,  but  and  depend  upon  the  measurement 

of  which  is  usually  small,  therefore  these  data  may  not  prove  very 
reliable.  The  term  G^'^o  is  small  compared  with  C at  least  for 
the  wings  tested  to  date  and  has  been  Ignored. 
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SOLUTION  0 7  THE  GENERAL  EQUATIOMS 


When  a  solution  to  the  differential  equations  of  notion  is  known, 
as  in  the  case  of  shells  or  cruciform  winged  missiles,  a  first  approximation 
to  the  stability  derivatives  can  usually  be  found  by  measuring  overall 
characteristics  of  the  trajectories,  such  as  frequencies  and  damping  rates. 

If  required,  a  set  of  differential  corrections  may  be  obtained  in  some 
way  such  as  the  following  .  A  Taylor*8  Series  expansion  of  the  equation  of 
motion  in  c*.  ,  ^  etc.  is  obtained  in  terms  of  increments  of  the 
stability  derivatives  and  initial  conditions.  The  zero  order  term  of  the 
Taylors  series  is  the  solution  of  the  differential  equation  using  the 
approximate  values  of  the  aerodynamic a.  This  term  is  computed  for  each 
station  and  compared  with  the  observed  values. 


The  difference  between  observation  and  computation  equals  the  sum  of  the 
remainder  of  the  Taylor  aeries  terms.  A  set  of  linear  equation  in  the 
differential  corrections  is  thus  set  up,  and  a  least  squares  solution  can  be 
obtained . 
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When  one  turns  to  the  present  situation  an  analytic  solution  to  the 
equations  does  not  appear  very  likely.  Pitch  and  yaw  modes  are  intercoupled , 
even  in  the  general  motion  at  2ero  steady  roll  rate.  Moreover  non-linear 
terms  such  as  rolling  moment  due  to  sideslip  and  non-linear  dynamic  terms 
are  very  important. 

Therefore  we  plan  to  employ  a  machine  solution  of  the  differential 
equations.  Some  scheme  of  variation  of  initial  conditions  and  stability 
derivatives  to  produce  a  set  of  differential  corrections  is  contemplated, 
somewhat  as  in  the  method  Just  discussed  for  the  trigonal  symmetry  case, 
but  as  yet  very  little  time  has  been  available  to  investigate  this  phase  of 
the  analysis.  The  method  requires  an  approximate  set  of  values  for  the 
stability  derivatives  and  initial  conditions,  and  our  analysis  effort  up  to 
now  has  been  spent  on  checking  the  suitability  of  one  way  of  finding  such 
approximations . 

The  method  used  is  the  obvious  one  of  treating  the  differential 
equations  of  motion  as  linear  equations  in  the  stability  derivatives.  The 
functions  of  pitch,  sideslip,  roll,  and  their  derivatives  which  occur  in  the 
differential  equations  may  be  measured  at  each  station,  from  the  trajectories. 
A  set  of  linear  equations  in  the  stability  derivatives,  one  for  each 
measurement  station  is  thus  obtainable*  and  can  be  solved  by  least  squares. 

The  original  equations  contain  the  angular  velocities  q,  r  and  their 
derivatives,  and  these  quantities  are  not  accurately  measurable  with  the 
present  CARDS  instrumentation,  because  accurate  data  in  lateral  C.G.  displace¬ 
ment  are  required  for  their  determine tion.  The  equations  must  therefore  be 
combined  to  eliminate  q,  and  r,  three  equations  resulting.  We  call  these 
the  pitch,  yaw  and  roll  equations,  from  the  main  aerodynamic  derivatives 
contained  in  each. 


250 


THE  EQUATION  IN  PITCH 

The  equation  in  pitch  is  shown  below.  The  derivatives  now 
expressed  with  respect  to  distance  down  range,  rather  than  to  time. 


-  +  -(i 


etc:  ratio  of  the  radius  of  gyration  (about  x  axis)  to  reference 
*  length 


^  -  air  density 

S  s  reference  area  (wing  area) 

c  z  n  length  (root  chord) 

The  quantities  on  the  right  hand  side  contain  no  aerodynamic 
coefficients  and  may  be  lumped  together  as  the  "dynamic  variable".  If  the 
dynamic  variable  can  be  derived  together  with  the  "dynamic  factors"  e<  , 
C&!  and  so  on,  then  at  each  measurement  station  we  have  a  linear 
equation  in  >  the  damping  derivatives  and  the  small  quantity  Cyp 

The  group  of  equations  for  all  stations  are  then  solved  by  a  least  squares 
procedure. 
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EQUATIONS  IK  YAW 


The  equation  in  yaw  appears  below.  It  will  be  observed  that  for  the 
flying  wing,  the  largest  aerodynamic  derivative  is  C^. 


However  if  the  method  were  applied  to  the  general  airplane  motion,  yawing 
and  damping  derivatives  would  be  of  significance  and  measurable  by  means 
of  this  equation. 

EQUATION  IN  ROLL 

The  equation  in  roll  is  as  follows. 


Pi 

r 
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From  this  equation  one  can  derive  the  coefficient  of  rolling  moment  due  to 
sideslip,  and  rolling  moment  due  to  roll  rate. 


MBASIffiEMEHT  OF  FIRST  &  SECOND  DERIVATIVES 

First  and  second  derivatives,  and  the  functions  oC.  ,  ^  and  tj> 

themselves,  have  been  calculated  at  points  midway  between  stations  by 
expressing  them  in  terms  of  first,  second  and  third  differences  in  the 
observed  data.  Measurement  errors  are  thus  carried  down  to  the  dynamic 
variables  of  the  differential  equations.  The  accumulation  of  errors  due  to 
differencing  is  probably  the  chief  reason  why  this  method  yields  only 
approximate  results.  Application  of  data  smoothing  methods  would  result  in 
partial  elimination  of  significant  variation  in  the  trajectory  functions, 
with  the  present  number  of  data  points  per  oscillation  wave  length. 


LIMITATIONS  OF  THE  METHOD 


One  of  the  limitations  of  the  direct  application  of  this  method  of 

* 

analysis  is  that  indeterminate  values  of  aerodynamic  derivatives  are  liable 
to  be  obtained,  if  the  derivatives  have  dynamic  factors  which  are  nearly 


proportional  to  one  another.  For  instance  the  factors  of  Q|,  and  '“L/C- 
in  the  roll  equation  have  proven  to  be  proportional  for  the  rounds  measured 
to  date.  Fortunately  here  the  quantity  turns  out  to  be 

small  compared  with  and  can  be  obtained  by  some  other 

means,  so  that  at  least  is  derivable. 

RESULTS 


c 
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This  approximate  analysis  has  been  applied  to  two  rounds  of  aspect  ratio 
one,  at  approximately  Mach  1.45*  One  round  was  of  aluminum,  the  other  of 
steel,  the  ratio  of  weights  about  ?  to  1. 
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Both  assumed  the  steady  state  roll  mode,  the  lighter  model  rotating  at  about 
0.14  rad/ft  as  compared  -with  0.098  rad/ft.  for  the  steel  wing. 

Derivatives  calculated  on  the  basis  of  flat  plate  linear  theory  for  sero 
angle  of  attack  and  sideslip,  using  the  methods  of  NACA  Report  971  by 
Malvestuto  and  Margolis,  are  presented  for  comparison  (see  table  X). 

PITCHING  MOMENT  &  LIFT  DUE  TO  ANGLE  OF  ATTACK 

The  internal  agreement  between  dynamic  and  steady  Btate  measurements  is 
2%  in  both  cases  for  pitching  moment.  It  had  been  intended  to  vary  by 

changing  CG  position  in  order  to  find  Ce*.  the  centre  of  pressure  but 
all  rounds  with  forward  CG  were  •'’allures.  These  wings  had  CG's  at  60.4  and 
6l.4$  root  chord,  -  too  near  one  another  to  give  good  results  for  cp  position. 
The  cp  as  calculated  is  65%  root  chord  but  a  change  of  one  percent  in  both 
rounds  can  alter  the  cp  to  67%.  On  this  basis  the  coranon  for  both 

rounds  is  as  shown.  The  result  is  somewhat  high  compared  with  linear  theory. 
The  values  of  calculated  from  the  yaw  equation  show  variations  which 

may  be  due  to  omission  of  terms  in  the  equation  for  asymmetry  moments 
caused  by  model  roughness  and  errors  in  construction.  The  steady  state  values 
of  Cioi  are  of  poor  dependability  since  the  steady  sideslip  angle  is  small 
compared  with  its  variation. 

ROLLING  .MOMENT  DUE  TO  YAW 

There  is  no  steady  state  solution  for  this  quantity.  The  roll  equation 
has  been  solved,  by  using  the  linear  theory  value  for  Ci*  ,  appearing  in  the 
equation. 


From  Steady  State  Equations  D  -  From  Dynamic  Equations 
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TABLE  I  -  STABILITY  DERIVATIVES 
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The  resulting  term  is  so  snail  that  the  inaccuracy  of  this  value  will  not 
appreciably  affect  the  result.  The  individual  rounds  differ  from  the  average 
by  about  \L,%  and  the  average  is  1056  lower  than  linear  theory  would  indicate. 
DAMPING  COEFFICIENTS  IK  FITCH 

A  large  variation  is  evident.  Even  signs  are  wrong  in  some  cases.  One 
possible  reason  for  this  result  is  that  the  dynamic  factor  multiplying  ( C-t,^ 

is  nearly  180°  out  of  phase  with  and  proportional  to 

the  dynamic  factor  of  Cm<^-  pitch  equation,  over  most  of  the 

trajectories.  The  results  may  therefore  be  somewhat  indeterminate  and  the 
coefficients  subject  to  considerable  error  • 

DAMPING  IN  ROLL 

The  dynamic  damping  in  roll  values  agree  fairly  well  with  each  other 
and  with  theory.  However  the  steady  state  results  are  unreliable,  again  due 
to  the  difficulty  of  measuring  steady  state  sideslip. 

No  correlation  of  the  sideslip  dynamic  variable  with  ^  has  been  found, 
indicating  that  the  effects  of  yawing  moment  due  to  sideslip  are  less  than 
experimental  errors. 

CONCLUSIONS 

The  quantity  of  data  obtained  to  date  is  too  scanty  to  determine  the 
accuracy  of  the  approximate  results.  In  any  case  the  final  test  must  be  the 
application  of  the  data  to  a  computer  solution  for  differential  coefficients. 
However  the  results  at  least  for  the  larger  stability  derivatives  appear 
promising  enough  to  warrant  further  work  on  the  problem. 

A  difficulty  has  arisen  in  the  indeterminate  solutions  for  damping 
derivatives  in  pitch,  which  will  apply  to  any  machine  solution  as  well.  One 
possible  means  of  overcoming  this  indeterminacy  is  to  fire  the  rounds  with 
high  initial  angles  of  attack  and  the  minimum  of  sideslip  angle  and  rate,  so 
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that  the  early  part  of  the  trajectory  will  consist  of  oscillations  in  pitch 
with  low  roll  rate.  The  damping  term  in  rate  of  change  of  angle  of  attack 
will  then  be  predominant. 

The  present  results  do  not  permit  one  to  predict  with  confidence 
whether  the  analysis  method  will  be  applicable  to  airplane-like,  low  roll 
rate  trajectories.  This  type  of  trajectory  will  produce  lower  peak 
amplitudes  in  total  angle  of  attack  resulting  in  higher  percentage  errors 
associated  with  the  dynamic  variables.  However  the  approximate  stability 
derivatives  so  derived  may  still  be  sufficiently  accurate  for 
application  to  differential  correction  methods.  It  is  hoped  that  further 
teats  with  emphasis  on  trajectories  with  low  roll  rate,  will  help  to  clarify 
this  situation. 
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SABOTS  USED  AT  THE  THOMPSON  AERCBALLISTICS  LABORATORY 


Wallace  H.  Allan 

U.  S.  Naval  Ordnance  Test  Station 

INTRODUCTION 

Aeroballistic  range  work  requires  some  means  for  projecting  the  missile 
through  the  range;  the  accepted  standard  for  this  work  is  a  gun.  The  use  of 
a  gun  entails  the  use  of  a  sabot  in  almost  all  cases,  the  principal  excep¬ 
tions  being  projectiles  designed  for  guns.  The  sabot  is  usually  a  piston  that 
holds  and  pushes  the  missile  and  seals  the  bore  against  the  propelling  gases. 
Sabots  used  at  the  Thompson  Aeroballi sties  Laboratory  fall  into  two  classes: 
those  intended  for  smooth-bore  guns,  usually  referred  to  as  f inner  sabots, 
and  those  intended  for  rifled  guns,  referred  to  as  spinner  sabots.  These  two 
classes  fall  into  two  more  classifications:  those  that  "break  apart"  after 
leaving  the  gun,  and  those  that  "follow  thru"  the  range  after  the  missile; 
and,  of  course,  there  are  those  that  do  a  little  bit  of  both. 

The  Thompson  Aeroballistics  Laboratory  has  no  systematic  program  of  sabot 
development.  Years  ago  we  did,  principally  to  give  us  a  feeling  for  the 
problems  involved.  Now,  each  missile  and  sabot  combination  is  treated  as  a 
separate  problem.  The  sabot  is  designed  on  the  basis  of.  past  experience  and 
the  first  attempt  at  the  design  is  usually  fired  through  the  Laboratory  for 
data,  although  an  outdoor  range  is  available  for  sabot  development  work.  If 
the  first  firings  are  successful,  then  the  program  can  be  couple  ted  in  the 
minimum  of  time  and  expense.  If  the  firing  is  unsuccessful  then  the  problem 
is  reappraised,  the  sabot  redesigned,  the  laboratory  repaired,  and  the  program 
resumed.  The  cost  of  sabots  and  models,  firing  time,  and  the  time  scale  of 
the  program  lead  us  to  believe  that  this  is  the  inexpensive  way  to  conduct 
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our  programs.  Damage  due  to  errant  rounds  has  been  small  compared  to  the  cost 
of  good  rounds  that  would  have  been  lost  under  a  more  conservative  program. 

SPINNER  SABOTS 

Let  us  first  consider  spinning  projectiles.  We  have  available  the  following 
twist  barrels:  40mm  1-30  twist;  3"/50  1-10,  1-15,  1-20,  1-25,  and  1-32  twist; 
37rO  1-22  twist;  5 ”/5^  1-10,  and  1-25  twist.  The  3"  diameter  spinning  projectile 
is  considered  optimum  for  our  range.  40mm  projectiles  are  usually  too  short 
for  accurate  orientation  readings  and  the  5"  projectiles  react  so  slowly  to 
external  forces  that  only  limited  portions  of  the  damping  arms  can  be  observed. 

Many  of  the  spinning  projectiles  fired  through  our  laboratory  do  not  re¬ 
quire  sabots  since  they  are  designed  for  gun  launching .  For  instance.  Fig.  1 
is  a  picture  of  a  full-scale  5 '7 54  projectile  and  two  subcaliber  models  of  this 
projectile.  The  models  can  be  fired  from  the  3”/50  and  40mm  guns,  and  of  course 
the  5"  projectile  is  fired  from  the  standard  5**/ 5^  gun.  In  the  case  of  the 
40ram  model,  our  only  available  barrel  for  the  40mm  gun  did  not  have  enough  twist 
to  stabilize  the  model  by  spin;  therefore,  a  subcaliber  sabot  was  designed  to 
allow  firing  the  40ram  model  from  a  3"  1-10  twist  barrel  giving  the  barrel  an 
effective  twist  of  1-18.  This  is  a  very  effective  way  of  obtaining  odd  twist 
ratios  not  readily  available  in  our  stock  of  gun  barrels.  Figure  2  illustrates 
this  type  of  subcaliber  sabot.  This  particular  sabot  consists  of  a  driving 
piston  containing  a  pre-engraved  rotating  band  and  an  expanded  foam-plastic 
body  (styrofoam)  to  support  the  round  before  firing.  The  base  of  the  missile 
is  fitted  with  pins  used  to  observe  by  photographic  means  the  spin  rate  of  the 
missile.  These  pins  fit  into  holes  in  the  piston  portion  of  the  sabot,  and 
the  torque  to  spin  the  missile  is  applied  through  the  pins.  Not  all  of  the 
torque  is  applied  through  the  pins,  since  the  coefficient  of  friction  between 
the  base  of  the  missile  and  the  piston  is  almost  sufficient  to  drive  the  missile 
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In  torque  In  most  cases.  Spin  and  air  drag  separates  the  foam  plastic  and 
piston  from  the  missile  after  leaving  the  gun.  This  type  of  sabot  is  also 
used  when  high  velocity  firings  are  required,  since  many  bore  size  rocket 
models  are  very  heavy,  limiting  the  ma-vimran  velocity  at  which  they  can  be 
fired.  This  type  of  sabot  has  been  useful  but  is  not  100$  reliable.  Often 
missiles  are  mallaunched  with  this  sabot.  Our  next  attempt  at  this  type  of 
sabot  will  eliminate  the  foam  plastic  and  replace  it  with  nylon  plastic. 
Subcaliber  sabots  are  frequently  driven  at  excess  spin  rates.  Quite  often 
the  copper  obturating  band  or  the  al  umlnum  or  copper  rotating  bands  burst 
under  the  high  spin  rates.  In  this  case  the  sabot  parts  are  caught  by  a 
blast  shield  in  front  of  the  gun.  Figure  3  illustrates  this  shield.  The 
large  circle  of  impact  marks  was  made  by  the  driving  bands  fired  from  the  3" 

1-10  twist  barrel.  The  circle  diameter  is  independent  of  velocity  but  depends 
mainly  upon  the  twist  ratio  of  the  gun.  If  a  sabot  is  apt  to  break  apart 
because  of  high  spin,  then  the  driving  band  is  deliberately  weakened  by 
partially  cutting  the  rotating  band.  Care  has  to  be  taken  that  the  driving 
pins  on  the  base  of  the  missile  are  not  damaged  by  the  breaking  of  rotating 
bands.  The  inner  portion  of  the  sabot  is  usually  made  of  75  ST  aluminum  and 
it  usually  "follows  thru"  the  laboratory  after  the  missile.  Figure  4  illus¬ 
trates  another  sabot  of  this  form.  This  particular  sabot  was  designed  by 
the  Aerojet-General  Corporation.  In  this  case  the  missile  is  held  to  the 
piston  by  three  pins  that  fly  out  under  centrifugal  force,  and  the  barrel 
prevents  the  pins  from  withdrawing  until  the  missile  leaves  the  barrel.  The 
front  support  ring  riding  on  the  ogive  is  cut  almost  in  half  so  that  centrifugal 
force  will  break  it  apart.  In  fact,  there  is  only  .020  of  an  inch  of  metal 
on  the  inside  diameter  holding  it  together.  The  sabot  is  blown  from  the 
missile  by  trapping  propelling  gases  in  a  chamber  between  the  sabot  and  missile. 
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An  earlier  version  of  this  sabot  had  a  copper  obturating  ring  located  in  the 
forward  groove  and  the  front  support  ring  was  held  together  on  the  outside 
diameter  by  l/l6  inch  of  metal.  Figure  5  Is  a  synchro-ballistic  picture  of 
the  earlier  version  of  this  sabot  in  action.  First  ve  see  that  the  front 
support  did  not  break  apart.  Next  the  forward  obturating  ring  expanded  under 
centrifugal  force  to  a  point  where  it  was  not  discarded  when  the  sabot  blew 
back.  Corrections  made  to  this  combination  were  simple  and  effective.  The 
forward  seal  was  discarded  completely  and  the  slot  in  the  front  support  was 
cut  to  only  .020  of  an  inch.  The  sabot  worked  fairly  well  then. 

Bore-diameter- spinning  rocket  models  generally  need  a  sabot  to  impart  the 
correct  spin  to  the  missile.  Figure  6  illustrates  a  piston  sabot  we  have  used 
for  many  years.  In  this  sabot  the  center  piston  fits  in  a  cylindrical  hole  in 
the  base  of  the  round.  A  fana/ll  hole  through  the  sabot  permits  pressure  to 
build  up  in  a  gas  chamber  located  at  the  end  of  the  piston.  Upon  leaving  the 
barrel  the  gas  reacts  against  the  piston  and  blows  the  sabot  free  of  the 
missile.  The  small  gas  hole  in  the  sabot  is  drilled  through  a  steel  insert 
to  prevent  gas  erosion  from  enlarging  and  filling  the  gas  chamber  with  molten 
metal  which  tends  to  jam  the  piston  in  the  hole.  The  driving  torque  is  applied 
from,  the  rotating  band  through  the  spin  pins  that  are  used  for  spin  measurements. 
Copper  rotating  bands  are  preferred  because  of  their  very  effective  obturating 
actions;  however,  the  spin  rates  are  usually  too  high  for  copper  so  that  75  ST 
aluminum  rotating  bands  are  used  with  a  small  copper  or  nylon  obturating  band. 
Another  version  of  this  type  of  sabot  is  illustrated  in  Fig.  7  •  Here  spherical 
balls  are  used  to  drive  the  missile  in  spin  and  then  upon  leaving  the  barrel 
centrifugal  force  acts  upon  the  balls  wedging  the  round  and  sabot  apart.  The 
balls  are  caught  by  the  blast  shield  in  front  of  the  gun.  The  principal 
advantage  of  this  sabot  over  the  earlier  style  is  that  a  center  pin  may  be 
used  in  the  tail  of  the  missile  without  interference  with  the  blow-off  system. 
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This  type  of  sabot  is  less  susceptible  to  sticking  than  the  piston  sabots.  A 
variation  of  this  system  is  displayed  in  Fig.  8.  This  system  uses  discs  instead 
of  balls.  There  seems  to  be  no  essential  difference  between  the  two  schemes  in 

performance.  Figure  9  shows  a  yaw- inducing  spinner  sabot.  Yaw  cam  be  induced 

\ 

in  these  sabots  by  stepping  the  base.  The  mass  eccentricity  and  the  reaction 
of  the  muzzle  blast  on  the  step  both  serve  to  induce  yaw. 

In  spinner  sabots  one  major  consideration  is  whether  a  pre-engraved 
rotating  band  should  be  used.  I  personally  prefer  this  style,  but  is  is 
expensive  to  produce.  At  high  pressures  there  seems  to  be  no  appreciable 
difference  between  the  pre-engraved  and  engraving  style  of  sabot;  however,  at 
low  pressures  the  action  of  the  engraving  type  cam  be  erratic  and  in  many  cases 
the  projectile  has  stuck  in  the  barrel  only  partially  engraved.  One  major 

problem  which  has  arisen  in  the  use  of  engraving-type  sabots  is  the  separation 

/ 

of  the  sabot  and  the  missile  in  the  barrel.  This  is  caused  by  the  sabot 
accelerating  the  missile  for  a  very  short  distance,  contacting  the  engraving 
band,  then  stopping  for  an  instant  while  the  pressure  builds  up  to  the  engraving 
pressure.  The  projectile  during  this  instant  separates  from  the  sabot  and  later 
when  the  sabot  contacts  the  missile  again  it  will  frequently  damage  the  base 
of  the  projectile  or  the  spin  pins.  An  effective  way  of  eliminating  this  dif¬ 
ficulty  is  to  pin  the  missile  and  sabot  together  as  previously  illustrated. 

Our  usual  way  of  overcoming  this  difficulty  is  to  take  great  care  that  the 
sabot  is  properly  seated  in  the  barrel  with  no  clearance  between  the  engraving 
band  and  the  lands  of  the  barrel.  Figure  10  illustrates  another  form  of 
spinner  sabot.  This  sabot  was  designed  by  Aerojet-General  Corporation.  The 
missile  fits  the  3”  gun  and  the  pre-engraved  band  on  the  missile  fits  the 
gun  rifling.  This  rotating  band  would  not  take  gun  accelerations  unassisted. 

The  sabot  actually  furnishes  the  main  driving  torque  through  its  frictional 
contact  with  the  base.  The  sabot  is  built  so  that  any  slippage  between  the 
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sabot  and  the  missile  does  not  harm  the  spin  pins.  This  sabot  and  model 

combination  was  quite  successful. 

/ 

FEWER  SABOTS 

Let  us  now  consider  finned  projectiles.  We  have  available  the  following 
smooth-bore  barrels:  lVl-208,  42mm,  3 "1-46,  3V2-65,  5'.'l-37,  and  8V25-53-  The 
l'.'l  and  42mm  guns  are  considered  undersize  for  our  range.  Most  missiles  are 
fired  through  the  3"  or  5"  gun,  with  the  8"  gun  rapidly  becoming  very  popular. 

We  have  used  a  standard  style  of  ’  "follow- thru"  sabot  for  years.  Figure  11 
shows  the  general  configuration.  In  this  sabot  the  missile  is  held  by  fingers 
in  correct  alignment  with  the  bore  before  firing  and  under  acceleration.  The 
sabot  is  very  stable  in  its  flight  and  usually  is  more  accurate  than  the  missile 
The  holes  in  the  tail  of  the  sabot  not  only  serve  to  lighten  the  tail  but  they 
prevent  the  tail  from  shearing  in  the  bore  under  the  expansion  due  to  chamber 
pressure.  The  obturating  ring  is  quite  far  forward  in  this  design.  Yaw  can  be 
induced  with  this  sabot  by  leaving  a  few  of  the  gas-relief  holes  located  near 
the  tail  undrilled.  The  muzzle  pressure  can  then  give  the  tail  a  kick  as  the 
missile  leaves  the  muzzle.  This  sabot  is  somewhat  expensive  to  make  but  very 
reliable  in  action  in  guns  of  40mm  and  3"»  5"  sabots  in  this  design  have  failed 

quite  badly.  The  after  tail  section  usually  blows  up  completely.  If  there  is 
a  complete  blowup,  the  rocket  model  is  usually  launched  fairly  well;  but  many 
times  the  blowup  is  not  complete  and  the  missile  is  mallaunched.  Our  first 
attempt  at  curing  this  mallaunch  was  a  completely  new  sabot  called  the  slug- 
styrofoam  follow-thru.  Figure  12  illustrates  this  sabot.  This  sabot  is  the 
simplest  we  have  ever  made.  The  slug  is  usually  made  of  aluminum  for  use  in  5” 
guns  and  micart a- aluminum  for  the  8"  gun.  The  missile  sits  on  its  tail  under 
acceleration  and  the  styrofoam  is  used  only  to  support  the  missile  in  bore 
alignment  before  firing.  The  styrofoam  blows  apart  upon  leaving  the  muzzle 


and  the  slug  separates  from  the  missile  because  of  its  greater  aerodynamic 
drag.  This  sabot  has  launched  some  missiles  very  well  and  it  has  apparently 
failed  in  other  cases.  At  the  present  time  we  hesitate  to  recommend  it. 

Figure  13  shows  a  slug  follow- thru  which  has  worked  quite  successfully  in  the 
case  of  a  missile  which  failed  to  fly  well  with  the  slug-styrofoam,  sabot.  This 
sabot  is  essentially  the  front  half  of  the  original  follow- thru.  It  looks  quite 
a  bit  out  of  a  balance  but  remember  that  the  sabot  is  of  8"  diameter  and  the 
model  only  weighs  4.5  lbs. 

The  original  break-apart  sabot  used  in  our  laboratory  is  shown  in  Fig.  14. 
This  style  of  'break  apart"  has  been  used  at  up  to  3"  diameter,  but  has  caused 
barrel  scoring  at  velocities  of  over  Mach  3*  This  sabot  is  expensive  to  make 
and  it  is  moderately  heavy  when  made  of  aluminum.  Figure  15  is  an  example  of  a 
break-apart  sabot  for  a  cone  fired  from  the  40mm  gun.  This  sabot  was  successful 
at  Mach  3.  Figure  16  shows  a  break-apart  sabot  for  a  rather  odd  object.  The 
sabot  is  for  a  5"  gun  and  is  made  of  nylon  with  a  rubber  obturating  disc.  The 
sabot  is  sawed  almost  completely  into  quarters  and  just  a  few  tabs  are  left  to 
hold  it  together  during  the  loading  process.  This  particular  design  feature 
was  shown  to  us  by  Dr.  Bull  of  CARDE.  This  sabot  has  been  quite  successful. 
Figure  17  shows  another  variation  of  this  style.  This  sabot  is  for  the  3"  gun. 
Figure  18  illustrates  the  design  in  which  a  base  plate  was  felt  to  be  necessary. 
This  is  a  sabot  for  a  5"  gun  and  the  ball  is  a  solid  2. 5 -Inch-diameter  sphere 
made  of  steel.  The  spherical  surface  prevents  the  use  of  a  sabot  which  com¬ 
pletely  breaks  apart,  since  under  acceleration  the  spherical  surface  would  tend 
to  wedge  the  sabot  quarters  in  the  barrel.  This  sabot  breaks  into  four  quarters 
and  a  base  plate.  We  were  able  to  drive  these  balls  at  up  to  Mach  5  with  this 
type  of  sabot.  The  acceleration  was  so  great  that  an  early  model  of  this  sabot 
failed  because  the  ball  was  driven  completely  through  the  aluminum  base  of  the 
sabot  during  the  initial  acceleration  phase. 
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Figure  19  illustrates  a  rather  unusual  sabot  for  an  unusual  model.  The 
model  has  a  tungsten  head  and  a  body  consisting  of  aluminum  tubing  .010  of  an 
inch  thick.  The  body  is  quite  fragile,  making  it  necessary  to  push  the  model 
at  the  base  of  the  tungsten  head.  The  sabot  therefore  has  a  rather  long  rod 
in  order  to  push  the  missile.  Unfortunately  the  missile  has  a  rather  high 
drag  so  that  the  sabot  cannot  be  expected  to  pull  free;  therefore,  the  sabot 
was  pulled  from  the  missile  by  a  device  that  is  called  the  sabot  stopper. 

Figure  20  shows  this  device.  A  tapered  steel  tube  is  fastened  onto  a  muzzle 
attachment  that  aligns  the  tube  with  the  bore.  Clearance  between  the  muzzle 
and  the  entrance  to  the  tube  is  maintained  so  that  gases  in  the  barrel  may 
escape.  The  tapered  tube  is  smaller  than  the  sabot  diameter  so  the  sabot 
is  retarded  passing  through  the  tube.  The  missile  is  of  a  still  smaller 
diameter  and  passes  through  the  "sabot  stopper"  unhindered,  the  retardation 
of  the  sabot  separating  it  from  the  missile.  In  this  particular  case  no 
attempt  was  made  to  capture  the  sabot  but  only  to  slow  it  down.  In  other 
cases,  complete  captures  have  been  made  at  velocities  of  the  order  of  1000  ft/sec. 
The  steel  tube  is  given  a  very  low  velocity  during  the  retardation  process  and 
it  pulls  free  of  the  gun.  The  sabot  in  this  particular  case  proved  to  fly 
stably  through  the  range.  Another  interesting  sabot  is  illustrated  in  Fig.  21. 
This  sabot  is  driven  in  spin  by  a  turbine  drive.  Escaping  gases  react  to 
impart  spin  to  the  sabot.  Mr.  Hitchcock  of  APG  wrote  a  paper  on  this  technique 
with  reference  to  imparting  spin  to  a  mortar  shell.  We  fired  four  models  of 
this  type  of  sabot  and  found  that  we  could  inpart  slow  spin  to  the  missile. 

Twist  rates  of  1-200  were  achieved  in  the  case  of  the  many- fluted  model.  The 
slugs  were  never  adapted  to  any  actual  model  studies  by  us. 

We  use  the  3"  j  and  8"  smooth-bore  guns  often  in  our  work.  The  5"  is 
the  favorite  gun  because  the  missiles  are  about  optimum  for  accurate  work  in 
our  laboratory.  The  8"  gun  is  used  particularly  for  problems  involving  models 
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that  are  structurally  weak.  The  great  length  of  the  barrel  permits  low  accelera¬ 
tions.  The  5"  and-  especially  the  8"  gun  are  regarded  as  one  of  our  most 
important  tools  allowing  us  to  fire  programs  considered  impractical  elsewhere. 
Sabot  design  is  considered  an  important  factor  in  the  successful  completion 
of  a  program.  Unfortunately,  the  design  of  a  sabot  is  still  in  the  stages  of 
an  art  rather  than  a  science  because  the  designer  has  to  work  with  small 
numbers  of  widely  varying  missile  configurations  usually  under  pressure  to 
make  the  first  round  fly  properly. 
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FIG.  2.  3 -INCH  SUBCALIBER  SABOT  FOR  4-Omra.  MODEL 


FIG.  3.  BLAST  SHIELD  LOCATED  NEAR  GUN  N 
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FIG.  5.  SYNCHRO-BALLISTIC  PICTURE  OF  MODEL  AND  SUBCALIBER  SABOT 
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FIG.  7.  BALL  SABOT  FOR  SPINNERS 
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FRICTION-DRIVE  SABOT  FOR  SPINNERS 


FIG.  11.  FOLLOW -THRU  SABOT  FOR  FIMERS 
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BREAK-APART  SABOT  FOR  BALLS 


FIG,  19.  FOLLOW -THRU  SABOT  FOR  RODS 
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FIG.  20.  SABOT  STOPPER  ATTACHMENT  AND  TUBE 
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Joseph  E.  Long 
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MODEL  LAUNCHING  TECHNIQUES  AND  OTHER  ITEMS 
RELATED  TO  RANGE  FIRINGS 

In  "one  early  days  of  the  NOL  ranges  the  techniques  used  in  launching 
models  were  quite  simple.  Finned  missiles  were  fired  from  smoothbore  guns 
using  simply-constructed  sabots.  Spinning  missiles  were  equipped  with  more 
or  less  conventional  rotating  bands  and  were  fired  from  rifled  guns.  Many 
range  techniques  have  recently  been  developed  to  aid  in  the  solution  of 
aeroballistic  problems.  Some  of  these  developments  vill  be  described  in 
this  paper. 

CONSTANT  SPIN  OF  A  FINNED  MISSILE 

When  finned  missiles  are  fired  from  smoothbore  guns  with  little  or  no 
roll  it  is  usual  to  determine  the  drag,  restoring  moment,  normal  force,  and 
damping  moment  coefficients.  An  accidental  cant  given  to  the  fins  of  a 
missile  during  manufacture  may  cause  a  model  to  roll  at  a  sufficient  rate 
to  produce  lunar  motion.  To  avoid  this  resonant  condition,  the  fins  may  be 
made  with  a  larger  cant.  This  increases  the  spin  and  makes  it  imperative  that 
the  Magnus  moment  coefficient  be  measured.  To  obtain  this  coefficient  a 
finned  model  is  madAto  travel  down  range  with  constant  spin  after  it  is  fired 
from  a  Jdfled  gun.  The  fins  of  the  model  are  canted  to  maintain  the  same 
spin  produced  by  the  rifling  of  the  gun.  A  four-finger  peel-away  sabot  is 
used  to  launch  the  model  as  shown  in  Figure  1.  The  inside  of  the  sabot  is 
contoured  so  that  the  model  will  fit  and  bear  along  each  finger.  The  outside 
of  the  sabot  is  pre-engraved  along  the  entire  length  of  the  fingers  and  the 
sabot  base,  the  engraving  matching  the  rifling  of  the  gun.  The  model  and 
sabot, must  fit  in  the  gun  precisely.  Otherwise  the  tail  section  of  the  model 
would  ■unwind  or  the  fins  would  bend.  Figure  2  shows  two  recovered  rounds. 

It  can  be  seen  that  there  was  no  fin  damage.  This  same  technique  has  been 
applied  in  guns  up  to  three  inches  in  diameter. 

LAUNCHING  SUBCALIBER  PROJECTILES 

It  has  been  customary  to  determine  the  drag  coefficient  by  firing 
scaled  models,  but  the  scaling  has  not  included  the  details  of  the  rifling 
marks  left  on  the  rotating  band  of  the  full-scale  projectile.  Instead,  the 
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models  were  manufactured  with  oversized  rotating  bands  and  fired  from  guns 
with  a  smaller  number  of  rifling  grooves.  It  was  subsequently  found  that 
the  drag  coefficient  differed  by  about  three  percent  between  models  fired 
with  an  oversized  rotating  band  and  models  fired  with  scaled  grooving  and 
rotating  bands.  Recently,  all  grooving  and  rotating  band  details  have  been 
copied  from  recovered  full-scale  rounds.  These  reproduced  details  on  a  sub¬ 
caliber  model  are  preserved  in  a  sabot  as  shown  in  Figure  3.  There  are  two 
pins  in  the  base  of  the  sabot,  which  engage  holes  in  the  model.  The  outside 
of  the  sabot  is  made  larger  than  the  rifling  of  the  gun.  A  section  of  the 
sabot  just  aft  of  the  position  of  the  model  base  is  weakened  by  reducing 
the  sabot  diameter,  the  thickness  of  this  section  depending  on  the  velocity 
to  be  attained.  At  firing  it  is  believed  that  the  part  of  the  sabot  in  front 
of  the  weak  section  breaks  up  in  the  gun.  As  the  sabot  base  and  model  leave 
the  gun  they  begin  to  separate  as  shown  in  Figure  4jand  the  model  is  free  to 
fly  down  the  range.  Through  this  technique  a  more  realistic  drag  curve  can 
be  obtained. 

LAUNCHING  SPHERES 

In  a  recent  program  very  small  spheres  were  fired  for  drag  determination. 
It  was  realized  that  they  were  too  smal 1  to  actuate  the  sparks  by  photoelectric 
triggering.  To  overcome  this  deficiency  it  was  decided  to  fire  a  larger, 
triggering  sphere  with  the  smaller  one.  Acutally  several  sizes  of  spheres 
were  fired  together.  Figure  5  shows  l/52,  l/l6,  3/32,  and  l/8  inch  diameter 
spheres  that  are  loaded  into  either  a  synthetic  rubber  cup  for  low  velocity 
or  a  nylon  cup  for  high  velocity  firings.  These  spheres  are  held  in  the  cup 
by  a  piece  of  lens  tissue.  This  assembly  is  loaded  into  a  30-caliber  cartridge 
case  and  fired  in  the  range.  It  was  discovered  later  in  the  program  that  the 
large  diameter  spheres  were  not  necessary  as  the  cups  containing  the  spheres 
would  travel  down  the  range  and  trigger  the  sparks.  By  varying  the  air 
pressure  in  the  range  and  shooting  at  different  velocities,  drag  coefficients 
can  be  obtained  for  a  wide  variety  of  Reynolds  and  Mach  numbers. 
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LAUNCHING  MODEL  AIRCRAFT 

A  great  deal  of  publicity  has  been  given  recently  to  the  "sonic  booms" 
produced  by  high-speed  aircraft.  Speculations  have  been  made  as  to  the 
number  of  booms  that  one  might  expect  to  hear.  A  simple  experiment  vas 
prepared  in  the  Aerodynamics  Range  to  study  this.  A  model  airplane  (obtain©! 
from  a  box  of  cereal)  was  put  into  a  40-mm  sabot  consisting  of  two  spring 
steel  fingers,  styrofoam,  and  a  cylindrical  sabot  base  as  shown  in  Figure  6. 

At  launching  the  spring  steel  fingers  would  pull  the  styrofoam,  away  frcm  the 
model  allowing  the  model  to  fly  down  range .  The  shadowgraph  in  Figure  7 
shows  the  shock  pattern  obtained  from,  this  model.  Frcm  this  pattern  it  is 
interesting  to  conjecture  as  to  the  number  of  booms  that  might  he  heard  under 
various  conditions  of  observation. 

JET  MODEL 

Not  all  problems  in  the  ranges  are  launching  problems.  Other  difficulties 
arise  in  trying  to  make  a  model  or  seme  of  its  components  do  a  particular  Job. 
For  example,  a  program  has  been  instituted  to  investigate  what  effects  a  Jet 
exhausting  frcm  the  base  of  a  missile  can  have  on  the  aerodynamic  coefficients. 
Some  means  had  to  be  employed  whereby  the  shadowgraph! c  plates  would  not  be 
exposed  to  a  flaming  Jet.  For  an  initial  study  a  "Sparklet"  bulb  containing 
liquid  CO2  at  835  psi  gauge  was  placed  in  a  cone  cylinder  model  as  shown  in 
Figure  8.  A  thin  steel  plate,  not  shown,  was  put  on  the  base  of  the  model  to 
prevent  the  cap  from  popping  out  due  to  setback  forces  and  allowing  the  CO^ 
to  escape  at  the  time  of  firing.  A  hole  was  drilled  in  the  center  of  the 
plate.  To  puncture  a  hole  in  the  cap  of  the  bulb  a  pin  was  fastened  to  the 
center  of  a  disk,  and  the  disk  was  placed  Just  behind  the  plate  in  the  gun. 

At  firing,  the  pin  would  puncture  the  bulb  allowing  the  C02  to  escape.  A 
shadowgraph  of  a  model  with  a  Jet  of  COg  is  shown  in  Figure  9*  These  models 
were  fired  at  a  Mach  number  of  0.5.  Preliminary  results  indicate  a  considerable 
change  in  the  drag  and  overturning  moment  coefficients  due  to  the  cold  Jets. 
Encouraged  by  these  results  it  is  planned  to  design  a  finned  model  with  a 
nozzle  and  compressed  gas  supply  for  launching  in  the  firing  ranges. 
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FUZE  AIMING 


The  ranges  can  perform  a  variety  of  tasks  not  necessarily  associated 
with  obtaining  aerodynamic  coefficients  for  a  specific  model.  A  30-mm  model 
was  fired  at  2800  ft/sec  in  the  Pressurized  Yaw  Card  Range  under  a  pressure 
simulating  an  altitude  of  40,000  feet  specifically  to  determine  whether  or 
not  the  fuze  armed  itself  under  actual  flight  conditions.  It  was  caught  in 
cotton  waste  to  prevent  its  destruction.  The  arming  device  consisted  of  a 
ball  with  a  cylindrical  hole  containing  a  powder  train.  In  flight  the  ball 
tends  to  rotate  about  the  axis  for  which  the  moment  of  inertia  is  a  maximum.. 
Upon  contact  with  a  foreign  object,  a  firing  pin  strikes  and  causes  the  powder 
train  to  burn.  The  powder  train  at  burning  causes  the  powder  in  the  projectile 
to  explode.  As  can  be  seen  in  Figure  10,  the  X-rayed  noses  of  the  models  prior 
to  firing  show  the  holes  with  the  powder  trains  in  the  balls  to  be  in  various 
positions.  After  these  inert  models  were  fired,  all  of  them  were  recovered 
from  the  cotton  waste  and  X-rayed  again.  The  models  penetrated  the  cotton 
waste  to  a  depth  of  twenty-five  feet.  After  firing  (Figure  ll)  the  holes 
were  lined  up  axially  (except  for  one)  which  would  allow  the  pin  to  arm  the 
fuze. 

PLASTICS  ROTATING  BANDS 

The  Plastics  Branch  of  the  Chemistry  Division  at  NOL  has  been  working  on 
a  plastics  material  to  replace  the  metal  rotating  bands  on  projectiles.  If 
plastics  can  be  used,  gun  wear  can  be.  decreased  and  the  life  of  a  rapid-fire 
gun  can  be  increased.  Models  with  plastics  rotating  hands  were  fired  from 
the  Block  Mount  Range  to  see  which  plastics  could  serve  satisfactorily  as  a 
rotating  band.  Upon  selecting  the  most  promising  plastics,  additional  models 
were  fired  in  the  Aerodynamics  Range  to  insure  that  the  proper  spin  was  being 
attained.  Some  models  were  heated  to  l60  degrees  Fahrenheit  while  other 
models  were  cooled  to  -7 6  degrees  Fahrenheit  to  see  that  they  would  perform 
satisfactorily  at  any  temperature  that  might  be  encountered.  A  maximum  of 
49  seconds  was  used  in  taking  any  one  projectile  from  the  furnace  or  freezer, 
putting  it  in  the  gun  and  firing  it.  The  Microflash  photographs  of  Figure  12 
show  that  both  the  plastics  on  these  models  held  intact  under  the  two  extreme 
temperature  conditions. 
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SPIN  TELEMETERING 


In  order  to  determine  the  spin  rate  of  a  small  projectile,  a  spin  sonde 
■was  developed  in  the  ranges  at  NOL.  larger  sondes  iLte  well-known,  but  this 
one  was  made  to  fit  into  a  20-nan.  shell.  The  construction  and  use  of  this 
spin  sonde  are  given  in  NavOrd  Report  3726.  The  sonde,  shown  in  Figure  13, 
contains  a  condenser,  resistor,  transistor,  and  coil  placed  in  a  plastic 
potting  compound.  A  battery  was  Used  which  will  supply  1.35  "volts  and  which 
will  operate  the  oscillator  for  200  hours.  The  oscillator  coil  serves  as  a 
loop- type  transmitting  antenna  with  the  longitudinal  axis  of  the  projectile 
lying  in  the  plane  of  the  loop  or  oscillator  coil.  The  frequency  of  oscillation 
lies  within  the  range  from.  455  to  465  kilocycles.  A  receiving  antenna  consists 
of  an  open  wire  transmission  line  running  parallel  to  the  projectile  trajectory 
with  the  projectile  trajectory  midway  between  the  wires.  As  the  projectile 
rotates  about  its  axis,  the  transmitting  antenna  emits  a  signal.  At  the 
receiver  this  signal  will  vary  in  Intensity  from  zero  to  some  maximum  value 
and  back  to  zero.  For  each  revolution  of  the  projectile  there  are  two 
positions  of  maximum  signal  and  two  positions  of  zero  signal.  The  transmitted 
signal  is  picked  up  by  the  receiving  antenna  and  is  displayed  on  a  cathode-ray 
oscilloscope.  The  oscilloscope  is  photographed  with  a  high-speed  camera. 

Several  models  of  the  type  shown  In  Figure  l4  were  fired  in  the  Aerodynamics 
Range  at  a  maximum  velocity  of  1200  ft/sec.  These  models  were  accelerated  at 
about  30,000  g*s.  A  sample  measurement  of  the  film  strip  indicated  that  over 
seven  spin  cycles  it  was  possible  to  determine  the  spin  rate  of  1152  EPS  with 
a  standard  deviation  of  5*7  RPS.  Another  sonde  under  development  at  NOL  is 
the  base  pressure  sonde.  A  few  shots  have  been  made.  Although  these  shots 
are  not  conclusive,  the  results  are  encouraging.  Ground  work  is  also  being 
made  toward  the  development  of  a  temperature  telemetering  device . 

HIGH-VELOCITT  POWDER  GUN  FIRINGS 

Improvements  and  modifications  are  constantly  being  made  to  the  ranges 
and  their  instrumentation.  Cylinders  and  spheres  were  recently  fired  in  the 
Aerodynamics  and  Pressurized  Ballistics  Ranges  In  an  attempt  to  determine  the 
highest  velocities  at  which  the  existing  circuitry  of  the  ranges  would  perform 
satisfactorily.  Normal  large  chamber,  smoothbore  powder  guns,  20,  J),  and  40-mm 
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in  diameter,  were  used.  The  models  were  made  of  nylon  and  weighed  from 
about  20  grams  to  50  grams .  The  maxi  mum  velocity  obtained  at  the  muzzle 
in  the  Aerodynamics  Range  at  atmospheric  pressure  for  a  HO-mm  model  was 
about  7000  ft/sec.  This  range  operated  satisfactorily  after  slight  changes 
were  made  in  the  circuits  so  that  the  models  would  trigger  the  sparks  at 
each  station.  The  maximum  velocity  obtained  on  a  30-mm  model  (Figure  15)  in 
the  Pressurized  Ballistics  Range  at  0.07  atmosphere  was  about  98OO  ft/sec 
at  the  muzzle.  No  changes  were  needed  in  the  circuits  of  this  range.  Major 
effort  is  now  being  devoted  to  the  development  of  models  which  will  permit 
1  anrif* tying  at  such  velocities  and  to  the  development  of  a  shorter  duration 
spark. 
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FIG.  2 

RECOVERED  20  MM  MODELS  FIRED  FROM  A  RIFLED  GUN 
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FIG.  4 

SUB  CALIBER  MODEL  FIRED  FROM  40  MM  SABOT 


FOR  LOW 

VELOCITY  FIRINGS 
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FIG.  5 

SABOT  FOR  MULTIPLE  LAUNCHING  OF  SPHERES 
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X-RAY  SAMPLE  OF  30  MM  FUZES  BEFORE  FIRING 
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COTTON  WASTE  IN  PRESSURIZED  YAW  CARD  RANGE. 
THESE  ROUNDS  WERE  FIRED  UNDER  CONDITIONS 
SIMULATING  AN  ALTITUDE  OF  40,000  FEET 
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FIG.12 

MODELS  WITH  PLASTIC  ROTATING  BANDS 
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SPIN  SONDE  AFTER  POTTING 


FIG.  14 
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FIG. 15 

A  NYLON  CYLINDER  FIRED  FROM  A  POWDER  GUN 
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APPENDIX 


Part  I 

The  following  publications  describe  ranges  presently  in  operation 
in  North  America. 

1.  Bull,  G.  V. ,  "Some  Aerodynamic  Studies  in  the  C.A.R.D.E. 
Aeroballistics  Range" ,  Canadian  Aeronautical  Journal, 

Vol.  2,  No.  5,  pp.  154-165,  May  1956. 

2.  May,  Albert  and  Williams,  T.  J. ,  Free-Flight  Ranges  at  the 
Naval  Ordnance  Laboratory,  NAVORD  Report  4063,  July  1955* 

5.  Rogers,  Walter  K.,  Jr.,  The  Transonic  Free-Flight  Range, 

BRL  Report  No.  849,  Feb.  1953. 

4.  Seiff,  Alvin,  A  Free-Flight  Wind  Tunnel  for  Aerodynamic 
Testing  at  Hypersonic  Speeds,  NACA  Report  1222,  May  1955* 

5.  Staff,  Aeroballistics  Laboratory,  Dynamic  Aerobalilstlc 
Evaluation,  NOTS  1152,  July  1955. 
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Part  II 


Additional  Information  on  North  American  range  facilities  can  be 
obtained  by  -writing  to  the  following. 

National  Advisory  Committee  for  Aeronautics 
Ames  Aeronautical  Laboratory 
Moffett  Field,  California 
Attn:  Mr.  H.  Julian  Allen 

Chief  Superintendent 

Canadian  Armament  Research  and  Development  Establishment 
P.  0.  Box  ibZJ,  Quebec,  Province  of  Quebec >  Canada 
Attn:  Dr.  Gerald  V.  Bull 

Commander 

U.  S.  Naval  Ordnance  Test  Station 
China  Lake,  California 

Attn:  Dr.  William  Haseltine,  Code  505 

Commanding  General 

Aberdeen  Proving  Ground,  Maryland 

Attn:  Dr.  Boris  G.  Karpov,  Ballistic  Research  Laboratories 
Commander 

V.  S.  Naval  Ordnance  Laboratory 
White  Oak.,  Silver  Spring,  Maryland 
Attn:  Dr.  Albert  May 
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Superintendent 
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Monterey,  California 

Attn:  Dr.  Head 
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Commander 
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M.  Fitzgerald 
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Chief  of  Naval  Research  2 
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Wright  Air  Development  Ctr. 
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Nellis  Air  Force  Base,  Nevada 
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Air  Force  Armament  Center 
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ACB  -  J.  E.  Stevens 

ACB  -  A.  S.  Galbraith 

Commander 

Air  Research  &  Development 
Command 
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Commander  10 

Air  Force  Missile  Test 
Center  (MTE) 

Patrick  Air  Force  Base,  Florida 

U.  S.  Atomic  Energy  Commission 
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Organization 

National  Advisory  Committee 
for  Aeronautics 
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National  Advisory  Committee 
for  Aeronautics 
Levis  Flight  Propulsion  Lab. 
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Attn:  E .  Reshotko 
F.  K,  Moore 

National  Advisory  Committee 
for  Aeronautics 
Ames  Aeronautical  Lab. 
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Attn:  H.  J.  Allen 
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L.  Neice 
A.  C.  Charters 

National  Advisory  Committee 
for  Aeronautics 
Langltey  Aeronautical  Lab. 

Langley  Field,  Virginia 
Attn:  J.  D.  Bird 
R.  Hopko 
P.  Huber 
A.  Sabol 
C .  E .  Brovn 
A.  Busemann 

Director 

Armed  Services  Technical 
Information  Agency 
Documents  Service  Center 
Knott  Building 
Dayton  2,  Ohio 
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Assistant  Secretary  of  Defense  (R&Dj) 
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Washington  25,  D.  C. 
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Attn:  Samuel  Feltman 

Ammunition  Labs. 

2  Commanding  General 

Frankford  Arsenal 
Bridge  and  Tacony  Streets 
Philadelphia  37#  Pennsylvania 

5  Commanding  General 

Redstone  Arsenal 
Huntsville,  Alabama 
Attn:  Technical  Library 
W.  D.  Murphree 
D.  H.  Newby 
C.  L.  Northrop 
N.  M.  Shapiro 


2  Commanding  Officer  1 

Diamond  Ordnance  Fuze  Labs. 
Building  92 
Washington  25,  D.  C. 

1  Commanding  Officer 

Watertown  Arsenal 
Watertown  72,  Massachusetts  2 

Attn:  R.  Muldoon 

1  Commanding  General 

Arnold  Engineering  Development 
Center 

Tullahoma,  Tennessee 

Attn:  Deputy  Chief  of  Staff,  R&D 

1  Commanding  Officer  and  Director 

David  W.  Taylor  Model  Basin 
Washington  7#  D.  C. 

Attn:  Aerodynamics  Lab. 


Johns  Hopkins  University 
Operations  Research  Office 
Department  of  the  Army 
7100  Connecticut  Avenue 
Chevy  Chase,  Maryland 
Washington  15,  D.  C. 
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Advanced  Development  Corp. 
20  South  Union  Street 
Lawrence,  Massachusetts 
Attn:  T.  R.  Munson 

W.  B.  Stephenson 

THRU:  District  Chief 

Boston  Ordnance  Dist. 
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Armour  Research  Foundation  of  1 
the  Illinois  Institute  of  Tech. 
Technology  Center 
Chicago  1 6,  Illinois 
Attn:  Mr.  W.  Casier 

THRU:  District  Chief 

Chicago  Ordnance  Dist. 

209  W.  Jackson  Blvd. 

Chicago,  Illinois 
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Applied  Physics  Laboratory 
Johns  Hopkins  University 
8621  Georgia  Avenue 
Silver  Spring,  'Maryland 
Attn:  Mr.  G.  L.  Seielstad 
R.  H.  Cramer 
L.  L.  Cronvich 
H.  A.  Wallskog 

THRU:  Naval  Inspector  of  Ord. 
Applied  Physics  Lab. 

Johns  Hopkins  University  3 
8621  Georgia  Avenue 
Silver  Spring,  Maryland 

Aerophysics  Development  Corp. 

P.  0.  Box  657 

Pacific  Palisades,  Calif. 

Attn:  Dr.  W.  Bollay 
D.  Bitonda 

THRU:  District  Chief 

Los  Angeles  Ordnance  Dist.  1 
55  South  Grand  Avenue 
Pasadena,  California 

Adalia  Limited 
1410  Stanley  Street 
Montreal  2,  P.  Q. ,  Canada 
Attn:  J.  B.  Reid 

Canadian  Army  Staff 
2450  Mass.  Ave.,  N,  W. 
Washington  8,  D.  C. 


Organization 

Aircraft  Armaments,  Inc. 
Cockeysville ,  Maryland 
Attni  N.  J.  La  Costa 

THRU:  Deputy  Dist.  Chief 

Baltimore  Regional  Office 
Philadelphia  Ord.  Dist. 
Main  Post  Office 
Baltimore  1,  Maryland 

ARQ,  Inc. 

Gas  Dynamic  Facility 
Tullahoma,  Tennessee 
Attn:  M.  K.  Kingery 
A.  J.  Zazzi 

THRU:  District  Chief 

Cincinnati  Ord.  Dist. 
Swift  Bldg. 

250  E.  9th  Street 
Cincinnati  2,  Ohio 

Bell  Aircraft  Corporation 
P.  0.  Box  1 
Buffalo  5,  New  York 
Attn:  K.  Pearce 

N.  F.  Meullen 
R.  J.  Whalen 

THRU:  Air  Force  Plant  Rep. 

Bell  Aircraft  Corp. 
Niagara  Falls ,  New  York 

Boeing  Airplane  Company 
Bomber  Weapons  Unit  (Mail 
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Physical  Research  Staff 
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Seattle  2k,  Washington 
Attn:  D.  Martin 

THRU:  Air  Force  Plant  Rep. 

Boeing  Airplane  Company 
Seattle,  Washington 
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Boeing  Airplane  Company  1 

Wichita,  Kansas 
Attn:  R.  E.  Wallace 

THRU:  Air  Force  Plant  Rep. 

Boeing  Airplane  Co. 

Wichita,  Kansas 

Boeing  Airplane  Company 
Plant  2 

Seattle  l4,  Washington  3 

Attn:  R.  E.  Bateman 

THRU:  Air  Force  Plant  Rep. 

Boeing  Airplane  Co. 

Seattle,  Washington 

Budd  Company 
Red  Lion  Plant 

Philadelphia  15,  Pennsylvania 

THRU:  Deputy  Diet.  Chief  1 

Philadelphia  Ord.  Di6t. 

128  N.  Broad  Street 
Philadelphia  2,  Pa. 

CONVAIR 

Div.  of  General  Dynamics  Corp. 

San  Diego,  California 
Attn:  C.  W.  Frick 
E.  Katz 

M.  F.  Romig  1 

S.  V.  Starr 

THRU:  Bureau  of  Aero.  Rep. 

CONVAIR 

San  Diego,  California 

Cornell  Aero.  Laboratory,  Inc. 

Wl-55  Genesee  Street 

Buffalo  21,  New  York 

Attn:  Miss  E,  T.  Evans,  Library  1 

THRU:  Bureau  of  Aero.  Rep. 

Cornell  Aero .  Lab . ,  Inc . 

P.  0.  Box  235 
Buffalo  21,  New  York 
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Ottawa,  Ontario,  Canada 
Attn:  J.  L.  Howland 
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J.  E.  Smith 

THRU:  Canadian  Army  Staff 
2I+5O  Mass.  Avenue 
Washington  8,  D.  C. 

Chance-Vought  Aircraft,  Inc. 

P.  0.  Box  5907 
Dallas,  Texas 

THRU:  Bureau  of  Aero.  Rep. 

Naval  Industrial  Reserve 
Plant  Aeronautical 
P.  0.  Box  5907 
Dallas,  Texas 

CONVAIR 

Division  of  General  Dynamics  Corp. 
Pomona  Division 
P.  0.  Box  1011 
Pomona,  California 

THRU:  Naval  Inspector  of  Ord. 

1675  West  5th  Street 
Pomona,  California 

Cornell  University 
Graduate  School  of  Aero.  Engrg. 
Ithaca,  New  York 
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CONVAIR  1 

Div.  of ‘General  DyBSmics  Corp. 

Fort  Worth  1,  Texas 
Attn:  Mr.  L.  W.  Bonnell 

THRU:  Air  Force  Plant  Rep. 

CONVAIR 

Fort  Worth  1,  Texas 

CONVAIR 

Div.  of  General  Dynamics  Corp.  1 

Ord.  Aerophysics  Laboratory 
Daingerfield,  Texas 
Attn:  J.  E.  Arnold 
R.  J.  Yolluz 
K.  L.  Goin 

THRU:  Assistant  Inspector  of 
Naval  Material 
Ord.  Aerophysics  Lab.  1 

Daingerfield,  Texas 

Canadair,  FJ..  Y 
P.  0.  Box  6087 
Montreal,  P.  Q.,  Canada 
Attn:  D,  A.  Jackman 
H.  J.  Luckert 

THRU:  Canadian  Army  Staff 

2450  Mass.  Avenue,  N.  W. 
Washington  8,  D.  C. 
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Canadian  Westlnghouse  Co. 

Electronic  Division 
Longvood  Road 
Hamilton,  Ontario,  Canada 
Attn:  J.  N.  Leavitt 

THRU:  Canadian  Army  Staff 

2^50  Mass.  Ave.,  N.  W. 
Washington  8,  D.  C. 

Chamberlain  Corp. 

Waterloo,  lova 
Attn:  E.  R.  Capon! 


Organization 

Defense  Research  Board 
Room  4755 
"A"  Building 
Ottawa,  Ontario,  Canada 
Attn:  J.  L,  Oat way 

THRU:  Canadian  Army  Staff 

2450  Mass.  Ave.,  N.  W. 
Washington  8,  D.  C. 

Directorate  of  Armament  Dev. 
Army  Headquarters 
Ottawa,  Ontario,  Canada 
Attn:  W.  B.  Snarr 

THRU:  Canadian  Army  Staff 

2450  Mass.  Avenue,  N.  W. 
Washington  8,  D.  C. 

Douglas  Aircraft  Co.,  Inc. 

El  SequndD  Division 
827  Lapham  Street 
El  Sequndo,  California 
Attn:  F.  C.  Newton 

THRU:  Air  Force  Plant  Rep. 

Douglas  Aircraft  Corp. 
3855  Lakewood  Blvd. 

P.  0.  Box  200 
Long  Beach  1,  Calif. 

Douglas  Aircraft  Co.,  Inc. 
Engineering  Department 
3000  Ocean  Park  Blvd. 

Santa  Monica,  California 
Attn:  R.  M.  Wood,  Library 
W.  S.  Cohen 
H.  Klein 

THRU:  District  Chief 

Los  Angeles  Ord.  Dist. 
55  South  Grand  Ave. 
Pasadena  2,  Calif. 
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1  Emerson  Electric  Manufacturing  1 

Company 

8100  W.  Florissant  Avenue 
St.  Louis  21,  Missouri 
Attn:  Mr.  G.  Hauser 

THRU:  St.  Louis  Air  Procurement 
Di strict 

Oklahoma  City  Air 
Materiel  Area 
1114  Market  Street 
St.  Louis,  Missouri 

2 

1  Firestone  Tire  and  Rubber  Co. 

Defense  Research  Division 
1200  Firestone  Parkway 
Akron  17,  Ohio 
Attn:  V.  E.  Lucas 

THRU:  District  Chief 

Cleveland  Ord.  Dist. 
Lincoln  Bldg. 

1567  E.  6th  St. 

Cleveland  l4,  Ohio 

6  General  Electric  Company  3 

3198  Chestnut  Street 
Philadelphia  h,  Pennsylvania 
Attn:  V.  Kebely 
R.  F.  Peck 
J.  Powers 
A.  M.  Smith 
W.  R.  Warren 
Y.  A.  Yoler 

THRU:  Philadelphia  Air 

Procurement  DiBt. 
l4ll  Walnut  Street  2 . 

Philadelphia  2,  Penna. 


General  Electric  Company 

Schenectady,  New  York 

Attn:  Mr.  F.  V.  Johnson,  A&QE 

THRU:  Rochester  Ord.  Dist. 
Branch  Office 
General  Electric  Co. 

Bldg.  23,  Room  232 
Schenectady,  New  York 
Attn:  Lt.  Col. 

Davies 

General  Mills,  Inc. 

Engineering  Research  and  Dev. 

2003  E.  Hennepin  Avenue 
Minneapolis  13,  Minnesota 
Attn:  R.  I.  Sakamaki 

R.  C.  Huntington, 

Mech.  Div. 

THRU:  Milwaukee  Air  Procurement 
District 

770  N.  Plank  in  ton  Ave. 
Milwaukee,  Wisconsin 

Grumman  Aircraft  Engrg.  Corp. 
Engineering  Plant  No.  5 
Betbpage,  New  York 
Attn:  R.  L.  Gustafson 
A.  E.  Munier 
R.  A.  Scheulng 

THRU:  Bureau  of  Aero.  Rep. 

Grumman  Aircraft  Engrg. 
Corp. 

Bethpage,  L»I. ,  New  York 

Hughes  Aircraft  Company 
Aerodynamics  Department 
Florence  Avenue  at  Teal  St. 

Culver  City,  California 
Attn;  I.  Naiman 

W.  L.  Phillips 

THRU:  Air  Force  Plant  Rep .  ,WEAPD 
Hughes  Aircraft  Co. 
Florence  Ave.  at  Teal  St. 
Culver  City,  California 
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1  M.  W.  Kellogg  Company 

Foot  of  Danforth  Avenue 
Jersey  City  3,  New  Jersey 
Attn:  Miss  E.  M.  Hedley 

THRU:  Inspector  of  Naval 
Material 

Naval  Industrial  Reserve 
Shipyard 

Bldg.  15,  Port  Newark 
Newark  5>  New  Jersey 

3  Lockheed  Aircraft  Corp. 

Missile  Systems  Division 
Van  Nuys,  California 
Attn:  E.  Ber shaded 

E.  T.  Cannon 
R .  S .  Swanson 

THRU:  Air  Force  Plant  Rep.  ,WEAD 
Lockheed  Aircraft  Corp. 
Factory  "A",  P.0.  Box  551 
Burbank,  California 

1  Lockheed  Aircraft  Corp. 

Factory  "A" 

P.  0.  Box  551 
Burbank,  Calif. 

Attn:  Mr.  Ed  Baldwin 

THRU:  Air  Force  Plant 
Rep.,  WEAPD 

Lockheed  Aircraft  Corp. 
Factory  "A”,  P.0*  Boat  551 
Burbank,  California 

4  Glenn  L.  Martin  Company 
Baltimore,  Maryland 
Attn:  J.  M.  Bidvell 

M.  L.  Coon 
L.  L.  Jackson 
T.  Reisert 

THRU:  Bureau  of  Aero.  Rep. 

Glenn  L.  Martin  Co. 
Baltimore  3,  Maryland 


a8 


1  Mass.  Institute  of  Tech. 

Naval  Supersonic  Lab.  80-208 
560  Memorial  Drive 
Cambridge,  Massachusetts 
Attn:  F.  H.  Durgin 

THRU:  Inspector  of  Naval 
Material 

Dev.  Contract  Dept. 

Mass.  Inst,  of  Tech. 
Cambridge  39 >  Mass. 

1  Massachusetts  Institute  of  Tech. 
Instrumentation  Laboratory 
Cambridge  39>  Massachusetts 
Attn:  L.E.  Wilkie  (41-203 

THRU:  Inspector  of  Naval 
Material 

Dev.  Contract  Dept. 

Mass.  Inst,  of  Tech. 
Cambridge  3 9>  Mass. 

2  McDonnell  Aircraft  Corporation 
P.  0.  Box  516 

St.  Louis  3>  Missouri 
Attn:  R.  M.  Flesh 
R.  E.  Rohtert 

THRU:  Bureau  of  Aero.  Rep. 

McDonnell  Aircraft  Corp. 

P.  0.  Box  516 

St.  Louis  3,  Missouri 

1  National  Research  Council  of 

Canada 

Defense  Research  Board 
Ottawa,  Ontario,  Canada 

THRU:  Canadian  Army  Staff 
2^50  Mass.  Avenue 
Washington  8,  D.  C. 

1  North  American  Aviation,  Inc, 

Columbus,  Ohio 
Attn:  W.  Simon 

THRU:  Air  Force  Plant  Rep. 

North  American  Aviation, Inc . 
Columbus  Plant 
Columbus,  Ohio 
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4  North  American  Aviation,  Inc.  1 

Missile  Development  Division 
12214  Lakewood  Blvd. 

Downey,  California 
Attn:  E.  Briggs 
J.  K.  Dew 
R.  B.  Oliver 
J.  Elms 

THRU:  Air  Force  Plant  Rep. 

North  American  Aviation,  2 
Inc . 

Los  Angeles  International 
Airport 

Los  Angeles  45,  Calif * 

1  North  American  Aviation,  Inc. 

Los  Angeles  International 
Airport 

Los  Angeles,  California 
Attn:  J.  Covert 

G.  Bussiere  5 

THRU:  Air  Force  Plant  Rep . ,WEAPD 
North  American  Aviation, 
Inc. 

Los  Angeles  International 
Airport 

Los  Angeles,  California 

1  North  Carolina  State  College 

Raleigh,  North  Carolina 
Attn:  Prof.  J.  W.  Cell  1 

THRU:  District  Chief 

Philadelphia  Ord.  Dist. 

128  North  Bg$>ad  St. 
Philadelphia  2,  Pa. 

1  Northrop  Aircraft,  Inc. 

Dept.  5485 

Ogden  Air  Material  Area 
Hawthorne,  California  2 

Attn:  D.  C.  Olmore 

THRU:  Air  Force  Plant  Rep./WEAFD 
Northrop  Aircraft,  Inc. 
Hawthorne,  California 
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Princeton  University 

Head,  Gas  ’Dynamics  Laboratory 

James  Forres tal  Research  Center 

Princeton,  New  Jersey 

Attn:  Prof.  S.  M.  Bogdonoff 

THRU:  Office  of  Naval  Research 
246  Broadway 
New  York  15,  New  York 

The  Ramo-Woolrtdge  Corporation 
5750  Arbor  Vitae  Street 
Los  Angeles  45,  California 
Attn:  A.  Ambrosio 
W.  D.  Hayes 

THRU:  District  Chief 

Los  Angeles  Ord.  Dist. 

55  S.  Grand  Avenue 
Pasadena,  California 

Republic  Aviation  Corporation 
Scientific  Staff 
Farming  dale ,  New  York 
Attn:  W.  Mcllroy 

C.  Renneinann 
R.  Senator 

THRU:  Air  Force  Officer -In-Charge 
Republic  Aviation  Corp. 
Farmingdale,L.I. ,  New  York 

Republic  Aviation  Corporation 
Guided  Missiles  Division 
255  Jericho  Turnpike 
Mineola,  New  York 
Attn:  R.  M.  Kennedy 

THRU:  Air  Force  Officer-In-Charge 
Republic  Aviation  Corp. 
Farmingdale,L.I. ,  New  York 

A.  V.  Roe  Canada,  Ltd. 

Aircraft  Division 
Malton,  Ontario,  Canada 
Attn:  J.  A.  Chamberlin 
W.  Taylor 

THRU:  Canadian  Army  Staff 

2450  Mass.  Ave.,  N.  W. 
Washington  8,  D.  C. 
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1  Sperry  Gyroscope  Company 

Div.  of  The  Sperry  Corp. 

Great  Keck,  L.I.,  Nev  York 
Attn:  J.  J.  Gallaghan, 

Federal  Dept. 

THRU:  Assistant  Inspector 
of  Naval  Material 
c/o  Sperry  Gyroscope  Co 
Div.  of  Sperry  Corp. 

Great  Neck,  L.I.,  N.  Y. 

1  University  of  Michigan 

Aero.  Research  Center 
Willow  Run  Airport 
Ypsilanti,  Michigan 
Attn:  J.  E.  Corey 

THRU:  Commander 

Central  Air  Proc.  Diet.  1 
W.  Warren  &  Lonyo  Aves. 
Detroit  32,  Michigan 

5  University  of  Minnesota 

Rosemount  Aero.  Laboratories 
Rosemount,  Minnesota 
Attn  i  R.  V.  Deleo 
R.  Herman 
F.  A.  Moynihan 

THRU:  Milwaukee  Air  Proc.  Dist. 
770  N.  Plank lnton  Ave. 
Milwaukee,  Wisconsin 

1  University  of  Southern  Calif. 

Engineering  Center 
5518  University  Avenue 
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